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Christopher Michael Harris 

Abstract 

The Standard Model of particle physics has been remarkably successful in describing 
present experimental results. However, it is assumed to be only a low-energy effective 
theory which will break down at higher energy scales, theoretically motivated to be 
around 1 TeV. There are a variety of proposed models of new physics beyond the 
Standard Model, most notably supersymmetric and extra dimension models. 

New charged and neutral heavy leptons are a feature of a number of theories of 
new physics, including the 'intermediate scale' class of supersymmetric models. Using 
a time-of-flight technique to detect the charged leptons at the Large Hadron Collider, 
the discovery range (in the particular scenario studied in the first part of this thesis) is 
found to extend up to masses of 950 GeV. 

Extra dimension models, particularly those with large extra dimensions, allow the 
possible experimental production of black holes. The remainder of the thesis describes 
some theoretical results and computational tools necessary to model the production and 
decay of these miniature black holes at future particle colliders. The grey-body factors 
which describe the Hawking radiation emitted by higher-dimensional black holes are 
calculated numerically for the first time and then incorporated in a Monte Carlo black 
hole event generator; this can be used to model black hole production and decay at next- 
generation colliders. It is hoped that this generator will allow more detailed examination 
of black hole signatures and help to devise a method for extracting the number of extra 
dimensions present in nature. 
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Preface 



This constitutes my research for a PhD degree in theoretical high-energy physics. Chap- 
ter ^ contains a review of the relevant theoretical background. It sets the context in 
which the work described in the rest of the thesis was carried out, and introduces the 
theoretical models on which the work is based. Chapter |21 describes work undertaken 
during the first year of my PhD on the prospects for detecting exotic heavy leptons at 
the Large Hadron Collider (LHC). This was published [1] as a study performed by mem- 
bers of the Cambridge Supersymmetry Working Group, but the contents of Chapter |21 
are entirely my own work. 

The remainder of the thesis concentrates on the study of miniature black holes at the 
LHC — a topic which is introduced in Chapter El Chapter |3 describes new theoretical 
work on black hole grey-body factors (building on previous work by Kanti and March- 
Russell); much of this work was published in [2] and contributions from my collaborator, 
P. Kanti, are clearly identified in Chapter 01 The final chapter gives details of the 
CHARYBDIS event generator written to simulate black hole production and decay. This 
was published in [3] and, as detailed in ChapterEl involved some significant contributions 
from P. Richardson (these were mainly involved in interfacing CHARYBDIS to existing 
Monte Carlo programs). The chapter also contains a brief discussion of some preliminary 
experimental studies using the event generator. Some of this work will be published in [4], 
together with more detailed analyses by members of the ATLAS collaboration. 
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Chapter 1 



Introduction 



1.1 The Standard Model 

The Standard Model (SM) of particle physics has been remarkably successful. Many 
of the parameters are measured to great accuracy, and the Standard Model has made 
many predictions which have been verified experimentally. However there are several 
reasons why it is widely believed that, whilst working well in the energy regimes which 
have been investigated to date, the Standard Model does not give the full picture. These 
issues, which are both technical and aesthetic, will be discussed further in the following 
sections. It is the belief that there is physics 'Beyond the Standard Model' (BSM) which 
motivates the continuation of large scale particle physics experiments, most notably the 
Large Hadron CoUider (LHC) at CERN which will collide protons with a centre-of-mass 
energy of 14 TeV. Of course, the other great aim of the LHC is to discover the famous 
Higgs boson, the only missing particle of the Standard Model. Although the Higgs 
mechanism is as yet unproved, much of the work in this thesis makes the implicit and 
sometimes explicit assumption that a Standard Model Higgs boson exists. 
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1.1 The Standard Model 



1.1.1 Particle content 

The particle content of the Standard Model is intimately related to the symmetries 
of the Lagrangian which describes the physics. The forces between the fermions are 
mediated by gauge boson fields which allow the Lagrangian to be invariant under gauge 
transformations — it is believed that all fundamental interactions are described by some 
form of gauge theory. These symmetries imply the existence of conserved currents and 
charges in the theory. 

Firstly the situation before the breaking of electroweak symmetry is reviewed. At 
this stage all particles are massless as the Higgs mechanism has not yet been invoked to 
give masses to the particles. 

There are three main symmetries of the Standard Model Lagrangian. Firstly there 
is the U(l) hypercharge symmetry; this is an Abelian symmetry and the corresponding 
gauge boson field is usually written as B^. The conserved charge, the hypercharge, is 
most commonly called Y. 

There is also an SU(2) weak isospin symmetry which is sometimes referred to as 
SU(2)l since the gauge field couples only to left-handed fermions. This is a non-Abelian 
symmetry since the Lie group generators are matrices. There are three accompanying 
gauge fields: VFf , W2 and and the conserved charge is I3, the third component of 
the weak isospin /. 

Finally there is the SU(3) non-Abelian symmetry of Quantum Chromodynamics 
(QCD) which has a conserved colour charge. The gauge bosons in this case are glu- 
ons and are described by the fields A^^ where the colour index a runs from 1 to 8. The 
gluons themselves carry colour charge — red, green and blue. 

The theory outlined above is clearly incomplete since we know that most observed 
particles are not massless. It is the introduction of a Higgs field (a complex scalar field 
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with a non-zero vacuum expectation value) to the theory which both gives most of the 
particles mass and causes the electroweak symmetry breaking. The from the SU(2) 
of isospin mixes with the B^^ from the U(l) of hypercharge to produce a massless photon 
field and a massive Z boson. At the same time the Wi and W2 states which can be 
identified with the bosons acquire mass, as do all the fermions with the exception 
of the neutrinos.'^ After this mixing there is a residual U(1)q electromagnetic gauge 
symmetry; the coupling of fermions to the photon is found to be proportional to the 
electromagnetic charge Q, which is related to the hypercharge and weak isospin through 
Q = Y + I^.^ The coupling to the Z boson is more complicated but is proportional to 
h ~ Q sin^ 6'w where is the Weinberg angle and sin^ 6'w ~ 0.23. Since I3 is different 
for the left and right-handed fermions, the coupling to the Z boson is also dependent on 
the fermion helicity. 

In Table im the fermionic particle content of the Standard Model is summarized and 
the various quantum numbers of the different particles are shown. From this table, in 
which L and R refer to left- and right-handed states respectively, it is clear why the 
isospin symmetry is often referred to as SU(2)l. 

Apart from the fermions (with spin s = 1/2) the SM particle content consists of 
the gauge bosons [s = 1) and the Higgs boson (s = 0). The Higgs field introduced in 
the Higgs mechanism has four degrees of freedom, but three of these are 'eaten' by the 
isospin gauge bosons to give mass to the W and Z bosons. Therefore after the symmetry 
breaking there is a massive neutral scalar with one degree of freedom — this is the particle 
known as the Higgs boson. Although the Higgs boson has not yet been discovered and 
the significance of hints during the final running of the Large Electron-Positron (LEP) 
collider has been reduced by more careful analyses [5], a low Higgs mass (~ 115 GeV) is 

''In the minimal Standard Model neutrinos are massless although recent evidence suggests that this 
is not the case. 

''It is not uncommon for T to be used for the isospin rather than /, or for hypercharge to be defined 
such that this relationship becomes Q = Y/2 + I3 . 



6 



1.1 The Standard Model 
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Table 1.1: Fermionic particle content of the Standard Model showing the rele- 
vant quantum numbers. 

still favoured. This is because of tight constraints coming from precision measurements 
of electroweak radiative corrections. 

In addition to the Standard Model gauge symmetries discussed above, there are two 
further 'accidental symmetries' in that all the terms in the Lagrangian are found to 
conserve both lepton number L and baryon number B. Lepton number is defined to be 
+1 for leptons, —1 for anti-leptons and zero for all other particles; baryon number is 
defined to be +1/3 for quarks, —1/3 for anti-quarks and zero for all other particles (this 
means that baryons have B = \). These symmetries were not imposed on the model, 
but appear naturally given the requirement that the terms in the Lagrangian should be 
renormalizable, Lorentz invariant and gauge invariant. In fact B and L are not exactly 
conserved because they are violated by non-perturbative electroweak effects [6, 7] , but 
they can normally be considered as conserved quantum numbers. Unfortunately for 
most theories of physics beyond the Standard Model the most general Lagrangian which 
can be written down includes terms which violate these symmetries. Potentially this 
causes various problems, most notably that it can allow protons to decay with a much 
shorter lifetime than the present experimental limit. 
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1.2 Motivations for new physics 

It has already been mentioned that there are various reasons for beheving there is physics 
beyond the Standard Model. The work in this thesis will address some, but by no means 
all, of these issues which are detailed below. 

1.2.1 Gravity 

A full model for particle physics would be expected to describe all the fundamental 
forces between particles — electroweak, strong and gravitational. However the Standard 
Model fails to describe any details of the gravitational force, and hence cannot be the full 
fundamental model. At energy scales of order the Planck mass, Mp, a theory of quantum 
gravitation will be required to describe the interactions between particles. This shows 
that the Standard Model will need to be replaced by an alternative theory at very 
high energy scales, and it is reasonable to expect that a more complete model than the 
Standard Model might be required even at energy scales only moderately higher than 
those which have already been investigated in detail. 

1.2.2 Hierarchy problem 

One of the most common arguments for physics beyond the Standard Model is what is 
known as the hierarchy problem. This can be expressed in two different ways which, 
whilst clearly related, are in some ways distinct. They are referred to here as the 
'aesthetic hierarchy problem' and the 'technical hierarchy problem'. 
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1.2 Motivations for new physics 



1.2.2.1 Aesthetic hierarchy problem 

As has been discussed above, the Standard Model has nothing to say about gravitational 
interactions. It also seems aesthetically unappealing that the mass scale with which 
particle physicists are most familiar (the electroweak mass scale, ~ 100 GeV) is so many 
orders of magnitude below the Planck scale (~ 10^^ GeV). If there is a fundamental 
'theory of everything' then it would seem more natural if all the energy scales were of 
the same order of magnitude. Another way of expressing this is just to say that it seems 
unnatural for gravity to be so much weaker than the other gauge forces. Why are there 
seemingly two different fundamental scales in nature? 

1.2.2.2 Technical hierarchy problem 

More technically, the hierarchy problem is the result of considering one-loop corrections 
to the Higgs mass. It is relatively easy to calculate the correction to the Higgs mass from 
diagrams like that in Figure fTTT| but since the diagram is quadratically divergent, a cut- 
off (A) must be introduced in the loop momentum integral to regulate the divergence. 
The correction to the 'bare' Higgs mass is 

^ 1^ t"^"^' + ^""f ln(A/mf) + ■■■]. (1.1) 

where Af is the Higgs-fermion coupling and mf the fermion mass. In itself this mass 
correction isn't necessarily a problem since the bare Higgs mass can have any value 
provided the physical Higgs mass is around the electroweak scale. The Standard Model 
is only considered to be a low-energy effective theory and so a natural value for A would 
be the Planck scale. However this means that the 'bare Higgs mass' would have a 
natural value of order the Planck scale and therefore, if there is a Higgs boson with 
mass of order 100 GeV, there would need to be cancellation at the level of 1 part in 
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f 
f 

Figure 1.1: Feynman diagram for the fermion one- loop correction to tlie Higgs 
boson mass. 

10^^. This is considered to require enormous 'fine-tuning' to the parameters of the bare 
Lagrangian. 

Assuming we beheve in the existence of the Standard Model Higgs boson and find the 
idea of fine-tuning unnatural this leave only one possibility: there must be new physics 
and/or new particle content at energy scales A ~ 1 TeV. In fact even if there is no 
SM Higgs boson, new TeV-scale physics is still required in order to unitarize the W-W 
scattering cross section. 

1.2.3 Unification of couplings 

If the strong and electroweak forces are to unify into a single gauge theory at a high 
energy scale then the gauge couplings must also unify. The particle content of a model 
determines how the couplings 'run' with energy and although in the Standard Model they 
almost unify at a scale of about 10^^ GeV, this unification isn't quite exact (as illustrated 
by Figure ITT^ . This can be seen as a motivation for additional particle content or new 
physics which will modify the running of the couplings to make the unification more 
exact. 
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1.3 Supersymmetry and intermediate scale models 
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Figure 1.2: Running of the gauge couplings assuming only Standard Model 
particle content. Values of the inverse couplings at energy /i are determined 
using the one-loop renormalization group equations, ai, a2 and as are the 
hypercharge, weak isospin and strong couplings respectively, but the conventional 
normalization means that ai = (5/3)ay. The width of the bands on the plot is 
due to experimental uncertainties. (Reproduced from [8]). 

1.3 Supersymmetry and intermediate scale models 

One type of new physics which is relevant to the work presented in this thesis is su- 
persymmetry (SUSY). This is one of the most popular choices for physics beyond the 
Standard Model and resolves several of the issues raised in the previous section. The 
basic idea is that an additional symmetry is introduced to the theory such that the La- 
grangian is invariant under transformations which convert fermions into bosons and vice 
versa. This immediately doubles the number of particles and more importantly helps 
to solve the technical hierarchy problem. This is because there will be cancellations 
between Feynman diagrams like that in Figure 11.11 and one-loop diagrams with S (the 
fermion's supersymmetric scalar partner, known as a sfermion) as shown in Figure IT751 
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go go 

Figure 1.3: Feynman diagram for the scalar one-loop correction to the Higgs 
boson mass. 



The Higgs mass correction in this case is 

'^^^s = ^""^ ln(A/ms) + ■■■], (1.2) 

with As the Higgs-scalar couphng and ms the scalar mass. Hence provided there are two 
scalar particles (with coupling As = |Afp) for every fermion, there will be cancellation 
between the fermion and sfermion loop diagrams. In fact the cancellation is only between 
the terms quadratic in A but the remaining Higgs mass corrections are only logarithmic 
in A and hence do not cause the same fine-tuning problems. Although we know that, 
if it exists, supersymmetry must be broken (since no superpartners of Standard Model 
particles have yet been observed) this supersymmetry breaking must occur in such a way 
that the leading-order cancellation remains. The requirement for the physical Higgs mass 
to be at the electroweak scale is now that the superparticles have masses of order 1 TeV 
or less. If this is not the case then fine-tuning problems start to resurface and one of 
the main motivations for SUSY has been lost. It is also difficult to construct SUSY 
models where the SM particle masses differ from those of their superpartners by more 
than about an order of magnitude. 

Various different ways to break supersymmetry have been studied. It has long been 
known that a viable way to break supersymmetry is in a hidden sector — the SUSY 
breaking in a particular model is then transmitted to the 'outside world' by some means 
(usually gravitational or gauge interactions). The supersymmetry breaking means that 
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1.3 Supersymmetry and intermediate scale models 



what is initially a very beautiful theory (and hence very attractive to particle theorists) 
becomes much less beautiful. Many new parameters (~103) are needed to describe a 
particular supersymmetric model, and only by making simplifying assumptions (e.g. 
regarding various high-scale masses and couplings) can much headway be made in con- 
sidering the phenomenology of these models. This makes supersymmetry much less 
attractive, particularly to SUSY sceptics who argue that a region of parameter space 
can always be found to ensure that supersymmetry cannot be ruled out by either exper- 
imental or cosmological constraints. 

However there are several other arguments in favour of supersymmetry. It can provide 
a dark matter candidate in the form of a lightest supersymmetric particle (LSP). This is 
the result of the R-parity symmetry which is often imposed to prevent proton decay but 
which also has the effect of making the LSP stable since it is forbidden from decaying 
into only SM particles. One of the other motivations for supersymmetry is that it 
improves the unification of gauge couplings discussed in ^1.2.31 (although the present 
value of as(^z) seems to be slightly too low to allow perfect unification — see [9] for 
a recent discussion). Figure ED shows the running of the gauge couplings with energy 
in the supersymmetric case. However the couplings do not unify at ~ 10^^ TeV in 
all supersymmetric scenarios. One scenario, motivated by string theory, is that of an 
intermediate scale [10-13]. In this type of model, gauge coupling unification is assumed 
to occur at an 'intermediate' energy scale (~ 10^^ GeV) which is the geometric mean 
of the electroweak and Planck energy scales. This choice is motivated by the scale 
of supersymmetry breaking in hidden-sector gauge-mediated scenarios and also by the 
usual solution of the strong Charge Parity (CP) problem. 

The strong CP problem is that there is an unnaturally small bound (< 10^^) on the 
B parameter of the QCD Lagrangian. The parameter is constrained to be this small 
because the B-term in the Lagrangian is CP-violating and would otherwise lead to a 
neutron electric dipole moment (on which there are strong experimental limits). The 
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Figure 1.4: Running of the gauge couplings assuming Minimal Supersymmetric 
Standard Model particle content. For more details of the quantities plotted, see 
Figure 11.21 (Reproduced from [8] ) . 

Peccei-Quinn mechanism [14] solves this problem by introducing an additional global 
chiral symmetry, U(1)pq, which is then spontaneously broken generating a Goldstone 
boson known as the axion field. The axion acquires a small mass due to non-perturbative 
effects and and has a lifetime which exceeds the age of the Universe by many orders 
of magnitude. The requirement that axion emission does not over-cool stars puts a 
lower bound bound on the axion decay constant /a, and an upper bound comes from 
dark matter constraints. Therefore /a is constrained to be in the approximate energy 
range 10^-10^^ GeV, as is naturally the case for scalars in intermediate scale string 
scenarios [11]. 

The intermediate scale is also consistent with some neutrino mass mechanisms, and is 
favoured by certain cosmological inflation models. High-energy cosmic ray observations 
and certain non-thermal dark matter candidates are also supporting arguments for an 
intermediate energy scale. For more details, see [10,11] and the references contained. 
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1.4 Extra dimension models 



Various intermediate scale models are motivated from different string models in ref- 
erences [10-13]. Some phenomenological constraints on intermediate scale models are 
provided by [15], and constraints from the muon anomalous magnetic moment have also 
been discussed [16,17]. One of the ways of achieving this intermediate scale unification 
is to include new particles (with ~ TeV masses) as part of extra supermultiplets added 
to the Minimal Supersymmetric Standard Model (MSSM). Extra heavy leptons added 
in such a scenario are the subject of the work described in Chapter El 

1.4 Extra dimension models 

Extra dimension models are another class of model which have recently become quite 
fashionable and also provide a way of solving the hierarchy problem.'^ The essence of 
these models is that the apparent weakness of gravity at macroscopic length scales is 
due to the presence of extra dimensions. 

All these models involve one or more (3 + l)-dimensional slices of space-time (known 
as branes or 3-branes) which are embedded in the full extra-dimensional space-time (the 
bulk). In most models the world as we know it and all the Standard Model particles are 
localized on the brane whilst gravitons can propagate in the bulk. There are, though, 
some models where SM particles are not localized on the brane — one will be introduced 
below. Another feature of extra dimension models is that the new dimensions are com- 
pactified or 'curled up' — necessary to explain why they have never been observed in the 
same way as the infinite dimensions we are familiar with. 

Different models involve extra dimensions of different sizes but it is initially sur- 
prising to learn that they can be relatively large (~1 mm) without being ruled out 

■^Many would argue, with good reason, that extra dimensions reaUy just reformulate the problem 
rather than solve it. 
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by experimental measurements of gravity. Extra dimensions of this size imply that the 
Newtonian force law will fail at these length scales. However the force law has, even 
now, only been verified down to distances of ~ 100 fim (see ^1.5.3. 1|) and when these 
extra dimension models were first proposed in 1998, the limit was only about 1 mm. 

As well as providing a solution to the hierarchy problem, extra dimension models 
can be motivated by string theory. When combined with supersymmetry, string theory 
is a candidate for a theory of quantum gravity but it requires the world to be 10- or 11- 
dimensional. String theory can also naturally account for some fields (e.g. SM particle 
fields) being fixed on (3 + l)-dimensional surfaces or branes whilst others, like gravitons, 
propagate in the extra dimensions. In both type-I and type-II string theories, matter 
and the Standard Model interactions are represented by open strings which start and 
end on a brane, whereas as gravity is described by a closed string which can propagate in 
all the spatial dimensions. It is possible that supersymmetry plays an important role in 
stabilizing the extra-dimensional geometry. Unfortunately, in common with the majority 
of proposals for BSM physics, one prediction of most of these theories is a very large 
value for the vacuum energy density or cosmological constant; this is in contradiction 
with the very small value calculated by cosmologists. 

In fact extra dimension models are not a new idea. As early as the 1920s a model with 
one extra compactified dimension was suggested by Kaluza and Klein in an attempt to 
unify electromagnetism and gravity. In such a scheme the photon is actually a component 
of the higher-dimensional graviton. Although their theory fell out of fashion due to its 
failure to explain the weak and strong forces or the relative strength of gravity, their 
ideas seem to have come back into fashion. 

The specific geometry of the extra dimensions varies in different models. The work 
in this thesis considers phenomenology in large extra dimension models, but two other 
extra dimension scenarios are briefly discussed in this introductory chapter. 
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1.4.1 Introducing different models and their phenomenology 

Before introducing the different extra dimension models we consider the wave equation 
for a scalar field in 5 dimensions; this is useful for demonstrating some of the main 
differences in phenomenology. We consider the fifth dimension to be compactified on a 
circle of radius R such that there is a periodic identification = + 2t:R. The scalar 
field equation for a particle of mass m propagating in the bulk is 

— -V^ + mM $(x,y,^,X5,t)=0, (1.3) 



where 



Qy2 Q^2 

The solution will be of the form 



^{x, y, z, X5, t) = <^n{x, y, z, t) exp ( j , (1.5) 



which gives the following equation of motion: 



— -Vl + m' + ^\<5,ix,y,z,t). (1.6) 

Viewed from a four-dimensional perspective there is clearly a tower of states of different 
masses, with mass splitting l/i?^; this is the so-called Kaluza- Klein (KK) tower of 
states. It is clear though that the Kaluza-Klein tower is simply a result of viewing a 
5-dimensional model in 4-dimensions. Taking a 5-dimensional view there is only one 
state of mass m which in the equivalent spin-2 case will be the massless graviton. 

At this point it is worth pointing out another feature of most of these models: there is 
no requirement for momentum conservation in the extra-dimensions. Energy-momentum 
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conservation in the four- dimensional world is a consequence of the translational invari- 
ance of the four space-time dimensions. However the 3-brane breaks translational in- 
variance in the extra dimension and hence momentum is these directions need not be 
conserved in interactions between brane and bulk states. In effect the brane will be 
considered as being infinitely heavy, or as being 'nailed' in place. 

1.4.1.1 'Large' extra dimensions 

The earliest extra dimensions models (at least in recent years) were those of Arkani- 
Hamed, Dimopoulos and Dvali (known as 'large' extra dimension,'' or ADD models 
[19-21]). These are perhaps the simplest of the extra dimension models and are closest 
to the ideas of Kaluza and Klein. 

Gravity is the only force which propagates in the full volume of the extra dimensions 
(the bulk). Hence the gravitational force in the four- dimensional world (on the brane) is 
'diluted' and gravity appears weak compared to the other forces which do not propagate 
in the extra dimensions. In these models there is only one fundamental energy scale, the 
electroweak scale, and so there is no longer a fine tuning problem — A is now of order a 
TeV. 

For large extra dimensions, equation ()1.6|1 shows that the mass splitting will be small, 
and all graviton modes (up to the mass scale involved in the process) will contribute 
to the process being considered. Although each mode has a negligible contribution, the 
number of gravitons involved can be very large which means care must be taken that 
the total contribution isn't inconsistent with precision SM measurements. 

'^There are models in which this approximation is not vahd, e.g. [18]. 

°The word 'large' may seem to be something of a misnomer, but the extra dimensions are large 
relative to the fundamental energy scale in the models. 
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There are two major classes of experimental signature: direct production of the KK 
excitations of gravitons (giving jets and missing energy) and virtual exchange of the 
KK excitations (modifying SM cross sections and asymmetries). The missing energy 
signatures are clearly different to those in other models with invisible particles (like 
SUSY) since from a four- dimensional point of view the graviton doesn't have a fixed 
mass — instead it appears to have a continuous mass distribution. 

The work in the later chapters of this thesis is almost entirely in the framework of 
ADD models; as a result there is much further discussion of the details of these models 
in grSI 

1.4.1.2 TeV~^ extra dimensions 

The models with extra dimensions of size TeV~^ [22] are somewhat different to other 
models like those of ADD. These models do not themselves attempt to solve the hi- 
erarchy problem; the smaller size of the extra dimensions means that the volume of 
the bulk is not be large enough to explain the weakness of gravity (although there are 
models which include a mixture of large and TeV^^ extra dimensions). In fact the extra 
dimensions are small enough to allow the gauge and Higgs bosons to propagate in the 
bulk without running into conflict with the constraints from the precision electroweak 
measurements made at LEP (for example the Z- width). 

Although they don't solve the hierarchy problem, there are various advantages of 
models like this which have the gauge and Higgs bosons propagating in the bulk. Gauge 
bosons in the bulk provide a mechanism for early unification of the coupling constants 
because the logarithmic running is modified by the extra dimensions [23]. If the Higgs 
field propagates in the bulk, the vacuum expectation value of the Higgs zero mode is 
able to generate spontaneous symmetry breaking. Models with two branes with a TeV~^ 
separation can also suppress proton decay if the quarks are localized on one brane and 
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the leptons on the other [24] (ahernatively there may be a 'thick' brane with quarks and 
leptons localized at different points but with the Higgs and gauge bosons propagating 
inside it [25]). Such models can also be useful to motivate the three generations of 
fermions and the hierarchy of fermion masses. 

In these models, R ~ TeV~^ and so, from equation ()1.6|) . the mass splittings are 
TeV-scale. Precision electroweak data requires that >2-5 TeV, but the LHC should 
be able to probe M^. up to around 14 TeV by searching for leptonic signatures of KK 
excitations of the photon and the Z boson. 

A subset of these models which are slightly different are those of universal extra 
dimensions where all fields are allowed to propagate in the bulk [26]. These can be 
constructed to evade current electroweak limits and might provide a dark matter can- 
didate [27]. One particular feature of universal extra dimension models is that KK 
excitations are pair-produced; this could lead to some distinctive phenomenology at 
collider experiments. 

1.4.1.3 'Warped' extra dimensions 

The 'warped' extra dimension models of Randall and Sundrum (RS) [28,29] solve the 
hierarchy problem in a different way to the ADD models. Gravity is not 'diluted' by 
the large volume of the bulk but by the strong curvature of the extra dimensions. RS 
models involve a second brane; gravity is strong on that brane but a 'warp factor' in the 
metric ensures that gravity on our brane is weak. These models may provide a way to 
explain why the cosmological constant seems to be so small. 

Equation ()1.6|1 cannot strictly be applied in RS models (since the geometry is non- 
factorizable) but the extra dimensions are again small which corresponds to a large mass 
splitting. This gives resonant graviton production where usually only the first graviton 
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resonance is considered since the LHC is unlikely to be of high enough energy to study 
higher resonances. Channels like G — > e+e^ and G — > 77 can be studied to find the 
mass of the first graviton resonance and to confirm its spin-2 nature [30]. 



1.5 Details of ADD models 
1.5.1 Gauss' law in extra dimensions 

Whilst it is relatively trivial to use Gauss' Law to derive the power dependence of the 
gravitational force in extra dimensions, slightly more care is required to determine con- 
stants of proportionality in the force law and to determine how these relate to the Planck 
scale. The difficulties arise both because there are a number of different conventions in 
use and because it is not trivial to visualize geometries involving more than three spatial 
dimensions. In all that follows in this thesis, n toroidally compactified extra dimensions, 
all of radius R, are assumed.^ 

In extra dimensions. Gauss' Law can be written as 

F • dS = f7„+2G(4+n)M, (1.7) 



where F is the force per unit mass, G(4+n) is the (4 + n)-dimensional Newton constant,^ 
M is the mass enclosed by the surface S and f2p (the area of a unit p-sphere) is given 

by 

Q„= . (1.8) 

P r(2±i) ^ ^ 



^Much of the extra dimension phenomenology discussed in this and later chapters could be modified 
by compactification on more complicated manifolds. 

^This is by no means unique as a way of defining the extra-dimensional Newton constant. 
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At short distances (<^ R) the extra compactified dimensions are effectively flat and so 
we obtain 



^■»«M = ^^. (1-9) 

whereas at long distances the usual 4-dimensional force law is obtained. 

The above suggests that G(4+„) ~ G(4)i?" in order to get a smooth matching between 
the two regimes. However it is worth taking a little more care in deriving the correct 
relationship between these constants. We follow the approach of ADD in one of their 
early papers [20] and use the analogy of a line of point masses. First consider an infinite 
two-dimensional line of masses, each of mass m and separated by a distance L. Provided 
we are at a large enough distance r from the line, it will appear like a continuous line of 
mass and the problem has become two dimensional. 

Using Gauss's law in three (spatial) dimensions and considering the flux through a 
cylinder of length d we obtain 

{2nrd)F = 47rG'(4)m i-j — > F{r) = ^^^(4)— . (1.10) 

So there is an effective three-dimensional Newton constant (applicable for two spatial 
dimensions) related to ^(4) by 

G(3) = |G(4) , (1.11) 



which is more usefully written as 



G(3) = ^^G(4) . (1.12) 



The situation with compactified extra dimensions is exactly analogous. The periodic 
identification means that the equivalent uncompactified theory has lines of mirror masses 
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m with separation L = 2ttR. So for r ^ L there is an effective 4-dimensional Newton 
constant related to G'(4+„) by 



1.5.2 Planck mass in extra dimensions 

The Planck mass is the scale at which gravity becomes strong, but there are varying 
definitions even in four dimensions. Initially we follow the conventions of ADD who use 
a reduced Planck scale, M^^^-^ given by 

MuXl^ = r. \ > (1-14) 



(4+n) G'(4+„)fi„+2 



and hence, using eq. ()1.13|1 . 



Mj) = (2vri?)"MS^„) = i?"Mp"+/^„) , (1.15) 

where Mp(4+„) is the Planck mass in 4 + n dimensions. This is introduced because it is 
found to be this scale which is most directly constrained by experimental bounds and 
which is required to be ~ 1 TeV in order to solve the hierarchy problem. However 
in many ways the reduced Planck mass is more natural since the interactions of the 



canonically normalized field in the Lagrangian are suppressed by y ^(4+n) • 

We can always view interactions in two different ways. From a 4-dimensional view- 



point the graviton coupling is very small (~ 1/y M(^4)) but many Kaluza-Klein modes 
contribute (for a process with energy scale E, there will be contributions from ~ {ER)"' 
modes). However from a (4 + n) -dimensional view-point the graviton coupling is much 



larger (~ 1/W M(4+„)) but there is no tower of graviton states. 
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Definition of M"+i , 


Comments 


a 
b 
c 
d 
e 


1 1 

47rG(4) R" 

1 1 

87rG(4) R" 

1 1 

G(4) R" 

1 1 
G(4) (27riJ)" 

1 1 
87rG(4) (27r_R)" 


Used by ADD and also GT [31] 

Used by GRW [32]; often known as M, 

Used by EOT- WASH [33] 

Used by Dimopoulos and Landsberg [34] 

Used by Han, Kribs and McElrath [35] 



Table 1.2: Different conventions for tiie fundamental Planck scale Mp(4+„). 
The last column indicates which conventions are used by the authors of some of 
the papers referenced in this thesis. 



As mentioned above there are a number of other ways of defining the Planck scale, 
and at the same time a number of different definitions of the extra-dimensional Newton 
constant. To compare most transparently with other conventions (which in part differ 
because of alternative definitions of the 4-dimensional Planck mass) equation p.l5|) is 
more conveniently written as 

For ease of reference in later parts of this thesis, a summary of many of the popular 
conventions is shown in Table IT!^ fin all cases n toroidally compactified extra dimensions, 
all of radius R, are assumed). 

1.5.3 Constraints on extra dimension models 

The differences between the conventions (particularly the first three) are usually rela- 
tively small since the conversion factors are close to unity for most allowed values of n. 
When quoting astrophysical and cosmological limits where there are large uncertainties 
it is not crucial to be careful about which convention is being used. However for collider 
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limits the convention should be clearly stated. For ease of comparison all the limits 
below are in the convention (convention 'b' in Table 11.2^ which is also used by the 
Particle Data Group. 

1.5.3.1 Limits from short-scale gravity experiments 

Short-scale gravity experiments provide limits on the size (i.e. the compactification ra- 
dius) of the extra dimensions and relations like those in Table lL2l can be used to convert 
these to limits on the fundamental Planck scale (in a particular compactification scheme). 
Seeking to find deviations from Newtonian gravity at short distances is difficult because 
at these distances gravity is not the dominant force, even between non-magnetic, electri- 
cally neutral interacting bodies. This means it can be necessary to fully understand the 
Casimir force (taking into account finite conductivity, surface roughness and non-zero 
temperature corrections) and van der Waals forces which become important at these 
length scales. There are also experimental considerations in seeking to isolate the ex- 
periments from mechanical, acoustic and thermal vibrations as well as electromagnetic 
fields. 

The best constraints at present come from modern versions of the Cavendish torsion 
experiment. These experiments now show no deviation from Newtonian gravity down 
to about 150 /im [33] which corresponds to M^,(4+2) > 1.8 TeV. 

1.5.3.2 Limits from collider experiments and future sensitivity 

The signature of real graviton production in the process e+e^ G7 is a photon and 
missing energy. Searches in this channel at the LEP collider were able to set limits on the 

'^These different conventions have unfortunately led to much confusion in the literature, and explain 
some apparent inconsistencies between the limits mentioned here and those in some other recent reviews, 
for example [36]. 
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fundamental Planck scale, usually given in the M^, convention. The best pubhshed limits 
range from M^4+2) > 1-1 TeV through M,(4+4) > 0.68 TeV to M^(4+6) > 0.51 TeV [37], 
but some reviews quote values from the later higher energy running at LEP which gives 
M,(4+2) > 1.38 TeV and M,(4+6) > 0.58 TeV (DELPHI collaboration, quoted in [38]) or 
even M^(4_|_2) > 1.45 TeV and M^,(4_|_6) > 0.61 TeV (L3 collaboration, in [36]). 

The Tevatron Run I limits only become competitive with the LEP II limits for the 
higher values of n. The best limits from the Tevatron so far come from the CDF 
collaboration; they find, at 95% confidence level, that M*(4^2) > 1-00 TeV but that 
M*(4+6) > 0.71 TeV [39]. These limits are obtained from a study of the 'one or two 
high-energy jets + missing tranverse energy' signature (from processes like gg —>■ Gg, 
qg — > Gq and qq — > Gg). This is found to give a better sensitivity than the 'photon + 
missing i^x' signature (from qq G7)' since rates are much lower than in the jet case. 
However this signature could be used as confirmation in the event of a discovery in the 
jet channel. It is expected that the Tevatron Run II will be able to improve slightly 
upon these limits. 

It should be possible to reach much larger values of by studying the single- 
jet signature at the LHC (the dominant sub-process is qg Gq). With 100 fb~^ 
of integrated luminosity, the ATLAS collaboration expects to be able to probe M* in 
the following ranges: M*(4+2) =4.0-9.1 TeV, M*(4+3) =4.5-7.0 TeV and M*(4+4) =5.0- 
6.0 TeV [40]. The lower limits reflect the fact that the theory is only an effective one 
and breaks down at energies ~ M*. If the fundamental Planck scale is so low that a 
significant number of events would have V§ > then it is not possible to accurately 
extract information. This also means that nothing can be said for the n > 4 cases — in 
efi^ect the LHG would have too high a centre-of-mass energy (this isn't an issue in the 
pessimistic case in which we are just setting exclusion limits on M^). The upper limits 

'The 'Z-h missing E't' signature is less significant because the Z can only be observed via its leptonic 
decay. 
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are due to the signal no longer being observable over the SM background. In the ranges 
given, both M^, and n could be extracted at the LHC, but this would require running at 
two different centre-of-mass energies. 

Some limits can also be placed on ADD models by considering virtual graviton emis- 
sion. However the sum of KK states involved is found to be divergent; this means a 
cut-off has to be introduced which makes the process sensitive to ultra-violet effects. 
Hence although it is possible to place limits on this cut-off scale (in various different 
conventions) it is not possible to convert them to direct limits on M^,. 

1.5.3.3 Limits from cosmic rays 

Models with large extra dimensions would allow the production of miniature black holes 
at TeV-scale energies — a topic which will be properly introduced in Chapter 01 Ultra 
high-energy cosmic rays interact with centre-of-mass energies up to ~ 400 TeV and 
so could produce black holes; their decays could then be observed by experiments like 
AGASA.j Although there are some uncertainties involved, ref. [41] finds limits on the 
fundamental scale which, for the larger values of n, are more restrictive than present 
collider limits (e.g. M,(4+4) >1.3-1.5 TeV and M,(4+7) >1.6-1.8 TeV). 

1.5.3.4 Astrophysical and cosmological constraints 

Some of the strongest constraints on these extra-dimensional models are astrophysical 
and cosmological. These put direct lower bounds on Mp(4+„) which cannot absolutely 
rule out extra dimension models but make them unattractive for two reasons: firstly the 
whole motivation for the theory is lost if the fundamental Planck scale is constrained to 
be significantly larger than the electroweak scale; secondly there would be no obvious 
phenomenology at the next generation of colliders. 

■'This provides reassurance to those who are unconvinced by the theoretical arguments and worry 
that black holes produced at the LHC might destroy the world. 
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There are however various ways of managing values of Mp(4+„) up to of order 10 TeV. 
This is because ahhough A ~ Mp(4_|_„), the exact form of this relationship depends on 
the theory being proposed. For example if the extra dimension model is implemented 
within a type-I string theory, it can be shown that the scale of new physics (A) is given 
by the string scale (Mg) but that this can naturally be an order of magnitude smaller 
than Mp(4+„). Hence the hierarchy problem is naturally resolved, and new physics is 
expected at energy scales accessible at the LHC. 

The first hmits we consider here are those from the supernova remnant SN1987A. The 
rate at which the remnant can lose energy into the extra dimensions by graviton emission 
must not violate the observed cooling rate which is assumed to be mainly due to neutrinos 
(the flux of which has been measured by Kamiokande and 1MB). The predominant 
mechanism for graviton emission is known as 'gravi-strahlung' (N + N^N + N + G) 
and is enhanced over other processes by strong interaction effects. The graviton emission 
rate is found to be proportional to 1/ (Mp(4+„))"'+^ so this can be used to set a hmit on the 
fundamental Planck scale. Using the convention once again this gives M=^(4+2) >25- 
45 TeV and M*(4_|_3) >2-4 TeV [42,43] depending on the assumptions made and exactly 
what value is used for the supernova temperature. 

There are also limits on the fundamental Planck scale from the distortion of the 
cosmic diffuse background radiation spectrum. The gravitons produced both by gravi- 
strahlung and by neutrino annihilation can decay into SM particles and for the low-mass 
KK modes G ^ 77 is kinematically favoured. Observations of the cosmic diffuse gamma 
radiation by EGRET and COMPTEL show that gravi-strahlung must account for less 
than about 1% of the total energy emitted by all the supcrnovac that have exploded 
during the history of the universe. This is calculated to set limits of M*(4+2) > 70 TeV 
and M*(4+3) > 5 TeV [44] (more recent calculations of the energy emitted by KK states 
trapped in neutron star haloes suggest M^(4_)_2) > 380 TeV and M^(4_|_3) > 24 TeV [45]). 
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Other authors have considered the contributions to the gamma radiation background 
from the gravitons produced by neutrino annihilation and obtained the hmits M*(4_|_2) > 
90 TeV and M*(4+3) > 4 TeV [46]. Tighter hmits can be obtained by assuming that the 
universe doesn't enter the radiation epoch until after it has been reheated by the decay 
of a massive scalar field (or by some other means of entropy production) . 

Other limits are derived by ensuring that the production of gravitons will not cause 
the universe to become matter dominated too early — this would lead to a value for 
the age of the universe which was too low. The limits are M*(4+2) > 70 TeV and 
M*(4_i_3) > 6 TeV [47] with the assumption that the universe can be considered as 'normal' 
up to temperatures of 100 MeV. 

Finally there have been Hubble Space Telescope observations that the surface temper- 
atures of some older neutron stars are higher than conventionally expected. The expla- 
nation for this could be that they are heated by the decay products of gravitons trapped 
in the neutron star haloes during the supernova collapse. However the requirement that 
this heating is not excessive puts limits on the fundamental scale: M*(4+2) > 1400 TeV 
and M*(4_,_3) > 50 TeV [45]. The equivalent hmit for n = 4 is of order a few TeV. These 
are clearly the most constraining of the bounds on the fundamental scale although they 
are subject to various cosmological uncertainties. Even so, n — 2 almost certainly seems 
to be ruled out and n = 3 is disfavoured. 



Chapter 2 

Exotic Heavy Leptons at the LHC 



2.1 Introduction 

We have already seen in Chapter ^ that supersymmetry is a popular theory for new 
physics beyond the Standard Model. It is seen as theoretically attractive, solves the 
hierarchy problem and improves the unification of the three gauge couplings at high 
energy scales. 

'Intermediate scale' models which can be motivated on various astrophysical and 
cosmological grounds (see ^l.Hj) are constructed so that the gauge couplings unify at 
~ 10^^ GeV, an intermediate energy scale. This intermediate scale unification can be 
achieved by including new leptons as part of extra supermultiplets added to the MSSM. 
Although this choice is not unique, the extra leptons in [15] are two left-handed SU(2) 
doublets and three right-handed singlets (all with vector-like copies). 

By assumption there are no new Yukawa couplings for the model in [15]. This 
means that the lightest heavy lepton is stable and the others will decay into it. This is 
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possible because the tree-level mass degeneracy of the charged and neutral leptons will 
be destroyed by electroweak symmetry breaking. 

As in the MSSM, the renormalization group equations can be used to determine the 
spectrum of supersymmetric particles which will have phenomenological implications 
for the LHC A study has also been performed [48] investigating how to distinguish 
intermediate scale models from other supersymmetric string scenarios at a future linear 
collider. However a characteristic feature of this model is the existence of the new 
leptons and so this work concentrates on the phenomenology due to the extra leptons 
themselves. More specifically it is considered whether it will be possible to detect the 
charged heavy leptons at the LHC, and for what range of masses. 

Previous limits on the masses of new quasi-stable charged leptons have been limited 
by the energy available for particle production in lepton colliders. The gauge unification 
arguments in intermediate scale models mean that the extra leptons are expected to have 
masses in the TeV range. At the LHC enough energy should be available to produce and 
detect these exotic heavy charged leptons. Since the charged leptons are expected to be 
quasi-stable they will decay outside the detector; therefore the experimental procedure 
is different to that in [49,50] where limits are obtained by considering the decay products 
of the leptons. 

This chapter starts f ^2.2j) by reviewing the existing cosmological and experimen- 
tal limits on heavy particles to ensure that leptons such as those incorporated in this 
intermediate scale model are not ruled out. In section the theoretical input to the 
HERWIG Monte Carlo event generator to take account of new heavy leptons is described. 
Then, in §2.4[ the use of a 'time-of-flight' technique to detect the leptons at the LHC 
is considered: heavy leptons will arrive at the detector significantly later than relativis- 
tic particles, and the prospects for using this time delay as a method of detection are 
studied, for a range of masses. There is also a study of how to distinguish such leptons 
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from scalar leptons on the basis of their different angular distributions. The results and 
conclusions are presented in ^2.51 and ^2.61 respectively. 

2.2 Cosmological and experimental constraints 

The mass degeneracy of the charged and neutral leptons will be broken by radiative 
corrections. The electromagnetic self-energy correction for the charged lepton ensures 
that its mass will be greater than that of the neutral lepton [51]. This means that 
the charged lepton can decay to its neutral partner producing either a real or virtual 
W, depending on the mass difference. The neutral heavy lepton can be assumed to be 
totally stable (i.e. it has a lifetime orders of magnitude longer than that of the universe) 
because of the assumption that there are no new Yukawa couplings. 

There are a number of cosmological and experimental limits on leptons (particularly 
charged leptons) as a function of both mass and lifetime: 

• The relic abundance of the leptons (calculated from the self-annihilation cross 
section) must not 'over-close' the universe, i.e. provide more than the critical energy 
density (~ 10~^ GeV cm^'^) which is presumed to be accounted for by dark matter 
[52,53]. 

• A stable, charged lepton must have a low enough relic abundance for it not to 
have been detected in searches for exotic heavy isotopes in ordinary matter (see 
e.g. [54]). 

• The massive lepton and any decays it may have must not significantly affect nu- 
cleosynthesis or the synthesized elemental abundances [55,56]. 

• If the lepton decays before the recombination era (at ~ 10^^ s) it must not distort 
the cosmic microwave background radiation [55]. 
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• If the lepton has a longer lifetime it must not contradict the limits from gamma-ray 
and neutrino background observations [57,58]. 

• The mass and lifetime of the new leptons must not be such that they would have 
been detected in a previous collider experiment [59, 60] . 

Most of these limits are summarized in [61] and no constraints on charged leptons 
are found for lifetimes less than ~ 1 s, even for masses up to the TeV scale. The mass 
splitting between the neutral and charged leptons is expected to be at least ~ 10 MeV [51] 
which means that the lifetime of the charged lepton is expected to be shorter than 1 s. 

This means the only relevant constraints on the mass of the heavy leptons being 
studied arc the experimental limits. These arc to some extent mo del- dependent, but 
the LEP collaborations put a lower limit (95% confidence) on the mass of a long-lived 
charged lepton at around 102 GeV [59, 60] (the Tevatron experiment has not improved 
upon the LEP limits). 

2.3 Theoretical input to herwig 

The general-purpose Monte Carlo event generator HERWIG 6.3 includes subroutines for 
both neutral and charged Drell-Yan processes in a hadron collider [62-64] . As described 
in [63] the initial-state parton showers in Drell-Yan processes are matched to the exact 
0{as) matrix element result [65]. However the neutral Drell-Yan processes all use the 
approximation that the two fermions produced can be treated as massless. It is reason- 
able that this is a valid approximation for the SM quarks and leptons but this will not 
necessarily be the case for the proposed new heavy leptons discussed above which may 
have masses of order 1 TeV. 
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Figure 2.1: Feynman diagram for Drell-Yan production of heavy leptons. 



For this reason, full Born-level expressions for the Drell-Yan cross sections were 
derived, taking into account the masses of the produced leptons. The final results (for 
both neutral- and charged-current cases) were expressed in terms of vector and axial 
couplings to keep them as general as possible, although in all that follows Standard 
Model couplings will be assumed. The details of the derivation are in Appendix 1X1 

The derived expressions for both differential and total cross sections for the neutral- 
current Drell-Yan processes (qq Z°/7 L^L^, see Figure ITTjl are shown below. The 
notation is influenced by that already used in subroutines within HERWIG 6.3. 

The differential cross section is given by 



where 
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In these expressions and are the charges (in units of the electron charge) of the 
initial- and final-state particles respectively. The angle 9* is the angle between the 
out-going lepton (L~) direction and the incoming quark direction in the centre-of-mass 
frame. and refer to the centre-of-mass energy and momentum magnitude of the 
produced lepton. Similarly the subscript 4 refers to the produced anti-lepton. E is the 
energy of both the colliding quarks in the centre-of-mass frame (i.e. s = AE"^). 

The couplings dy and are related to the normal vector and axial coupling constants 
to the Z° (cv and ca) by relations like 



Equation ()2.1|) is a general expression which includes the Z°/7 interference terms; the 
massless case is retrieved by setting p^^ = E^ = E^^ = E (as well as = 7714 = 0). The 
normal charged-current Drell-Yan case (qq' L^L°) will not be relevant in this 

work but can be obtained by setting cy = ca = 1 and replacing gz, mz, Tz and q^^^ by 
(7w/v^5 mw, Tw and respectively. 

Integration of eq. ()2.ip gives the following expression for the parton-level cross section: 



In order to distinguish these heavy leptons from heavy supersymmetric partners of 
SM particles it would be necessary to consider the angular distribution of the produced 
particles. For the neutral-current production of left /right-handed scalars the differential 
cross section is also derived in Appendix ^ and is found to be 
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where 



D 



(s - mlf + m|r| 




(2.8) 



The couphngs and are defined as above and /il/r = 5'L/Rfi'z/2e where (yfL/R is the 
coupling to left/right-handed sleptons at the gauge boson-slepton-slepton vertex. 

The sin^ 6* angular distribution contrasts with the asymmetric distribution for the 
heavy leptons seen in eq. (j2IH); this will be studied in §2. 5. HI 

2.4 Implementation 

2.4.1 Modifications to herwig 6.3 

The above formulae were incorporated into a new subroutine added to HERWIG 6.3 
together with new particle entries (for the heavy leptons) and a new process code (IPROC) 
for the new process. 

Increasing the mass of the new leptons allowed the variation of cross section with 
mass to be studied — the mass was varied over a range from of order the top quark mass 
to a maximum of order 1 TeV. The lifetimes of the charged leptons produced were set 
to 1 s so that the number of decays occurring inside the detector will be negligible. 

The angular distribution of the produced leptons was studied over the same mass 
range to investigate the possibility of distinguishing heavy leptons from MSSM sleptons. 
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2.4.2 Time-of-flight technique 

The main aim was to consider the possible detection of charged heavy leptons at the 
LHC and for what mass ranges this might be practical. This work refers specifically to 
the technical specifications of the ATLAS detector [66-68]; however it is believed that 
similar results will be obtained for the CMS experiment [69]. 

The method used was a time-of-fiight technique as discussed in ref. [70] . This method 
utilizes the fact that, when compared to relativistic particles, there is a considerable time 
delay for heavy particles to reach the muon system. Heavy charged leptons like those 
being considered in this work will be detected in both the central tracker and the muon 
chambers, and from the measured momentum and time delay it is possible to reconstruct 
the mass. 

Imperfections in the time and momentum resolutions will broaden the mass peak. 
Uncertainty over which bunch crossing a particular detected particle comes from may also 
provide a background signal (from muons produced in Drell-Yan processes or in heavy 
quark decays). Plots from [68] show some differential cross sections for the processes 
most likely to produce muons which might be mis- identified as heavy leptons. These are 
reproduced in Figures EUl and EISl 

The studies presented here do not attempt to model the background, which it is not 
thought will be significant (there is no 'physics' background, only 'detector' background 
as discussed above), and use a simple cut [70] on the range of time delays allowed. 

The calculation of a reconstructed mass m from the time delay At and momentum 
is straightforward. For a lepton hitting the radial part of the muon spectrometer, the 
time delay with respect to a relativistic {P = v/c = 1) particle is given by 



At=—{E-p), 
Pt 



(2.9) 
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where r is the radius of the outer layer of the muon system, px the transverse momentum, 
and E and p the total energy and momentum respectively. Substituting for E in the 
energy-momentum invariant E"^ — = m? gives the result 



2 2 

m = p 



p^At f^^ ^ pjAt 



(2.10) 



It is also necessary to check for occasions when the lepton hits the end-cap of the detector 
by considering the magnitude of Pt/Pz- The muon system is modelled approximately 
as a cylinder of radius 10 m and with a half-length of 20 m [68]. A pseudo-rapidity cut 
requiring that \ri\ < 2.7 was applied to take account of the region close to the beam 
where particles cannot be detected. 

From equation ()2.10j) . the time delay and measured momentum for any particle 
detected in the muon system can be used to calculate its mass. As in [70] the time 
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delay was smeared with a Gaussian (although a width of 0.7 ns [68] rather than 1 ns was 
used) and a cut on the smeared time delay was applied such that 10 ns < At < 50 ns. 
Increasing the lower limit on At would reduce the efficiency but improve the mass 
resolution by removing many of the high-/5 leptons for which the time resolution is 
poor. The upper limit eliminates very slow particles which lose most of their energy 
in the calorimeter. An alternative upper limit of 25 ns, reflecting practical concerns, is 
discussed in §2.5.21 

The most important contribution to the momentum uncertainty Aj9 is due to the 
measurement error on the sagitta (the deviation from a straight line of a charged particle 
in a magnetic field). Apsag was taken into account by using 

^^ = l.lxlo-^ (2.11) 

where p is the total momentum in GeV. The constant in eq. ()2.11|) was obtained from [66] 
as an average over the different parts of the ATLAS detector. This expression should be 
valid near the discovery limit for heavy leptons when multiple scattering is insignificant — 
the maximum time delay cut is found to remove most of the low-/5 leptons for which 
multiple scattering would have more effect. However to allow investigations over a full 
mass range a multiple scattering term (Apms) was also incorporated with 

Apm, = 2 X IQ-^^p^ + m^. (2.12) 

The constant depends on the material distribution in the spectrometer — the above value 
is given in [71] for the ATLAS muon detector. 

The two errors in eqs. ()2.11|) and ()2.12|) were taken to be independent and hence 
combined in quadrature. A third possible contribution to momentum resolution — due 
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to fluctuation in the energy loss in the calorimeter — is neglected as it is always dominated 
by one of the other terms. 

Using the momentum and time resolution as described above it was possible to 
satisfactorily reproduce the mass resolution as a function of 13 for 101 GeV particles 
as given in [66]. 

The efficiency of the muon detection system is approximated to be 85%, independent 
of the particle momentum [66]. 

2.5 Results 
2.5.1 Backgrounds 

Any 'detector' background is greatly dependent on the details of the muon detection 
system. Such a background would come from muons produced either in the decays of 
heavy quarks or directly via the Drell-Yan process (both these processes have very high 
cross sections compared to that for heavy lepton production). If timing inadequacies 
mean that some of these muons are mistakenly thought to have come from an earlier 
interaction, i.e. they appear to have a large time delay, they could obscure the exotic 
lepton mass peak. For a background signal to look like a heavy lepton neutral-current 
signal, however, two opposite charge muons would have to be mis-identified at the same 
time which makes it very unlikely that background could be significant. 

To eliminate any possible background from heavy quark decays, a cut requiring the 
transverse momentum to be greater than 50 GeV was made to produce the plots in 
§2.5.21 Figures 12.41 and 12.51 show that a cut at 50 GeV should remove most of the 
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Figure 2.4: pT spectrum for muons produced in b decays. 




50 100 150 200 250 300 
Transverse momentum of heavy lepton (GeV) 




50 too 150 200 250 300 350 400 
Transverse momentum of heavy lepton (GeV) 



Figure 2.5: spectra for 175 and 250 GeV heavy leptons from Drell-Yan 
processes. 



muons from bottom quark decay without significantly reducing the exotic heavy lepton 
signal. It can be seen that as the mass of the exotic heavy lepton increases, the peak 
in the px spectrum becomes less well defined and shifts to higher values. Top quark 
decays can give rise to muons at high but the small top cross section, combined with 
the improbability of mis-identifying two muons, should suffice to eliminate this source 
of background. 
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2.5.2 Heavy lepton mass peaks 

The resuhs presented in this section show the reconstructed mass peaks for an integrated 
luminosity of 100 fb~^ after the cuts outhned above have been apphed. 

The possibihty of detection at the LHC is found to be entirely cross section limited — 
the decrease of cross section with increasing lepton mass is shown in Figure 12.61 The 
fraction of leptons passing the cuts actually slightly increases with mass due to the longer 
time delays. 

Figure 1213 for the neutral-current Drell-Yan process makes it clear that for masses up 
to about 1 TeV it should be possible to detect new charged heavy leptons at the LHC 
(particularly if a larger integrated luminosity can be obtained). A discovery criterion of 
10 leptons (from 5 events) in the mass peak is found to give a mass limit of 950 GeV for 
the detection of such leptons using an integrated luminosity of 100 fb^^. 

To reduce possible backgrounds the plots are based on events from which both leptons 
produced passed the cuts imposed. 
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Figure 2.7: Reconstructed mass peaks for exotic leptons with masses of 
175 GeV, 500 GeV, 750 GeV and 900 GeV. 



The long drift times in the muon chambers used by ATLAS allow the recording of 
events with time delays at arrival much greater than the 25 ns beam crossing interval of 
the LHC [68]; following [70] a maximum time delay of 50 ns is assumed here. However, 
to correctly associate the delayed track in the muon system with the event recorded in 
the other detector systems would require a specialized trigger based, for example, on 
the presence of high-px tracks in the inner detector (these would identify the true event 
time). Without such a trigger, the maximum allowed time delay would be 25 ns, and the 
discovery reach would be reduced to a lepton mass of 800 GeV. The statistical sample 
available for the angular distribution analysis in ^2 .5. 31 would also be reduced by a factor 
of about 2, depending on the lepton mass. 
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2.5.3 Distinguishing leptons from sleptons 

As mentioned in ^2. HI it should be possible to distinguish these new heavy leptons from 
heavy scalar leptons by studying the angular distribution with which they are produced. 
The two different angular distributions in the centre-of-mass frame are shown in eqs. ()2.H) 
and ()2.7j) . The application of the cuts described in ^2.4.2l as well as the effects of the mass 
and momentum resolution mean that the observed angular distributions are somewhat 
different to these. 

The forward-backward asymmetry for the exotic heavy leptons was not initially con- 
sidered because of the difficulty in a p-p collider of distinguishing the quark and anti- 
quark in the centre-of-mass frame. Hence only | cos 9* \ was used when the angular 
distributions were compared. Figure 12.81 shows the angular distribution for both exotic 
leptons and sleptons of mass 400 GeV. In an experimental situation, the exact quark 
direction is not known and so the angle used in this plot is not exactly 9* — instead 
it is the angle between the out-going lepton and the incoming proton. Examining the 
event records before and after the initial-state parton shower (in which the quark can 
acquire some transverse momentum) shows that the difference between the two angular 
distributions is small. 

A chi-squared test was applied to find the probability that a scalar model describes 
the heavy lepton data. The model was the angular distribution of a large Monte Carlo 
sample of scalars, and the total number of detected scalar pairs was rescaled to the 
expected number of detected pairs of heavy leptons. An 'average' data-set was obtained 
from a rescaled Monte Carlo sample of leptons, generated according to their angular 
distribution. Both these rescaled distributions had negligible theoretical errors compared 
to the (Poisson distributed) statistical errors on the model for 100 fb~^ of integrated 
luminosity. 
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Figure 2.8: Angular distributions of detected pairs of 400 GeV particles. 9* is 
the angle between the out-going particle and the incoming quark in the centre- 
of-mass frame. 



Figure 1231 shows that it will be possible to rule out scalar leptons at a confidence level 
of better than 90% up to a mass of 580 GeV. The angular distribution does not change 
significantly over the mass range of this plot — the decreasing probability of ruling out 
the scalar hypothesis as the mass increases is mainly because of the decreasing statistics. 

An alternative approach to ruling out the scalar hypothesis is to try to measure the 
forward-backward asymmetry of the lepton pair. To do this it is necessary to reliably dis- 
criminate between the quark and anti-quark directions; this is usually based (e.g. in [72]) 
on the assumption that the quark is a valence quark (with relatively high momentum) 
whereas the anti-quark quark (with lower momentum). A method based on 

the boost of the lepton pair was investigated for 500 GeV leptons and although the 
discriminator was found to be reliable for large boosts, the original asymmetry (~ 0.3) 
is not large enough to improve on the distinction discussed above. Cuts which remove 
events where the lepton pair has only a small boost (and hence the quark direction is 
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Figure 2.9: Probability at which it will be possible to rule out the scalar 
hypothesis by studying the angular distribution of the detected leptons. 



more uncertain) provided only a modest improvement because the increase in measured 
asymmetry is accompanied by a significant decrease in statistics. 



2.6 Conclusions 

This work was motivated by some of the new models of physics beyond the Standard 
Model in which the gauge couphngs unify at an intermediate energy scale (~ 10^^ GeV). 
These intermediate scale models can include additional leptons added to the Minimal 
Supersymmetric Standard Model and which are expected to have TeV scale masses. 

It is concluded that, assuming Standard Model couplings and a long enough lifetime, 
it should be possible to detect charged heavy leptons up to masses of 950 GeV with 
100 fb~^ of integrated luminosity. Background from mis-identified muons should be 
negligible and can be reduced by the application of a p™'" cut at about 50 GeV. It will 
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be difficult to use the angular distribution of the produced particles to distinguish them 
from scalar leptons for masses above 580 GeV because the falling cross section limits 
the available statistics. 

The absence of a specialized trigger in the inner detector could reduce the discovery 
hmit to 800 GeV and would reduce the statistics available for the angular distribution 
analysis. 

No present cosmological or experimental limits were found to rule out additional 
exotic leptons at the masses considered in this work. 



Chapter 3 

Introduction to Black Holes 



3.1 Black holes in four dimensions 

The work in this thesis relates to extra-dimensional black holes, but it is useful to first 
summarize some of the properties and features of black holes in four dimensions, as well 
as the history of the theory describing them. Black holes (or 'frozen stars' as they were 
originally called) can only be correctly described by Einstein's 1915 theory of general 
relativity (GR). Einstein himself was never completely convinced by the idea of black 
holes but, in 1916, Karl Schwarzschild apphed GR to a static non-spinning massive 
object and obtained his well-known metric showing that an object of mass M has a 
singularity at the Schwarzschild radius rs = 2GM/c^. Therefore an object with radius 
smaller than rs is a black hole with an event horizon at this radius. The Schwarzschild 
radius can be calculated for any object but, since most objects do not have the density 
required to be black holes, it is usually much smaller than the object's physical size — for 
example, the Schwarzschild radius of the earth is easily calculated to be approximately 
1 cm. 
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3.1 Black holes in four dimensions 



It is entirely coincidental that the exact horizon radius can be derived using an 
'escape velocity' argument which would be expected to be only qualitative since it is 
entirely Newtonian and fails to take into account the relativistic effects. However both 
Wheeler (1783) and Laplace (1796) independently used such an argument to postulate 
the existence of black holes as follows. An object of mass m propelled vertically upwards 
(with velocity v) from the surface of a mass M with radius r can only escape from the 
gravitational pull if 

1 9 GMm , , 
-mv^ > . 3.1 

2 r ^ ' 

If r decreases, v must be larger if the object is to escape. The radius at which escape is 
only just possible even if travelling at the speed of light is found by setting v — c. Hence 
we obtain 

2GM , , 

exactly as Schwarzschild derived using the later relativistic theory. 

The above argument — that not even light can escape from a black hole — would seem 
to suggest that black holes can only grow. However Hawking showed that the com- 
bination of general relativity and quantum mechanics predicts that black holes can 
evaporate by emitting Hawking radiation [73] . The mechanism for this is that a particle 
anti-particle pair is created at the event horizon; one has positive energy and escapes 
the black hole's gravitational attractions but the other negative energy one falls back 
into the black hole leading to a net decrease in its mass. 

Hawking found that the radiation spectrum is almost like that of a black body, and 
can be described by a characteristic Hawking temperature given by 

7h = , (3.3) 

where k is the Boltzmann constant. Throughout the rest of this work natural units are 
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used with the following definitions: 

fl = c = k = l. (3.4) 

This makes it possible to use energy units (usually GeV) to express distances (1 GeV~^ = 
1.97 X 10"^^ m), time intervals (1 GeV~^ = 6.58 x 10"^^ s) and temperatures (1 GeV = 
1.16 X 10^3 K). 

It is important to remember that, despite the wide acceptance of the theory, Hawking 
radiation has not been experimentally verified.^ In a four dimensional world we would 
only be able to observe Hawking radiation by studying astrophysical black holes. In fact 
although there are many astrophysical black hole candidates, it is difficult to uniquely 
identify them as black holes (for a review of some of the difficulties see [74] and the 
references therein) . The two main classes of black hole candidate are stellar black holes 
(with masses of order a few solar masses) in binary star systems, and super-massive 
black holes (~ 10^^ solar masses) at galactic centres which may power quasars and other 
active galactic nuclei. 

The masses of black hole candidates in binary systems can be estimated by studying 
the Doppler shift of emission lines from the secondary star and, since neutron stars have 
a maximum mass above which they will collapse to a black hole, this provides indirect 
evidence for the existence of astrophysical black holes. Additionally, there seems to be 
confirmation that the primary stars are black holes from the emission spectra of Doppler- 
shifted X-rays. These come from metals (e.g. iron) which are stripped off the outer layer 
of the secondary star and are then ionized due to the high temperatures of the accretion 
disc around the black hole candidate. 

'^Compare this with the comments about the Higgs boson in i il.ll 
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Observation of intense emission from super-massive black hole candidates is con- 
sistent with matter accretion by relatively small objects and hence also supports the 
black hole hypothesis. The most significant difference between black holes and other 
massive compact objects is the presence of an event horizon. Accretion models can 
provide evidence for both stellar and super-massive black holes by predicting how much 
energy falling on the central object from the accretion disc is 'lost' (as for a black hole 
event horizon) and how much is re-radiated and observed (as in the case of an ordinary 
compact object). 

Observing astrophysical Hawking radiation would absolutely verify the existence of 
black holes, but unfortunately the huge masses of these objects make this impossible. 
The large masses mean that even the highest black hole temperatures would be ~ nK — 
lower than the lowest temperatures achieved in the laboratory (~ fiK) and much lower 
than the temperature of the cosmic microwave background radiation (2.7 K). This means 
that the observation of astronomical Hawking radiation is impossible: it would take 
~ 10^^ years for a single photon from the closest black hole candidate to hit the earth 
and these photons would be swamped by the X-ray emission from matter being accreted 
into the disc surrounding the black hole. This is one of the reasons why the idea of 
observing Hawking radiation from miniature black holes produced in particle colliders 
is so exciting. 

In spite of the difficulties in studying black holes there is a huge body of work in 
the mathematical literature considering their properties including the emitted radiation, 
the singularities, and the effect of charge and angular momentum. A full review of this 
literature is not attempted here, but instead a few points particularly relevant for work 
on extra-dimensional black holes will be highlighted. 
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3.1.1 Grey-body factors 

Black holes are often spoken about as if the Hawking radiation is exactly like that 
of a black body. However this is not the case — Hawking's original work showed that 
the energy spectrum is in fact 'grey-body'. The spectrum is modified from that of a 
perfect black body by an energy-dependent grey-body factor; these factors are also spin- 
dependent and the main effect in four dimensions is to suppress the low-energy part of 
the spectrum (see Figure 1 of [75]). 

Perhaps counter-intuitively the grey-body factors for particle emission from black 
holes are equal to the absorption cross sections for the same particles incident on the 
black hole. However it is this property which ensures that the grey-body factors do not 
destroy the thermal nature of the black hole; equilibrium with a heat bath at the same 
temperature is still possible. Of course this equilibrium is unstable — black holes have 
a negative heat capacity so if the the temperature of the black hole is slightly below 
(above) that of the heat bath, there will be a net absorption (emission) of energy so the 
black hole mass will increase (decrease) with a corresponding further decrease (increase) 
in Hawking temperature. 

A substantial part of the work in this thesis (see Chapter ^ describes the numerical 
calculation of extra-dimensional grey-body factors. 

3.1.2 Super-radiance 

As mentioned above, the grey-body factor for a particular particle can be calculated 
from the absorption cross section of that particle incident on the black hole. However in 
the case of rotating black holes it was found that for incident bosons with low energies 
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it is possible for the absorption cross section to be negative [76-78]. This is equivalent 
to the reflection coeflicient being greater than unity which means that energy incident 
on the black hole can be amplified. This amplification was found to be greatest when 
the black hole was maximally rotating, and for a particle of spin s was found to occur in 
the mode with angular momentum quantum numbers i = m = max(s, 1). Perhaps most 
significant is the increase in magnitude of the super-radiant effect as s increases; peak 
amphfications quoted in [79] are 0.3% for scalars, 4.4% for gauge bosons but 138% for 
gravitons. In theory, the phenomenon of super-radiance makes it possible to construct 
a 'black hole bomb' [78] by surrounding a rotating black hole with a spherical refiector; 
the trapped electromagnetic energy would grow exponentially and eventually the mirror 
would explode. 

3.1.3 Back reaction 

The calculation of the Hawking radiation emitted by a black hole relies on the assumption 
that the back reaction of the black hole metric can be ignored. Essentially this is a 
requirement that the black hole mass Mbh (and equivalently the Hawking temperature 
Th) can be considered constant during the emission of each particle. For this to be a good 
approximation requires that the energy of the emitted particle satisfies u <S Mbh-^ This 
will be true for the majority of particle emissions provided that Th ^ Mbh (equivalent 
to the requirement that Mbh ^ Mp or that S'bh ^ 1 [80], where Mp and S'bh are the 
Planck mass and black hole entropy respectively). Regardless of the initial black hole 
mass, the radiation will no longer be described by Hawking's formula in the later stages 
of the decay. 

''The spectrum for the emission of higher energy particles (with u) ^ Mbh) must clearly be modified 
from Hawking's grey-body spectrum if for no other reason than that w > Mbh/2 is kinematically 
forbidden. 
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3.1.4 Black hole charges 

The 'no-hair' theorem states that the only meaningful characteristics of black holes are 
the charges related to continuous local symmetries. This means that we expect the 
electric charge of a black hole to be a meaningful quantity, but the same is not true 
for baryon and lepton number, for example. Black hole production and decay cannot 
violate Abelian gauge symmetries so electric charge must be conserved. Coloured black 
holes have to be treated slightly differently to electrically charged holes because of the 
non- Abelian nature of the gauge symmetry (of course colour charges are irrelevant for 
very massive 4D black holes, but they will be more significant for the extra-dimensional 
holes to be introduced in ^3.2^ . Although coloured black hole solutions can be found [81] 
the charge can never be observed by an observer at infinity, and so is not considered to 
characterize the black hole. 

The Hawking emission process occurs just outside the event horizon where the metric 
is described only by the charge, mass and angular momentum of the black hole. Hence 
the emission from, and decay of, the black hole cannot be dependent on any other 
properties of the black hole (or the matter from which it was originally made up). This 
is equivalent to saying that these properties (e.g. baryon number) are meaningless and 
will be violated in the decay. 

3.1.5 Information paradox 

There are many unanswered questions regarding the black hole information loss paradox, 
and these are intimately related to how the black hole decays, particularly in its final 
stages. The essence of the problem is that since a black hole evaporates to a mixture 
of thermal radiation (according to Hawking) any information which has fallen into the 
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black hole appears to be lost. Since mass, angular momentum and charge are the only 
observable characteristics of a black hole, any other information carried by matter which 
becomes trapped would appear to be lost. 

Many of the suggested solutions to this problem relate to what happens once the 
black hole mass becomes comparable with the Planck mass. Hawking's semi-classical 
calculations are no longer valid but without a quantum theory of gravity it is unclear 
what will happen to the remnant. Some argue that information is not lost but might be 
stored in a stable remnant of the black hole. A more popular idea is that information 
actually comes out of the black hole as non-thermal correlations within the Hawking 
radiation, whilst other theories accept and incorporate the loss of quantum information 
(for reviews of all these issues and possible resolutions, see [82-86]). 

3.2 Black holes in extra dimensions 

The extra- dimensional models of Arkani-Hamed, Dimopoulos and Dvali (ADD) [19- 
21] and Randall and Sundrum (RS) [28, 29] which were introduced in Chapter Q were 
motivated by the desire to explain the hierarchy problem — that is, the sixteen orders of 
magnitude difference between the electroweak energy scale and the Planck scale. 

In the standard version of those works the Standard Model fields are localized on a 
3-brane, which plays the role of our 4-dimensional world, while gravity can propagate 
both on the brane and in the bulk — the space-time transverse to the brane. In theories 
with large extra dimensions the traditional reduced Planck scale, M(4) ~ 10^^ GeV, 
is only an effective energy scale derived from the fundamental higher- dimensional one, 
-^p(4+n)5 through the relation 

M(4) - M^ll^^ R\ (3.5) 
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The above relation involves the volume of the extra dimensions, V ~ i?", under the 
assumption that R is the common size of all n extra compact dimensions. Therefore, 
if the volume of the internal space is large (i.e. if i? ^ £p, where ip = 10~^^ cm is the 
Planck length) then Mp(4_|_„) can be substantially lower than M(4). Only if Mp(4_(_„) can 
be as low as ~ TeV does such a model provide a solution to the hierarchy problem (and 
permit study at the next generation of colliders). 

In the regime r <^ R, the extra dimensions 'open up' and gravity becomes strong. 
Hence, Newton's law for the gravitational interactions in this regime is modified, with the 
gravitational potential assuming a l/r"+^ dependence on the radial separation between 
two massive particles. As a result there are limits on the fundamental scale from short- 
scale gravity experiments, but the most stringent constraints (at least for n = 2 and 3) 
were found in §1.5.31 to be cosmological, specifically from observations of neutron stars. 
Although there are many uncertainties these tend to exclude even the n = 3 case, while 
allowing models with Mp(4+„) ~ 1 TeV for n > 4. 

If extra dimensions with R ^ ip exist, then black holes with a horizon radius rh 
smaller than the size of the extra dimensions R are higher- dimensional objects centred 
on the brane and extending along the extra dimensions. It is these black holes which 
are the subject of the remaining chapters of this thesis. It has been shown that extra- 
dimensional black holes have modified properties, e.g. they are larger and colder than 
four-dimensional black holes of the same mass [87]. One striking consequence of the 
theories with large extra dimensions is that the lowering of the fundamental gravity 
scale allows for the production of such miniature black holes during scattering processes 
with centre-of-mass energy ^ Mp(4+„) [31,32,34,88-92]. 

There has already been a significant amount of work on production of black holes at 
future particle colliders [35, 93-123] and that is also the focus of the work in this thesis. 
However, as mentioned in ^1.5.3.31 a low fundamental energy scale would also make it 
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possible for miniature black holes to be created in the earth's atmosphere [41,124-138]. 
For the highest energy cosmic neutrinos (> 10'* TeV) the black hole production cross 
section from neutrino-nucleon interactions would be larger than the Standard Model 
interaction rate. There are several different techniques involved in detecting cosmic ray 
particle showers (for a review see [139]) but showers from black hole production and 
decay are expected to be characterized by particular particle content and multiplicity, 
and a favoured quasi-horizontal direction. Before the LHC becomes operational cosmic 
ray experiments should be able to improve their present limits [41] on the Planck scale or, 
more optimistically, observe black hole production. Recently it has been suggested [138] 
that black hole production could be one explanation for the hadron-rich 'Centauro' 
cosmic ray events. The theoretical results derived in Chapter HI should also be useful for 
the study of black hole decay in a cosmic ray context. 

3.2.1 Black hole production 

There has been much discussion in the literature (e.g. [31,89-91,140]) about what the 
cross section for black hole production is, but the consensus opinion that the geometrical 
cr ~ vrrg is valid seems to have been confirmed by the numerical work in [141]. More 
formally we can write 



In order to obtain the coefficient it is necessary to consider the black hole production 
process in detail. Various analytic techniques are available but for collisions with non- 
zero impact parameter in more than four dimensions, the collision must be modelled 
numerically [141]. This allows the determination of both F„ and femax; the maximum 
impact parameter for which black hole formation will occur. It is found that F„ is of 
order 1 (for example Fo = 0.647 and Fi = 1.883) and that femax/'^s ranges from 0.804 
(n = 0) to 1.37 (n = 7). 



a = 




(3.6) 
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Although a full numerical relativistic approach was necessary to obtain these values, a 
simple model of the collision can give good agreement. In four dimensions Thome's hoop 
conjecture [142] predicts the creation of a black hole in the case where any two partons 
from the colliding particles pass within the horizon radius corresponding to their centre- 
of-mass energy. In the extra-dimensional case there is probably an equivalent 'volume 
conjecture' [143,144], but a similar prediction is expected. 

Here we follow the approach in [145] and assume that a black hole will form if 
b < 2rh(MBH, J), where r-^ is the horizon radius for a black hole with mass Mbh and 
angular momentum J. This differs from that in [41] which assumes b < r^^lM^n, J), but 
provides less good agreement with the numerical results.'^ The angular momentum J 
can be usefully related to a dimensionless rotation parameter a* by 

{n + 2)J 
2rhMBH 

There will be more discussion of this in ^4.51 but rh is found to be related to rs by 

rh = (3.8) 

(l + a2)„+i 

It is therefore possible to determine 6max by setting b = 2rh and assuming that J = 
6Mbh/2. This gives 

1 

n\ 2 n+l 

"■"'^ r.. (3.9) 



2 

which is in excellent agreement with [141] for n > 1. It is then possible to estimate Fn 
(assuming cr = vrfo^g^^) and again the agreement is very good. 

Above we have assumed that the mass of the black hole will be equal to the centre-of- 
mass energy of the colliding partons. It is possible to put a lower bound (dependent on 



'^Reference [41] was published before the numerical results in [141], whereas [145] was published 
afterwards. 
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the ratio 6/&max) on the fraction of the energy which is trapped by the black hole. Lower 
bounds were derived for head-on {b = 0) collisions in [146-150] (for four dimensions) 
and in [140] (for the extra-dimensional case). However the value when the impact factor 
is close to its maximum is most important, since this is expected to give the largest 
contribution to the cross section. The 4D case was studied in [140] and then the extra- 
dimensional case in [141] showing that the bound on the fraction of trapped energy can 
be as low as ~ 0.1 for larger values of n. A significant amount of angular momentum 
could also be lost at this stage. More work is required in this area since if this bound is 
saturated the black hole decay discussed in this work would be significantly affected [151] 
since only if the black hole mass is larger than a few times the fundamental Planck mass 
can these objects be treated semi-classically. 

At this point it is worth stressing the difference between black holes as a signature 
for extra dimensions and the graviton emission signatures briefly mentioned in §1.4.1.11 
and §1.5.H.2I The graviton emission signatures are necessarily in the 'cis-Planckian' or 
'sub-Planckian' regime (V^ <^ Mp(4+„)) where the low-energy effective theory is valid 
and graviton emission can be reliably calculated with a perturbative expansion. This 
means that if Mp(4_|_„) is too low, the ability of the LHC to extract information from these 
signatures will be severely limited. In contrast, black holes would be produced in the 
'trans-Planckian' or 'super-Planckian' regime (v^ ^ Mp(4+„)) where a semi-classical 
approach is valid. If Mp(4_(_„) is too high it will not be possible to investigate this 
regime at the LHC, although it is hard to make a good estimate of what the minimum 
value of V^/Mp(4+„) should be if we are to be confident of examining trans-Planckian 
physics. Black holes are not the only feature of such a regime — for colliding partons 
with larger impact parameter (i.e. large compared to the horizon radius corresponding 
to their centre-of-mass energy) there will instead be gravitational elastic scattering with 
small momentum transfer. This can be calculated using the eikonal approximation 
and would result in jet-jet production close to the beam and with high centre-of-mass 
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energy [32]. There have been some preliminary experimental studies of this by the 
ATLAS collaboration which suggest that it might be possible to extract information on 
n and Mp(4+„) from these signatures [152]. 

Another possibility in the trans-Planckian regime is the production of 'black' p- 
branes (p-dimensional, spatially extended solutions of gravitational theories in extra 
dimensions). The production and decay of p-branes would exhibit some similar features 
to black holes (which can be considered as 0-branes) but the details of the decay are 
more model-dependent. Although the cross section for the production of p-branes can 
dominate that for black holes this is only found to be the case in compactifications 
with a mixture of small and large extra dimensions [100,104,153,154]. Qualitatively 
the explanation for this is that symmetric compactifications (with all extra dimensions 
of a similar size) favour the production of 0-branes (spherically-symmetric black holes) 
whereas non-spherically-symmetric p-branes are more likely to be produced in models 
with an asymmetric compactification, e.g. [155]. There will be no further discussion of 
p-brane production in this work, although it is an area in need of more study since string 
theory would seem to favour the existence of more than one compactification scale. 

The trans-Planckian and cis-Planckian energy regimes in which it might be possi- 
ble to extract parameters in a relatively model-independent way are separated by the 
Planckian regime in which quantum gravity effects become important. A theory of 
quantum gravity is required to predict cross sections and experimental signals in this 
regime. One possibility is string theory; this would introduce another energy scale, the 
string scale Mg, which could naturally be slightly smaller than Mp(4_|_„) (see ^1.5.3.4|) . 
Therefore, in the Planckian regime it might be possible for the colliding partons to be 
excited into string modes or even very excited 'string balls'. 
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3.2.2 Black hole decay 

Once produced, these miniature black holes arc expected to decay almost instantaneously 
(typical lifetimes are ~ 10~^^ s).'^ According to refs. [31,91], the produced black holes 
will go through a number of phases before completely evaporating. 

• Balding phase : The black hole emits mainly gravitational radiation and sheds the 
'hair' inherited from the original particles, and the asymmetry due to the violent 
production process.^ 

• Spin-down phase : The typically non-zero impact parameter of the colliding partons 
leads to black holes with some angular momentum about an axis perpendicular to 
the plane. During this phase, the black hole loses its angular momentum through 
the emission of Hawking radiation [73] and, possibly, through super-radiance. 

• Schwarzschild phase : A spherically-symmetric black hole loses energy due to the 
emission of Hawking radiation. This results in the gradual decrease of its mass 
and the increase of its temperature. 

• Planck phase: The mass and/or the Hawking temperature approach the Planck 
scale. A theory of quantum gravity is necessary to study this phase in detail but 
it is suggested that the black hole will decay to a few quanta with Planck-scale 
energies [31]. 

As in the 4-dimensional case [156] , it is reasonable to expect that the Schwarzschild 
phase in the life of a small higher- dimensional black hole will be the longest one, and will 
account for the greatest proportion of the mass loss through the emission of Hawking 
radiation. 

'^This is only true in the ADD model; in the RS model the black holes can be stable on collider time 
scales [97]. 

''This phase can also be considered as a production phase since it is predominantly the emission of 
the 'junk energy' which is not trapped by the event horizon as the black hole forms. 
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3.2.3 Experimental signatures 

The phases of black hole decay described in the previous section combine to produce 
distinctive experimental signatures which make it unlikely that black hole events would 
be mistaken for many other processes. Some of these features are outlined below: 

• There is a very large total cross section, particularly at high centre-of-mass energies 
(because the parton- level cross section grows with energy); 

• The ratio of hadronic to leptonic activity is roughly 5:1 and as a result the amount 
of energy visible in the detector is large (the relative emission into gravitons is also 
relatively small); 

• Most events have a relatively high multiplicity with many hard jets and hard 
prompt leptons; 

• As the total event transverse energy increases the average multiplicity increases 
and the average energy of each primary emitted parton decreases — a manifestation 
of the infra-red ultra-violet connection of gravity; 

• Events have a high sphericity since most black holes are produced almost at rest; 

• Hard perturbative scattering processes are suppressed by the non-perturbative 
black hole production at high energies. 

It has been suggested [112] that the high temperature of the black hole (» Aqcd, 
the scale above which physics becomes perturbative) might cause a quark-gluon 'chro- 
mosphere' to form around the decaying black hole; the Hawking radiation would then 
thermalize leading to a suppression of the number of hard hadronic jets. However this 
would require ~ 10 quarks to be emitted as primary partons [112, 151] which is unlikely 
to be a problem for many black hole events at the LHC because of the lower entropy 
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of higher-dimensional black holes (see the multiplicity discussion in ^4.3.5p . For future 
higher energy particle colliders, the black hole signatures are more likely to be modified 
although there would still be hard hadronic jets from W, Z and Higgs bosons (as well 
as r leptons) which would penetrate the chromosphere before decaying [157]. 

Experimental cuts proposed to isolate black hole events with negligible SM back- 
ground include requiring the total energy deposited in the calorimeter to be > 1 TeV 
and there to be > 4 jets with energies above 100 GeV (including a high-energy lepton 
or photon which is useful for triggering) [34]. The theoretical work in Chapter 0] and the 
event generator described in Chapter El will eventually make it possible for more quanti- 
tative predictions to be made about typical black hole events, and for experimental cuts 
to be tailored accordingly. 

3.2.4 Proton decay 

Proton decay is always a potential problem for models which describe physics beyond the 
Standard Model. The present limits on proton lifetime require that, at a 90% confidence 
level, the partial lifetime for the decay p e^vr^ is greater than 1.6 x 10^^ years [158]. 
In general introducing new physics at a scale of A can allow the decay of the proton 
(mass nip) with a width F ~ {m^/A'^) unless the operators in the Lagrangian which 
would allow this are explicitly forbidden or at least heavily suppressed. 

To stop present experimental limits on proton decay being violated in models of new 
physics it is usually necessary to artificially impose symmetries which forbid the prob- 
lematic higher dimensions operators (as in the Standard Model this does not completely 
forbid proton decay since B and L symmetries can be broken by radiative effects). 

In extra dimension models there is new physics at the fundamental Planck scale (so 
A ~ Mp(4+„) ~ 1 TeV) and the possible production of black holes in extra dimension 
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models means that gravity itself could be responsible for proton decay. This could occur 
through the production of virtual black holes in processes like q + q — > Black Hole —>■ 
q + i + . . ., where the ellipsis represents possible additional particles (e.g. gravitons, 
gluons, photons or neutrinos). The consequence of the 'no hair' theorem is that the 
only requirement for such processes is the conservation of charge, energy and angular 
momentum. 

The usual method of suppressing proton decay in models of new physics is by im- 
posing symmetries on the high-energy theory — for example, R-parity in supersymmetry. 
However in the extra-dimensional case the high-energy theory is an unknown theory of 
quantum gravity so it is not clear whether it is possible to impose such a symmetry to 
suppress but not totally forbid proton decay. The safest approach is to explicitly forbid 
proton decay by imposing a symmetry on the low-energy effective theory but again this 
is not trivial. A global symmetry for baryon number does not seem to be viable since 
the black hole decay is not guaranteed to conserve quantum numbers associated with 
global symmetries. This leaves discrete gauge symmetries involving combinations of B 
and L as apparently the only possibilities, since it is argued in [159] that quantum num- 
bers associated with discrete gauge symmetries as well as continuous gauge symmetries 
will be conserved in the decay of black holes. Care must be taken to ensure that any 
such symmetries introduced are anomaly free, and that they still allow the generation 
of baryon asymmetry and neutrino masses. 

Of course the above argument has made the usual assumption that quantum numbers 
associated with global symmetries can be violated, i.e. information is lost in the black 
hole decay process. Some would argue [160] that the solution of the Hawking information 
loss paradox discussed in ^3.1.51 might involve the conservation of global charges like 
baryon number in which case many of the above problems would be avoided. 
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Some of the TeV~^ extra dimension models introduced in ^1.4.1.21 provide another 
way of avoiding proton decay issues. The models involve physically separating the quarks 
and leptons [24,25] which exponentially suppresses their wave-function overlap and hence 
also the possibility of proton decay. To obtain normal SM physics the Higgs and gauge 
bosons must propagate in the bulk (or the 'thick' brane) between the lepton and quark 
fields. This requires the mass of the lowest KK modes for the bosons to be high enough 
to avoid the experimental constraints from the precision electroweak measurements, and 
translates to a lower limit on Mp(4+„) (at least of order 20 TeV [35]). As a result such 
models are disfavoured phenomenologically because the high value of Mp(4+„) rules out 
black hole production at the next generation of colliders. 

3.3 Basic formulae and assumptions 

Although extra dimension models have only recently become popular among particle 
physicists, Myers and Perry [161] had previously worked on black holes in higher- 
dimensional space-times. They considered the form of the gravitational background 
around an uncharged (4+r;,)-dimensional black hole. In the non-rotating (Schwarzschild) 
case the line-element is found to be given by 



ds'^ = -h{r) df + h{r)-^ dr^ + d^l 



'2+n 5 



(3.10) 



where 




(3.11) 



and 



dnl 



'2+n ~ 



dol^^ + sin^ e, 



n+1 



del + sin^6'„ ( ... + sin' 
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In the above, < if < 27i and < < vr, for i = l,...,n + 1. As shown in [161], 
the extension of the usual 4-dimensional Schwarzschild calculation gives the following 
horizon radius: 



where Mp(4+„) is the fundamental (4 + n)-dimensional Planck scale in convention 'd' of 
Table Ol 



The black holes being considered in this work are assumed to have horizon radii 
satisfying the relation £p <C rs <^ R. The former inequality guarantees that quantum 
corrections are not important in the calculations, while the latter is necessary for the 
black holes to be considered as higher-dimensional objects (i.e. the curvature of the 
extra dimensions can be ignored on the scale of the black hole). The tension of the 
brane on which the black hole is centred is assumed to be much smaller than the black 
hole mass which means it can be neglected in this analysis, and a zero bulk cosmological 
constant is also assumed. This is essentially equivalent to the assumption of the ADD 
scenario although much of what follows can also be applied to the RS models provided 
the bulk cosmological constant is small; this corresponds to a small warp factor and so 
the space-time can be considered as almost spherically symmetric. If this is not the case, 
singularity problems emerge and the black hole solutions lose their spherical symmetry 
(the horizon becomes flattened on the brane in a pancake shape [162]). 

A black hole of a particular horizon radius is characterized by a Hawking temperature 
related by 




(3.13) 



Th 



(n + 1) 
47rrs 



(3.14) 
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and the black hole entropy is given by 



4.rsMBH ^/n + l\i^_ (3.15) 
n + 2 \n + 2 Tn ^ ^ 



The emitted Hawking radiation is almostlike that of a black body at this temperature. 
In fact the flux spectrum, i.e. the number of particles emitted per unit time, is given 
by [73] 

-^r- = exp (./Th) T 1 ' ^'-''^ 

while the power spectrum, i.e. the energy emitted per unit time, is 



^^-Z.^^ exp {u/T^) T 1 • ^^-^^^ 



In the above, s is the spin of the degree of freedom being considered and I is the 
angular momentum quantum number. The spin statistics factor in the denominator is 
— 1 for bosons and +1 for fermions. For massless particles |p| = and the phase-space 
integral reduces to an integral over uj. The term in front, a][\uj)^ is the grey-body factor 
which encodes valuable information about the structure of the surrounding space-time 
including its dimensionality. Chapter 0] will give more detail on how the grey-body 
factors are calculated. 

As the decay progresses, the black hole mass decreases and the Hawking temperature 
rises. The work performed here will allow some comments to be made on whether the 
usual 'quasi-stationary' approach to the decay is valid; if such an approach is correct 
then the black hole has time to come into equilibrium at each new temperature before 
the next particle is emitted. 

An important point to stress is that eqs. ()3.16|) and ()3.17|) refer to individual degrees 
of freedom and not to elementary particles, like electrons or quarks, which contain more 
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Particle type 


do 


di/2 


di 


Quarks 





72 





Gluons 








16 


Charged leptons 





12 





Neutrinos^ 





6 





rhoton 








2 


ZO 


1 





2 


W+ and W- 


2 





4 


Higgs boson 


1 








Total 


4 


90 


24 



Table 3.1: Degrees of freedom for different particle types emitted from a black 
hole. 



than one polarization. Combining the necessary degrees of freedom and their correspond- 
ing flux or power spectra, the relative numbers of different elementary particles produced, 
and the energy they carry, can be easily computed. Since the higher-dimensional black 
holes which might be produced at the LHC have relatively high Hawking temperatures 
(~100 GeV) in most cases it is possible for all elementary particles to be produced. The 
numbers of degrees of freedom dg are summarized in Table 13.11 These take account of 
the fact that for each massive gauge bosons one of the degrees of freedom comes from the 
Higgs mechanism. The numbers of degrees of freedom will be important in calculating 
total flux and power emission in §4.3.5[ and for the construction of the black hole event 
generator (see Chapter EJ. 

In reference [163] it was argued that the majority of energy during the emission of 
Hawking radiation from a higher-dimensional black hole is emitted into modes on the 
brane (i.e. Standard Model fermions and gauge bosons, zero-mode gravitons and scalar 
fields). This argument was based on their result that a single brane particle carries as 
much energy as the whole Kaluza-Klein tower of massive excitations propagating in the 



It is possible that di/2 = 12 should be used for the neutrinos since right-handed neutrinos, if they 
exist, will also be emitted by the black hole. 
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bulk. Essentially the reasoning is that the KK tower is simply the result of the projection 
onto four dimensions; from the extra-dimensional view-point there is only one graviton 
mode. Since there are many brane modes for the Standard Model particles the majority 
of energy is expected to be emitted on the brane. 

Chapter |3] of this thesis is concerned with the numerical calculation of the extra- 
dimensional grey-body factors. These will be a key part of the state-of-the-art black 
hole event generator described in Chapter |S| The event generator simulates both the 
production and decay of small black holes at hadronic colliders and, by using the new 
results for the grey-body factors provides estimates for the spectra and relative numbers 
of the different types of elementary particles emitted. 



Chapter 4 



Grey-body Factors in (4 + n) 



Dimensions 



4.1 Introduction 

The idea of grey-body factors, both for four- dimensional and extra-dimensional black 
holes, was introduced in Chapter IHl Extra-dimensional black holes are of interest be- 
cause, if scenarios like that of ADD [19-21] are realized in nature, they could possibly be 
produced at high-energy colliders. To understand the particle emission probabilities and 
spectra for these black holes it is necessary to calculate the extra-dimensional grey-body 
factors. 

For ease of reference, equations ()3.16|1 and ()3.17|) describing the flux and power 
spectra are reproduced below: 



1 c/"+3 



p 



(4.1) 



dt 



exp {u/Tn) T 1 (2vr) 
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and 



(4.2) 



dt 



exp (cu/Th) T 1 (27r)"+3 ' 



{(jj) is the grey-body factor'^ which can be determined by solving the equation of 
motion of a particular degree of freedom in the gravitational background and computing 
the corresponding absorption coefficient A^l\ The absorption coefficient is related to 
the grey-body factor by [164] 



As already mentioned, the grey-body factors modify the spectra of emitted parti- 
cles from that of a perfect thermal black body even in four dimensions [73]. For a 
4-dimensional Schwarzschild black hole, geometric arguments show that, in the high- 
energy {uj ^ Th) regime, E^l^j,''-'!^ oc (curh)^. Therefore the grey-body factor at high 
energies is independent of uo and the spectrum is exactly like that of a black body for 
every particle species [156,165-167]. The low-energy behaviour, on the other hand, is 
strongly spin-dependent. A common feature for fields with spin s = 0, 1/2 and 1 is 
that the grey-body factors reduce the low-energy emission rate significantly below the 
geometrical optics value [75,156]. The result is that both the power and fiux spectra 
peak at higher energies than those for a black body at the same temperature. As s in- 
creases, the low-energy suppression increases and so the emission probability decreases. 
These particle emission probabilities have been available in the literature for over twenty 
years [156,166,167]. 

These characteristics of the four- dimensional grey-body factors demonstrate why it is 
important to calculate the equivalent extra-dimensional factors. Studying their depen- 

'^The quantity (t^*' {ui) is alternatively called the absorption cross section. It is also common in the 
literature to refer to the absorption probability 1^^"^^, related to (T'f\ijj) through eq. (|4.3|) . as the 
grey-body factor. The notation F^^-* is sometimes used for jyl^^-'p. 




2n^(n+i)/2 Y\{n + l)/2] (2£ + n + 1) {I + n)\ 



(s)\2 
'I I 



(4.3) 
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dence on the dimensionahty of space-time is important if detected Hawking radiation 
emitted by small black holes is to be used to try to determine the value of n. Both 
the energy spectra and the relative emissivities of scalars, fermions and gauge bosons 
will be potentially useful in this task. The main motivation for this work was to obtain 
grey-body factors which could be used in the Monte Carlo event generator described in 
Chapter thereby allowing the possible determination of n to be studied in detail. 

The procedure for calculating grey-body factors needs to be generalized to include 
emission from small higher- dimensional black holes with rh < -R; there is an added com- 
plication in this case because black holes can emit radiation either on the brane or in the 
bulk. The emission of particle modes on the brane is the most phenomenologically inter- 
esting effect since it involves Standard Model particles. However gravitons will certainly 
be emitted into the bulk (since it is the propagation of gravity in the bulk which accounts 
for its apparent weakness and hence solves the hierarchy problem) and it is possible that 
a (4 + n)-dimensional black hole will also emit bulk scalar modes. These are potentially 
important because in some extra dimension models there are additional scalars propa- 
gating in the bulk; for example, supersymmetric models inspired by string theory involve 
graviton supermultiplets which can include bulk scalars (see e.g. [168, 169]). Numerical 
results for bulk and brane scalar modes are necessary to allow comparison of the total 
bulk and brane emissivities. This would allow the usual heuristic argument for brane 
dominance (because there are many more SM brane degrees of freedom than there are 
degrees of freedom in the bulk — see ^3.3|) to be confirmed or disproved. 

There have already been some analytic studies for extra-dimensional black holes [145, 
170,171] which have suggested that the grey-body factors can have a strong dependence 
on the number of extra dimensions. However it is known that results from the power 

'^Scalars propagating in the bulk would affect many of the constraints described in m.5.31 although 
in model-dependent ways since the couplings to SM particles would not be universal (unlike those of 
the graviton). 
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series expansions in these papers are only accurate in a very limited low-energy regime 
(especially in the case of scalar fields) and even the full analytic results are not necessarily 
reliable above the intermediate-energy regime. 

The work described in this chapter focuses on the emission of energy from uncharged 
(4 + ?2)-dimensional black holes; initially only Schwarzschild-like black holes are con- 
sidered although there are also some results for rotating holes. Section 14.21 gives some 
details of a master equation describing the motion of scalars, fermions and gauge bosons 
in the 4-dimensional induced background — the basic steps of Kanti's calculation [2] are 
given in Appendix El This is particularly useful in that it helps to resolve ambigui- 
ties present in similar equations which have previously appeared in the literature. The 
brane grey-body factors and emission rates are also defined in ^4.2| and analytic and 
numerical methods for computing these quantities are discussed. Exact numerical re- 
sults for emission on the brane (for scalars, fermions and gauge bosons) are presented 
in §4.31 This allows a comparison with the earlier analytic studies of the Schwarzschild 
phase as well as a discussion of the relative emissivities and their dependence on n. 
A calculation of the total power and fiux emitted (and hence the lifetime and number 
of emitted particles expected for extra-dimensional black holes) is also possible in this 
section. The emission of bulk scalar modes is studied in ^4.41 and the results for the 
relative bulk-to-brane emissivities are also presented in this section. The analysis pro- 
vides, for the first time, exact results for the relative bulk and brane emissivities in all 
energy regimes and for various values of n. Finally section 14.51 includes a discussion of 
the emission of scalars from rotating black holes; this allows some comments to be made 
on super-radiance for extra-dimensional black holes. The conclusions of the work in this 
chapter are summarized in §4.61 
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4.2 Grey-body factors for emission on the brane 

This section and ^4.HI focus on the emission of brane-localized modes, leaving the study 
of bulk emission and of relative bulk-to-brane emissivity until §4.41 

The brane-localized modes propagate in a 4-dimensional black hole background which 
is the projection of the higher-dimensional one, given in eq. ()3.10|) . onto the brane. The 
induced metric tensor follows by fixing the values of the extra angular co-ordinates 
{6i = 7i/2 ioT i > 2) and is found to have the form 



The grey-body factors are determined from the amplitudes of in-going and out-going 
waves at infinity so the essential requirement is to solve the equation of motion for a 
particle propagating in the above background. For this purpose, a generalized master 
equation was derived by Kanti (see Appendix EI) for a particle with arbitrary spin s; this 
equation is similar to the 4-dimensional one derived by Teukolsky [172]. For s =1/2 and 
1, the equation of motion was derived using the Newman-Penrose method [173,174], 
while for s = the corresponding equation follows by the evaluation of the double 
covariant derivative g^'^D^D^ acting on the scalar field. The derived master equation is 
separable in each case and, by using the factorization 



ds^ = -h{r) dt^ + h{r)-^dr^ + {dO^ + sin^ 9 dip^) . 



(4.4) 



^-iut ^im^ R^{r) sS'Pie) , 



(4.5) 



we obtain the radial equation 




dRg 



dr 



) 



+ 




+ 2iujsr — 



isujr'^h' 



h 



+ s(A"-2)-,A, i?,(r) = 0, (4.6) 
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where A = hr"^ . The corresponding angular equation has the form 

1 d ( ^d,Sf\ ( 2ms cote o 2. ^ 

smfc* dQ \ dU J \ smO sm 6 

(4.7) 

where e^™''^ sS^{9) = sY^™'(i72) are known as the spin-weighted spherical harmonics and 
sXe is a separation constant which is found to have the value ^A^ = £{£+!) — s{s + l) [175]. 

For s = 0, equation ()4.(i|l reduces as expected to eq. (41) of ref. [170] which was used 
for the analytical study of the emission of brane- localized scalar modes from a spherically- 
symmetric higher-dimensional black hole. Under the redefinition Rs = A^^Pg, eq. ()4.()j) 
assumes a form similar to eq. (11) of ref. [171] which was used for the study of brane- 
localized fermion and gauge boson emission. The two equations differ due to an extra 
term in the expression of the latter one, which although vanishing for s =1/2 and 1 
(leading to the correct results for fermion and gauge boson fields) gives a non-vanishing 
contribution for all other values of s. Therefore, the generalized equation derived by 
Cvetic and Larsen [176] cannot be considered as a master equation valid for all types 
of fields.^ Hence the derivation of a consistent master equation was imperative, before 
addressing the question of the grey-body factors in the brane background. This task 
was performed by Kanti and led to eqs. ()4.6|) and ()4.7|) above. 

For the derivation of grey-body factors associated with the emission of fields from 
the projected black hole, we need to know the asymptotic solutions of eq. ()4.6|) both as 
r — Th and as r cxd. In the former case, the solution is of the form 

i^i*^) = a;;) a- e-^-'-* + 42 e^^^* , (4.8) 

'^A similar equation was derived in ref. [145] but due to a typographical error the multiplicative 
factor s in front of the A"-term was missing, leading to an apparently different equation for general s. 



\sr{e) = o, 
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where 

dr* _ 1 

We impose the boundary condition that there is no out-going solution near the horizon 
of the black hole, and therefore set A^^^^ = 0. The solution at infinity is of the form 



and comprises both in-going and out-going modes. 

The grey-body factor af''^ (uj) for the emission of brane- localized modes is related to 
the energy absorption coefficient Ae through the simplified relation 



ai'\uj) = ^^i2i+l)\Af^\^ (4.11) 



where from now on a 'hat' will denote quantities associated with the emission of brane- 
localized modes. The above relation follows from eq. ()4.3|) by setting n = since the 
emission of brane-localized modes is a 4-dimensional process. The absorption coefficient 
itself is defined as 

^(oo) ^(h) 

|^(.)|2^^_ Apt (4.12) 

''in -'in 

in terms of the out-going and in-going energy fiuxes evaluated either at infinity or at 
the horizon. The two definitions are related by simple energy conservation and lead to 
the same results. The choice of which is used for the determination of the absorption 
coefficient may depend on numerical issues. 

The fiux and power spectra (equations 1)4.11) and ()4.2j] respectively) of the Hawking 
radiation emitted on the brane can be computed, for massless particles, from the 4- 
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dimensional expressions: 

= y (2£ + , I, - ; (4.13) 

dt ^' ^ ' exp (cu/Th) T 1 277 ' ^ ^ 

Note, however, that the absorption coefficient A^^"* stiU depends on the number of ex- 
tra dimensions since the projected metric tensor of eq. ()4.4|1 carries a signature of the 
dimensionahty of the bulk space-time through the expression of the metric function h{r). 



4.2.1 Analytic calculation 

The grey-body factors have been determined analytically for the 4-dimensional case 
in refs. [156, 166, 167] both for a rotating and non-rotating black hole. In the (4 -|- 
rz) -dimensional case, refs. [170, 171] have provided analytic expressions for grey-body 
factors for the emission of scalars, fermions and gauge bosons from a higher-dimensional 
Schwarzschild-like black hole. Both bulk and brane emission were considered for an 
arbitrary number of extra dimensions n. Reference [145] presented analytic results for 
brane-localized emission from a Kerr-like black hole, but only in the particular case of 
n = 1. All the above results were derived in the low-energy approximation where the 
procedure used is as follows: 

• Find an analytic solution in the near-horizon regime and expand it as in-going and 
out-going waves so that the A^^"* coefficients can be extracted; 

• Apply the boundary condition on the horizon so that the wave is purely out-going; 

• Find an analytic solution in the far-field regime and again expand it as in-going 
and out-going waves; 
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• Match the two solutions in the intermediate regime; 

• Extract {Ael"^ and expand it in powers of ur^. 

The solution obtained by following the above approximate method is a power series 
in ccTh, which is only valid for low energies and expected to significantly deviate from 
the exact result as the energy of the emitted particle increases (this was pointed out 
in [75,156] for the 4-dimensional case). In [171], the full analytic result for the absorption 
coefficient (before the final expansion was made) was used for the evaluation of the 
emission rates for all particle species."^ The range of validity of these results, although 
improved compared to the power series, was still limited since the assumption of small 
coTh was still made during the matching of the two solutions in the intermediate regime. 
As the equations of motion are too complicated to solve analytically for all values of oor-^, 
it becomes clear that only an exact numerical analysis can give the full exact results for 
grey-body factors and emission rates. 

4.2.2 Numerical calculation 

There are various numerical issues which arise while trying to do the full calculation of 
the absorption coefficient and the complexity of these problems depends strongly on the 
spin s of the emitted particle. The usual numerical procedure starts by applying the 
relevant horizon boundary condition (a vanishing out-going wave) to Rs{r). Then, the 
solution is integrated out to 'infinity' and the asymptotic coefficient A^^^ is extracted 
in terms of which Af^ can be calculated. 

By looking at the asymptotic solution at infinity in eq. ()4.1()|1 we can see that, for the 
scalar case, the in-going and out-going waves are of comparable magnitude. Therefore 

'^This full analytic approach was pursued in [171] as a result of some of the numerical results presented 
in this thesis. 
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it is relatively easy to extract the coefficients at infinity and thus determine the grey- 
body factor. However for fields with non-vanishing spin, i.e. fermions, vector bosons and 
gravitons, this is not an easy task. First of all, different components carry a different 
part of the emitted field: the upper component \E'+s consists mainly of the in-going wave 
with the out-going one being greatly suppressed, while for the lower component \&_<( the 
situation is reversed (for a field with spin s 7^ 0, only the upper and lower components 
are radiative). Distinguishing between the two parts of the solution (in-going and out- 
going) is not an easy task no matter which component is used, and it becomes more 
difficult as the magnitude of the spin increases. In addition, the choice of either positive 
or negative s to extract the grey-body factor (i.e. using either the upper or the lower 
component) affects the numerical issues. 

If s is negative, then the horizon boundary condition is easy to apply because compo- 
nents of the in-going solution on the horizon will be exponentially suppressed, as can be 
seen from eq. ()4.8j) . However eq. ()4.1()|1 shows that negative s also means the out-going 
solution at infinity is enhanced by r~^^ with respect to the in-going one; this makes 
accurate determination of ^4^^^ very difficult. Hence the negative s approach is not used 
in this work. 

On the other hand if s is positive, it is easy to accurately extract A^^^ because the 
out-going solution at infinity is suppressed by r~^*. However, close to the horizon the 
out-going solution is exponentially smaller than the in-going one. Therefore the solution 
for Rs{r) can be easily contaminated by components of the out-going solution. This 
problem becomes worse for larger s and first becomes significant for s = 1. 

For n = 0, various methods (see e.g. [79]) have been used to solve the numerical 
problems which arise in the gauge boson case. The approach of Bardeen in [79] is also 
applicable for n = 1; some details of this method are briefiy outlined below. 
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The method starts from a mathematicaUy equivalent form of the radial equation in 
Kerr in-going co-ordinates:*^ 

A-'4- ( A^+i^^ - 2ir''uj^ - \2{2s + l)iujr - sA" + ,A^] i?, = , (4.15) 
dr \ dr J dr 

where ^A^ = ^{^ + 1) — s(s — 1). The asymptotic solutions are now different from those 
of eqs. ()4.8|) and ()4.10p : they are given by 

i?(^) = A^) + i?i;;A-e2'-^*, (4.16) 



and 

R'r' = B^'^,^Bt:^'--. (4.17) 



We can follow the method in Appendix A of [79] and apply the co-ordinate transfor- 
mation 

x=—-l. (4.18) 
rh 

This transformation is easily applicable for either n = or n = 1 since in both cases 
A" = 2. In the case s = — 1 (which is mathematically equivalent to the s = 1 discussed 
earlier) applying this transformation for n = gives 

xix + 1) , - 2iturh [x^ + 2a; + l) — ^ - \-2iur^{x + 1) + A] = , (4.19) 
dx^ dx 

where A = iiJL +1). The asymptotic solutions can now be written in terms of a;, with 
more care being taken over sub-leading terms in the in-going solution at the horizon: 

= bI^\i + aix + a2X^) + B^^.r^^^' x{l + hx + . . .) . (4.20) 



°This was derived, for general n, in [177] and for n = reduces to the result found in [79]. 
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It is easy to obtain tlie coefficients ai and 02 by substituting the in-going solution 
into the radial equation and comparing terms in different powers of x. This leads to 

ai = l + 2-^, (4.21) 

and 

a2 = 2(i_2ia;rh) ^ 2icjrh)ai - 2icurh] . (4.22) 

Instead of solving CR = 0, the equation CY = / is solved, where Y = R — {1 + aix). 
This means that 

f = 2{1 - 2iurh)a2X . (4.23) 

Close to the horizon the leading-order terms in Y and Y' are a2x'^ and 2a2X — these are 
used to apply the horizon boundary condition. Bardeen found that the equation in this 
form can now be stably integrated from a; = to a; = 00. The Bi^ and B^ut coefficients 
extracted can be used to obtain the grey-body factors using similar formulae to those 
discussed later for A^^ and Aout- 

For n = 1, equations ()4.19|1 and ()4.21|) - ()4.23|) are slightly modified and become 

cPR dR 
x(x + 2) , - 2icurh (x^ + 2a; + 1) — ^ - [-2iujrJx + 1) + A] = ; (4.24) 
aa;^ ax 

a, = l + t-^; (4.25) 
2c<jrh 

02 = -7- — r [(A + 2iuri,)ai - 2iujrh] ; (4.26) 

4(1 — -JCoTh) 

/ = 4(1 -icjrh)a2a;. (4.27) 



Unfortunately no analogous method has been found for n > 2 because there is then 
no simple equivalent of eq. ()4.19|) . Instead, an alternative transformation of the radial 
equation ()4.6|) was used. Writing y = r/ry^ and -Ri = y F{y) e~^'^^'* , the wave equation 
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becomes 

<PF . dF 



(hy^) ^ + 2y{h- lur^y) :i- - + 1)F = . (4.28) 
Since on the horizon y = 1 and h = 0, the boundary conditions become F{1) = 1 and 



d£ 
dy 



y=i 2a;rh 



^ (4.29) 



Using both this method and the Bardeen method for the n = 1 case provided a useful 
cross-check on the spin-1 numerical results. 

For fermions no special transformation was necessary and therefore the radial equa- 
tion ()4.6|) was used. However the application of the boundary condition at the horizon 
is made slightly easier by the transformation Pg = A'^Rs so that the asymptotic solution 
at the horizon becomes 

pf ) = e"*""* + 42 e'""* . (4.30) 
Since we require Aq^]. = 0, a suitable boundary condition to apply as r — >• fh is that 

Ps = l, (4.31) 



while using eq. (j4.9|l we also obtain 

dPs . dr* iu 
dr dr h{r) 

The above boundary conditions ensure that |^[n''P = 1- The asymptotic form for Pg at 



infinity now looks like 



since A — as r — oo. 



= ^ ^ + 4ut^ — (4.33) 
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The numerical work described in this chapter was almost entirely performed using 
the NDSolve package in Mathematica. However in some cases Fortran programs calling 
NAG routines (specifically D02EJF and its associated routines) were used as a check on 
the numerical robustness of the results obtained in Mathematica. The NAG routines 
implement the Backward Differentiation Formulae (BDF) whilst in Mathematica either 
the BDF method or the implicit Adams method is used. 

There are various considerations which must be taken into account in order to obtain 
results to the desired accuracy (at least three significant figures). Firstly, although the 
horizon boundary condition cannot be applied exactly at rh (due to singularities in the 
boundary condition and the differential equation) the error introduced by applying the 
condition at r = rn (where Arh = Tn — ''"h '"h) must be small. This can be investigated 
by studying changes in the grey-body factors for order-of- magnitude changes in Arh. 
Similarly it must be verified that the value Too used as an approximation for 'infinity' 
does not introduce errors which will affect the accuracy of the result. This is illustrated 
in Figure 14.11 in which a^^^ = Yli (^f^^ is plotted in units of vrr^. The oscillations in the 
numerical result die away as the asymptotic form given in eq. ()4.10|) is obtained. A value 
of Too = 10000 was used in the work that follows to ensure the grey-body factor was 
accurate to at least three significant figures. 

Care must also be taken that the numerical integration procedure is sufficiently ac- 
curate that significant integration errors are avoided even at large values of r. Plots like 
those in Figure 14.11 are also useful for investigating this because they can be used to 
confirm that there is no significant gradient on a straight line drawn through the centre 
of the oscillations. 

Finally, for each energy being considered, enough angular momentum modes must 
be included in the summation so that only modes which do not contribute significantly 
are neglected. For any value of i the absorption coefficient will approach unity at high 
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Figure 4.1: Stability of the n = 0, s = 0, ujr^ = 0.5 grey-body factor as r^o is 
changed. The right plot is a magnified version of part of the left plot to illustrate 
the accuracy obtained. 



enough energies and the (2£+l) factor in eq. ()4.1H] means that higher angular momentum 
modes will start to dominate the grey-body factor. It was found that for the highest 
values of curh considered in this work it was necessary to include contributions from in 
excess of ten angular momentum modes. 



4.3 Numerical results for brane emission 

In this section, results are presented for grey-body factors and emission rates for brane- 
localized scalar, fermion and gauge boson fields, as obtained by numerically solving the 
corresponding equations of motions. The relationship between the absorption coefficient 
Af^ and the A coefficients in equations ()4.H()|1 and ()4.3Hj) is different in each case, so 
each will be considered separately. 
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4.3.1 Spin fields 

The numerical integration of eq. ()4.fjj) for s = yields the solution for the radial function 
Ro{r) which smoothly interpolates between the asymptotic solutions of eqs. ()4.H()|1 and 
fl4.33|l in the near-horizon and far-field regimes respectively. The absorption coefficient 
is easily defined in terms of the in-going and out-going energy fiuxes at infinity, or equiv- 
alently by the corresponding wave amplitudes and A^^^) in the same asymptotic 

regime. Therefore eq. ()4.12j) may be written in the form [170] 

l4°)p = l-|7tf p = l- 

where TZe is the corresponding reflection coefficient. 

The plot presented in Figure 14.21 shows, for several values of n, the grey-body factor 
for the emission of scalar fields on the brane (for completeness, the plots include values 



A 



out 



A. 



(oo) 



(4.34) 
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of n ruled out on astrophysical grounds, i.e. n = 1 and 2).^ The grey-body factor is 
obtained by using eq. 1)4.111) and summing over the angular momentum number £. 

For n = and turh 0, the grey-body factor assumes a non-zero value which is equal 
to 47rr^ — that is, the grey-body factor for scalar fields with a very low energy is given 
exactly by the area of the black hole horizon. As the energy increases, the factor soon 
starts oscillating around the geometrical optics limit = 277rr^/4 which corresponds to 
the spectrum of a black body with an absorbing area of radius r,, = 3-\/3^h/2 [165-167]. 
The n = result agrees exactly with Page's result presented in Figure 1 of [75] (this 
is also found to be the case for the fermion and gauge boson grey-body factors shown 
later). 

In the extra-dimensional case, the grey-body factor has an asymptotic low-energy 
value which is the same for all values of n, and at high energies it again starts oscillating 
around a limiting value. This is always lower than the 4-dimensional geometrical optics 
limit because the effective radius depends on the dimensionality of the bulk space-time 
through the metric tensor of the projected space-time in which the particle moves. For 
arbitrary n, it adopts the value [163] 



The above quantity is a strictly decreasing function of n which causes the asymptotic 
grey-body factor, = vrr^ , to become more and more suppressed as the number of 
extra dimensions projected onto the brane increases. The values of o"g are tabulated in 
Table 113 for different values of n.^ 



^Throughout the numerical analyses the horizon radius is an arbitrary input parameter which 
remains fixed. 

^For the larger values of n it was found that these asymptotic values are not approached until 
relatively high energies, typically cjrh ^ n. 




(4.35) 
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n 





1 


2 


3 


4 


5 


6 


7 




6.75 


4 


3.07 


2.60 


2.31 


2.12 


1.98 


1.87 



Table 4.1: High-energy limits of grey-body factors for brane emission, given in 
units of 7rr^ . 

The power series expression of the grey-body factor determined in [170] matches 
the exact solution only in a very limited low-energy regime. In the limit ior\^ —>■ 0, 
the asymptotic value Anr^ is recovered as expected; however, as the energy increases 
the exact solution rapidly deviates from the behaviour dictated by the dominant term 
in the wrh expansion. This was confirmed in [171] where full analytic results for the 
grey-body factors were plotted, although still in the low-energy approximation. The 
behaviour depicted in Figure 3 of ref. [171] is much closer to the exact one, shown here 
in Figure 14. 2[ and successfully reproduces some of the qualitative features including 
the suppression of the grey-body factor as the dimensionality of the bulk space-time 
increases. Nevertheless, as previously stated, even this result rapidly breaks down as the 
energy increases; it is particularly unreliable for the n = case and fails to reproduce 
the high-energy oscillations for any value of n. This means that the exact numerical 
solution obtained in this work is the only reliable source of information concerning the 
form of the extra-dimensional grey-body factor throughout the energy regime. 

The numerical solution for the grey-body factor allows the computation of the energy 
emission rate for scalar fields on the brane. Using eq. ()4.14|) it is found that the suppres- 
sion of the grey-body factor with n does not lead to the suppression of the emission rate 
itself. The behaviour of the differential energy emission rate is given in Figure As n 
increases, the increase in the temperature of the black hole, and therefore in its emitting 
power, overcomes the decrease in the grey-body factor and leads to a substantial en- 
hancement of the energy emission rate. Figure 14.31 clearly shows that the enhancement 
of the peak of the emission curve can be up to several orders of magnitude compared to 
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n=0 
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Figure 4.3: Energy emission rates for scalar fields on the brane from a (4 + n)D 
black hole. 

the 4-dimensional case. In addition the peak in the spectrum is shifted to higher values 
of the energy parameter ujr^ since for fixed rh an increase in n corresponds to an increase 
in the temperature of the radiating body. 

In order to quantify the enhancement of the emission rate as the number of extra 
dimensions projected onto the brane increases, the total flux and power emissivities were 
computed, for a range of values of n, by integrating eqs. 1)4.11) and ()4.2p with respect to 
uj. The results obtained are displayed in Table ^21 The relevant emissivities for different 
values of n have been normalized in terms of those for n = 0. From the entries of this 
table, the order-of-magnitude enhancement of both the flux and power radiated by the 
black hole as n increases is again clear. 



4.3.2 Spin 1/2 fields 



Unlike the case of scalar fields, the study of the emission of fields with non-vanishing spin 
involves, in principle, fields with more than one component. Equation ()4.6)1 depends on 
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n 



Flux 



Power 




1 
2 
3 
4 
5 
6 
7 



1 

4.75 
13.0 
27.4 
49.3 
79.9 
121 
172 



1 

8.94 
36.0 
99.8 

222 
429 
749 
1220 



Table 4.2: Flux and power emissivities for scalar fields on the brane. 



the helicity number s which can be either +1/2 or— 1/2 and this leads to radial equations 
for the two different components of the field. As mentioned in ^4.2.2^ the upper (+1/2) 
and lower (—1/2) components carry mainly the in-going and out-going parts of the 
field respectively. Although knowledge of both components is necessary to construct 
the complete solution for the emitted field, the determination of either is sufficient to 
compute the absorption coefficient A^j^^ {j is the total angular momentum number). For 
example, if the in-going wave is known in the case of the emission of fields with spin 
s =1/2, eq. ()4.12|) may be directly written as [171,176] 



The above follows by defining the incoming flux of a fermionic field as the radial compo- 
nent of the conserved current, = V2(7'^^'^'^^^, integrated over a two-dimensional 
sphere and evaluated at both the horizon and infinity. 

The grey-body factor is again related to the absorption probability through eq. 1)4.111] 
with i now being replaced by the total angular momentum j. Numerically solving the 
radial equation ()4.6j) and computing A^^"^^ gives the behaviour of the grey-body factor, 
in terms of the energy parameter ur^ and the number of extra dimensions n. This is 
shown in Figure 14.41 for four different values of n. 




(4.36) 
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0.2 0.4 0.6 0.8 1 1.2 

Figure 4.4: Grey-body factors for fermion emission on the brane from a (4-|-n)D 
black hole. 

As in the scalar case, at low energies the grey-body factor assumes a n-zero asymptotic 
value; this depends on the dimensionality of space-time and increases with n. The 
enhancement of d'l^^^\uj) with n in the low-energy regime continues up to intermediate 
values of curh where the situation is reversed: as n becomes larger, the high-energy grey- 
body factor becomes more and more suppressed as was found in the scalar case. The 
full analytic results derived in [171] provide a reasonable description of the low-energy 
behaviour except in the n = case; however, as expected, they fail to give accurate 
information for the high-energy regime. As for scalar fields, the high-energy grey-body 
factors for fermions are shown to oscillate around asymptotic values determined by the 
effective radius of eq. ()4.H5|1 . 

The energy emission rate for fermion fields on the brane is shown in Figure 14.51 for 
various values of n. As n increases, the power emission is found to be significantly 
enhanced at both low and high energies. The emission curves show the same features 
as for the emission of scalar fields, i.e. the peak increases in height by several orders 
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of magnitude and shifts towards higher energies (due to the increasing black hole tem- 
perature). Quantitative results regarding the enhancement with n of both the flux and 
power fermion emission spectra were obtained by integrating the data in Figure |331 they 
are shown in Table 14.31 Once again, the enhancement is substantial, and in fact is even 
more significant than for scalar emission. 



n 


Flux 


Power 





1 


1 


1 


9.05 


14.2 


2 


27.6 


59.5 


3 


58.2 


162 


4 


103 


352 


5 


163 


664 


6 


240 


1140 


7 


335 


1830 



Table 4.3: Flux and power emissivities for fermions on the brane. 
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Figure 4.6: Grey-body factors for gauge boson emission on the brane from a 
(4 + n)D black hole. 



4.3.3 Spin 1 fields 

In the case of the emission of gauge boson fields, the incoming flux can be computed from 
the energy-momentum tensor T^*^ = 2(7^^,0"^^,^"^^ the time-radial component is 

integrated over a two-dimensional sphere and then evaluated at the horizon and at 
infinity. By making use of the solution for the in-going wave, the following expression 
for the absorption probability of eq. ()4.12|) . is obtained [171,176]: 



1^ 



(1)|2 



A 



(h) 



A 



(oo) 



(4.37) 



The exact results for the grey-body factors and energy emission rates for gauge boson 
fields are given in Figures 14.61 and 14.71 respectively. A distinct feature of the grey-body 
factor for gauge fields, previously observed in the 4-dimensional case, is that it vanishes 



92 



4.3 Numerical results for brane emission 



1-1 



n=6 



10' 




0.5 1.0 







1.5 



2.0 



Figure 4.7: Energy emission rates for gauge fields on tiie brane from a (4 + n)D 
black hole. 

when 0. The same behaviour is observed for all values of n. This leads to the 

suppression of the energy emission rate, in the low-energy regime, compared to those 
for scalar and fermion fields. Up to intermediate energies the grey-body factors exhibit 
enhancement with increasing n as in the case of fermion fields. An asymptotic behaviour, 
similar to the previous cases, is observed in the high-energy regime with each grey-body 
factor assuming, after oscillation, the geometrical optics value (which decreases with 
increasing n). This result establishes the existence of a universal behaviour of all types 
of particles emitted by black holes at high energies. This behaviour is independent of 
the particle spin but dependent on the number of extra dimensions projected onto the 
brane. 

Again the total flux and power emissivities for gauge fields on the brane can be 
obtained for different values of n. The exact results obtained by numerically integrating 
with respect to energy are given in Table I^Til As anticipated, the pattern of enhancement 
with n is also observed for the emission of gauge bosons. It is worth noting that the 
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n 


Flux 


Power 





1 


1 


1 


19.2 


27.1 


2 


80.6 


144 


3 


204 


441 


4 


403 


1020 


5 


689 


2000 


6 


1070 


3530 


7 


1560 


5740 



Table 4.4: Flux and power emissivities for gauge fields on the brane. 

relative enhancement observed in this case is the largest among all particle types — this is 
predominantly due to the large suppression of gauge boson emission in four dimensions. 

4.3.4 Relative emissivities for different species 

It is interesting to investigate how the relative numbers of scalars, fermions and gauge 
bosons emitted on the brane change as the number of extra dimensions projected onto 
the brane varies. In other words, this means finding out what type of particles the black 
hole prefers to emit, for different values of n, and what fraction of the total energy each 
particular type of particle carries away during emission. 

Comparing the energy emission rates for different types of particles and for fixed n, 
gives the qualitative behaviour. Figures [4. 8r a) and I4.8r b) show (with a linear scale) the 
a power spectra for n = and n = 6 respectively; these two figures demonstrate very 
clearly the effects already discussed, namely the orders-of-magnitude enhancement of 
the emission rates and the displacement of the peak to higher energies, as n increases. 
Figure |3IHt a) reveals that, in the absence of any extra dimensions, most of the energy 
of the black hole emitted on the brane is in the form of scalar particles; the next most 
important are the fermion fields, and less significant are the gauge bosons. As n increases. 
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0.2 0.4 0.6 0.8 1 1 2 3 4 5 6 



Figure 4.8: Energy emission rates for the emission of scalars, fermions and 
gauge bosons on tlie brane witii (a) n = and (b) n = 6. 



the emission rates for all species are enhanced but at different rates. Figure IHHtb) clearly 
shows that, for a large number of extra dimensions, the most effective 'channel' during 
the emission of brane-localized modes is that of gauge bosons; the scalar and fermion 
fields follow second and third respectively. The changes in the flux spectra are similar 
as n increases. 

The above behaviour can be more helpfully quantified by computing the relative 
emissivities for scalars, fermions and gauge bosons emitted on the brane by integrating 



n 


s = 





^ = 1 


s = 1 




n 


s = 


^ = 1 


s = 1 





1 




0.37 


0.11 







1 


0.55 


0.23 


1 


1 




0.70 


0.45 




1 


1 


0.87 


0.69 


2 


1 




0.77 


0.69 




2 


1 


0.91 


0.91 


3 


1 




0.78 


0.83 




3 


1 


0.89 


1.00 


4 


1 




0.76 


0.91 




4 


1 


0.87 


1.04 


5 


1 




0.74 


0.96 




5 


1 


0.85 


1.06 


6 


1 




0.73 


0.99 




6 


1 


0.84 


1.06 


7 


1 




0.71 


1.01 




7 


1 


0.82 


1.07 


Black body 


1 


0.75 


1 




Black body 


1 


0.87 


1 



Table 4.5: Flux emission ratios. Table 4.6: Power emission ratios. 
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the flux and power emission spectra. The relative emissivities obtained in this way are 
shown in Tables 1^31 and (they are normalized to the scalar values). The ratios for 
n > 1 are available for the flrst time in the literature as a result of this numerical work, 
while the n = results would appear to be the most accurate ones available. From [75], 
flux and power ratios are found to be 1 : 0.36 : 0.11 and 1 : 0.56 : 0.23 respectively, in good 
agreement with those ratios shown in Tables 14.51 and 14.61 In [31] the power ratio for 
= is given as 40 : 19 : 7.9 i.e. 1 : 0.48 : 0.20 which is in less good agreement. More 
careful examination shows that the relative power emitted by the s = 1/2 and s = 1 
degrees of freedom agrees exactly with Table 031 and so the disagreement seems to come 
from the scalar value. This is unsurprising since the scalar value quoted in [31] was only 
estimated from a plot in [166]. 

The entries in these tables reflect the qualitative behaviour illustrated above for some 
extreme values of the number of extra dimensions. For n = 0, the scalar flelds are the 
type of particle which are most commonly produced and the ones which carry away 
most of the energy of the black hole emitted on the brane; the fermion and gauge flelds 
carry approximately 1/2 and 1/4, respectively, of the energy emitted in scalar flelds, 
and their fluxes are only 1/3 and 1/10 of the scalar flux. For intermediate values of n, 
the fermion and gauge boson emissivities are considerably enhanced compared to the 
scalar one and have become of approximately the same magnitude — e.g. for n = 2 the 
amount of energy spent by the black hole in the emission of fermions and gauge bosons 
is exactly the same, although the net number of gauge bosons is still subdominant. For 
large values of n, the situation is reversed: the gauge bosons dominate both flux and 
power spectra, with the emission of fermions being the least effective channel both in 
terms of number of particles produced and energy emitted. The reader is reminded that 
the above results refer to the emission of individual scalar, fermionic or bosonic degrees 
of freedom and not to the elementary particles. However it is easy to combine the ratios 
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in Table 14.61 with the relevant numbers of degrees of freedom (see Table 13. 1|) to obtain 
relative emissivities for different particle types. 

These numerical results confirm that the relative emissivities of different particles 
produced by small, higher-dimensional black holes depend on the number of extra di- 
mensions projected onto the brane. Therefore, if Hawking radiation from such objects 
is detected, this could provide a way of determining the number of extra dimensions 
existing in nature. 



4.3.5 Total flux and power emitted 

Whilst it is often only necessary to know the relative emissivities of the fields of different 
spins, sometimes (for example, when estimating the black hole lifetime) the absolute 
values are required. For the case of brane emission from a non-rotating black hole, the 
flux emitted in a field of spin s can be written as 

io exp(^/TH)Tl27r J, exp (^) ^ 1 27r ' ' 
This expression is obtained by integrating equation ()4.13|1 using the definition 

r(^)(u;rh) = X;(2^+l)l'^f I'- (4-39) 

e 

Since the integrand in eq. ()4.38p is entirely a function of ur^, the substitution x = ury, 
can be used to express F^^^ as 

where the only dependence on the radius (or equivalently the mass) of the black hole is 
contained in the factor outside the integral. 
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n 






r _p(i) 





00133 


000486 


000148 


1 


00631 


00439 


00283 


2 


0.0173 


0.0134 


0.0119 


3 


0.0364 


0.0283 


0.0301 


4 


0.0655 


0.0499 


0.0596 


5 


0.106 


0.0789 


0.102 


6 


0.160 


0.116 


0.159 


7 


0.229 


0.163 


0.231 



Table 4.7: Total flux emitted in fields of different spins. 



n 


^2 p(0) 
n 


^2p(l/2) 


2p(l) 
n 





0.000298 


0.000164 


6.72x10-^ 


1 


0.00266 


0.00232 


0.00182 


2 


0.0107 


0.00973 


0.00971 


3 


0.0297 


0.0265 


0.0296 


4 


0.0661 


0.0575 


0.0686 


5 


0.128 


0.109 


0.135 


6 


0.223 


0.187 


0.237 


7 


0.362 


0.299 


0.386 



Table 4.8: Total power emitted in fields of different spins. 



A similar procedure, starting from equation ()4.2j) . shows the power emitted in a field 
of spin s to be 

rlJo exp (i^) T 1 27r • ^ ' ' 

These results of these integrations are given in Tables 14.71 and 14.81 The values in 
these tables can be combined with the total numbers of degrees of freedom for each 
spin. These are given in the last row of Table ItTD and are 4, 90 and 24 for s = 0, 1/2 
and 1 respectively. Therefore the total flux F and power P emitted by black holes of 
different numbers of dimensions can be calculated — the results are shown in Table 
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n 




rlP 





0526 


0.0175 


1 


0.489 


0.263 


2 


1.56 


1.15 


3 


3.41 


3.21 


4 


6.18 


7.09 


5 


9.98 


13.5 


6 


14.9 


23.4 


7 


21.1 


37.6 



Table 4.9: Total flux and power emitted from black holes of different numbers 
of dimensions. 



It is possible to use these values for the total flux emitted from extra-dimensional 
black holes to comment on the usual assumption that black hole decays can be considered 
as quasi- stationary. The simplest way of trying to verify this is to compare the typical 
time between emissions (given by F~^) with rh (the time for light to cross the black hole 
radius in the natural units being used). To be sure that the quasi-stationary approach 
is valid, we would require ^ rh or equivalently r^F <^ 1. From Table 031 this can 
be seen to be true only for the cases n = and n = 1. 

Another requirement for the validity of a semi-classical description of black hole 
production and decay is that the lifetime r ^ 1/Mbh so that the black hole is a well- 
defined resonance. This assumption can also be tested here by using the values of P in 
Table 1331 The relationship 



P=-^?§!1 = 4, (4.42) 

at r'^ 



h 



can be integrated as follows: 



JMbh ^^P(4+n) \ n + 2 j Jq V^P(4+n) 



2 

n + 1 



(4.43) 
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tMbh 



n 



Mbh = 5Mp(4+„) 



Mbh = 10Mp(4+n) 







47500 



761000 



1 



202 
23.3 
6.60 
2.77 
1.43 
0.846 
0.544 



1610 
148 
37.4 
14.6 
7.23 
4.12 
2.59 



2 



3 
4 
5 
6 
7 



Table 4.10: Values of tMbh for different values of n. 



In the above, the final Planck phase and any complications caused by kinematically 
forbidden emissions have been ignored. Proceeding with the integration gives 



Values of tMbh for different values of n are shown in Table 14.101 both in the case 
-^BH = 5Mp(4+„) and in the case Mbh = 10Mp(4+„). It should be noted that switching 
to convention 'a' for the definition of Mp(4+„) the high-n values in the table would be 
significantly altered (for example, the n = 7 values would be multiplied by a factor of 
13.2). The 'long'-lifetime requirement is clearly different to the quasi-stationary issue 
discussed above since it depends on the initial mass of the black hole. However given 
the limits on Mp(4_|_„) and the energy available at the LHC, it would again seem that 
this requirement may not be satisfied for the higher values of n. 

The results presented in Table 14.101 differ significantly from equivalent results which 
can be obtained from an expression given in [178]. Partly this is due to inclusion of the 
grey-body factors here, but this does not account for order-of-magnitude differences. It 
would appear that the lifetime is incorrect in [178] because of the failure to take all the 
degrees of freedom of the emitted particles into account. 




(4.44) 
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n 


(iV) 


s 


rSrl 




Mrh = 10Md^^_i_ ^ 
-'"^BH -'-^-'''^F(4+n) 


Muu = 5Afp/'/i_L^i 


M-RU = 10MpM4-^1 





75 3 


301 


314 


1260 


1 




36 1 


43.1 


122 


2 


6.82 


17.2 


21.0 


52.9 


3 


4.78 


11.4 


14.2 


33.7 


4 


3.80 


8.75 


11.0 


25.2 


5 


3.22 


7.23 


9.13 


20.5 


6 


2.80 


6.20 


7.92 


17.5 


7 


2.52 


5.50 


7.06 


15.4 



Table 4.11: Values of (N) and 5bh (calculated using eqs. (jTlB]) and (|XT3|l 

respectively) for different values of n. 



The numerical results presented in Table 1^711 can also be used to calculate the average 
particle multiplicity, again completely ignoring the complications of the Planck phase and 
kinematic constraints on the decay. Defining rh_F = f, the average particle multiplicity 
can be expressed as 



P Jo A/7rMp(4+„) \ n + 2 j Jq V^P(4+n) 

Performing the integration gives 



n + l 



(4.45) 



where the final identification has been made by referring back to equation ()3.15|) . Ta- 
ble 14.111 shows the average multiplicities obtained for different values of n along with 
the (initial) entropy S'bh- The average multiplicity values given improve on those in [34] 
(where (A^) ~ Mbh/2Th was assumed based on the constant temperature approxima- 
tion) and [179] (where a more careful approach was taken but still without the full 
numerical grey-body factors). 



Chapter 4: Grey-body Factors 



101 



The average muhiphcities in the first column clearly do not satisfy (A^) ^ 1 for 
n > 3 which means that a significant amount of the decay will be affected by kinematic 
and Planck phase considerations. However using the alternative definition of Mp(4+n) 
these results would again be modified (for n = 7 they would be multiplied by a factor of 
3.64). There is, of course, no fundamental difference between these conventions — setting 
-^p(4+n) = 1 TeV in convention 'a' simply corresponds to a lower value of the Planck 
mass in convention 'd' and hence for a fixed black hole mass the Hawking temperature 
is lower and the average multiplicity higher. 

The tabulated entropies allow comparison with the requirement that 1 / v^Sbh should 
be much less than unity if statistical fluctuations of the number of micro-canonical 
degrees of freedom are to be small (as is necessary for a semi-classical description of the 
black hole decay to be valid). 

4.4 Emission in the bulk 

An important question regarding the emission of particles by higher- dimensional black 
holes is how much of this energy is radiated onto the brane and how much is lost in 
the bulk. In the former case, the emitted particles are zero-mode gravitons and scalar 
fields as well as Standard Model fermions and gauge bosons, while in the latter case 
all emitted energy is in the form of massive Kaluza-Klein gravitons and, as discussed 
in ^4.11 possibly also scalar fields. In [163], it was shown that the whole tower of KK 
excitations of a given particle carries approximately the same amount of energy as a 
massless particle emitted on the brane. Combining this result with the observation that 
many more types of particles live on the brane than in the bulk, it was concluded that 
most of the energy of the black hole goes into brane modes. The results obtained in [163] 
were only approximate since the dependence of the grey-body factor on the energy of 
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the emitted particle was ignored and the geometric expression for the area of the horizon 
(valid in the high-energy regime) was used instead. 

In order to provide an accurate answer to the question of how much energy is emitted 
into the bulk compared to on the brane, it is necessary to take into account the depen- 
dence of the grey-body factor on both energy and number of extra dimensions. There 
is a similar discussion in [179] but it is incomplete because full numerical results for the 
grey-body factor in all energy regimes were not available. In this section the emission 
of scalar modes in the bulk is thoroughly investigated. Exact numerical results are pro- 
duced for the grey-body factors and energy emission rates in terms of the energy and 
number of extra dimensions. This allows the questions mentioned above to be addressed 
by calculating the total bulk-to-brane relative emissivities for different values of n. 

4.4.1 Grey-body factors and emission rates 

The analysis in this section is relevant for gravitons and possibly scalar fields, and 
requires knowledge of the solutions of the corresponding equations of motion in the 
bulk. Only the case of scalar fields is considered here since the bulk equation of motion 
is known — the emission of bulk scalar modes was previously studied analytically, in the 
low-energy regime, in [170].*^ 

A scalar field propagating in the background of a higher-dimensional, non-rotating 
(Schwarzschild-like) black hole, with line-element given by eq. ()3.10|) . satisfies the fol- 
lowing equation of motion [170] 



''After this work was completed, some plots of grey-body factors for bulk scalars were discovered in 
Figure 1 of [109]. It is unclear how these were calculated and only the £ — Q contributions for n =1-3 
are shown; however, they provided a useful cross-check on some of the results presented here. 




(4.47) 
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This radial equation was obtained by using a separable solution similar to equation ()4.5|) 
but with sSriO) replaced by S'^iOi, 62,..., On). The functions e*™^ S^{9i, 62,..., 6^) 
can be written as Yp{VL2+n) '■, these are generalizations of the usual spherical harmonic 
functions to the case of (3 + n) spatial dimensions [180] and account for the + n + 1) 
term in eq. ()4.47|1 . 

As for the emission of particles on the brane, the determination of the grey-body 
factor for bulk emission requires the above equation to be solved over the whole radial 
domain. The exact solution for the radial function must interpolate between the near- 
horizon and far-field asymptotic solutions, given by 

= (t'"^''-' + e*--* , (4.48) 



and 

R^-^ = ^ + 4~t^ ^ , (4.49) 

respectively. Again the boundary condition that no out-going solution should exist near 
the horizon of the black hole is imposed which means that A^^}^ = 0. The solution at 
infinity comprises, as usual, both in-going and out-going modes. 

The absorption probability \Ai\'^ may then be calculated either by using eq. ()4.34p or 
from the ratio | A^^'' p (this latter approach is found to be more numerically stable 
for low energies and/or larger values of n). A tilde will now denote bulk quantities, 
as opposed to brane quantities which carry a hat. The grey-body factor (Tf {uj) may be 
determined by using eq. ()4.3|1 . The dimensionality of the grey-body factor changes as 
the number of extra dimensions varies; therefore, in order to be able to compare its 
values for different n, it is normalized to the horizon area of a (4 -|- n)-dimensional black 
hole. The discussion of the normalization procedure and the subsequent comments on 
the geometrical optics limit are due to Kanti [2,177]. 
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First equation ()4.3p is re-written in the form 



2" 

ai{uj) = — r 

TT 



n 



A. 



(u;rh)" 



+2 



Np 



(2£ + n + 1) {^ + n)\ 
(n + 1)! 



and the horizon area in the bulk is defined as 



2-K n+l 



A=<+W d'^W j sin'^e^+i d(sin^„+i 



k=l 



n + 3 



(4.50) 



where is the multiplicity of states corresponding to the same partial wave i. For a 
(4 + n)-dimensional space-time it is given by 



(4.51) 



(4.52) 



Equation ()4.50p allows the computation of the grey-body factor's low-energy limit 
once the corresponding expression for the absorption coefficient is determined. Analytic 
results for Ai were derived in [170] by solving eq. ()4.47|1 in the two asymptotic regimes 
(near-horizon and far-field) and matching them in an intermediate zone. It was found 
that the low-energy expression of the £ = absorption coefficient has the form 



l^ol 



a;rh 



n+2 



r[(n + 3)/2]^ 



+ 



(4.53) 



where the ellipsis denotes higher-order terms in the uor^i power series. These terms, as 
well as the corresponding expressions for higher partial waves, vanish quickly in the limit 
coTh 0, leaving the above term as the dominant one. Substituting into eq. ()4.5()j) it is 
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clear that, in the low-energy regime, the grey-body factor is given by the area A-^ of the 
black hole horizon. This behaviour is similar to the 4-dimensional case except that the 
area of the horizon now changes with n. 

In the high-energy regime, it is anticipated that an equivalent of the geometrical 
optics limit discussed in ^4.3.11 will be recovered. In four dimensions, the low- and 
high-energy asymptotic limits are 47rr^ and vrr^ respectively. This has led to the naive 
generalization that, in an arbitrary number of dimensions, the high-energy expression 
for the grey-body factor will be approximately f2n+2r"+^/4, where is defined as in 
eq. ()1.8|) . It will shortly be shown that this is in fact an over-estimate of the high-energy 
limit. 

As in four dimensions it is assumed, for high energy particles, that the grey-body 
factor becomes equal to the area of an absorptive body of radius r^, projected on a plane 
parallel to the orbit of the moving particle [165]. According to ref. [163], the value of 
the effective radius remains the same for both bulk and brane particles and is given 
by eq. ()4.35|1 . The area of the absorptive body depends strongly on the dimensionality 
of space-time and its calculation requires one of the azimuthal angles to be set to tt/2. 
A careful calculation reveals that the 'projected' area is given by 

^'- (.T2)ri(;+2)/2| -""-;rT^ («^) 

The above relation reduces to the usual 4-dimensional result {Ap = 7rr^) for n = 
but, compared to the naive fi„+2r^+^/4, leads to values reduced by 50% for higher 
values of n. Assuming that the grey-body factor at high energies becomes equal to the 
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n 









1.69 


1.69 


1 


1.70 


2 


2 


1.77 


2.36 


3 


1.85 


2.72 


4 


1.93 


3.09 


5 


2.01 


3.45 


6 


2.08 


3.81 


7 


2.16 


4.17 



Table 4.12: High-energy limits of grey-body factors for bulk emission, given in 
units of the (4 + n)-dimensional area A\^. 



absorptive area Ap, it can be explicitly written as 
""^"^ + 2 a„+2 

(4.55) 

_ 1 r[(n + 3)/2] fn + ?>\ /n + 3\ ~ 
~ (n + 2) r[(n + 2)/2] J 

In the above the same normalization is used as in the low-energy regime — that is, the 
normalization is in terms of Ah, the area of the (4 + n)-dimensional horizon. The 
values predicted by eq. ()4.55|1 are tabulated in Table 14. 12^ along with the more naive 
VLn+2 r"+^/4 prediction. 

Turning now to the numerical analysis, the grey-body factors can be found by using 
eq. ()4.5()j) and the exact numerical results for the absorption coefficients. Their behaviour 
is shown in Figure 14.91 As it was anticipated after the above discussion, the normalized 
grey-body factors, in the low-energy regime, tend to unity for all values of n as each 
one adopts the value of the black hole horizon area to which it has been normalized. 
As for the emission of scalar fields on the brane, the grey-body factors are suppressed 
with increasing n in the low-energy regime, start oscillating at intermediate energies and 
then tend to their asymptotic high-energy limits. A simple numerical analysis shows that 
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Figure 4.9: Grey-body factors for scalar emission in the bulk from a (4 + n)D 
black hole. 



the naive expression f2„_|_2 r"+^/4 fails to describe the high-energy asymptotic limits for 
all values of n larger than zero. In contrast, eq. ()4.55j) gives asymptotic values which 
are verified by the numerical results. As in the brane emission case, it is found that 
relatively large values of toTh are required before the asymptotic values are approached 
for the higher values of n (this can be seen in Figure in which the n = 4 and 6 curves 
clearly do not attain the asymptotic values given in Table . 

In general, the suppression of the grey-body factor for bulk emission at low energies 
is milder than the one for brane emission. However this does not lead to higher emission 
rates for bulk modes compared to those for brane modes: the integration over the phase- 
space in eq. ()4.2|1 involves powers of turh which cause an increasingly suppressive effect 
in the low-energy regime as n increases. Nevertheless, the increase in the temperature 
of the black hole (still given by Th = (n + l)/47rrh) eventually overcomes the decrease 
in the grey-body factor and causes the enhancement of the emission rate with n at high 
energies (as well as shifting the spectrum peak to higher energies). The behaviour of 
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the differential energy emission rates as a function of the energy parameter curh is shown 
in Figure 14.101 for some indicative values of n. Comparing these numerical results with 
the analytic results of ref. [170] shows that the earlier results were reasonably successful 
in describing the low-energy behaviour of both the grey-body factors and the energy 
emission rates for scalar fields in the bulk. 



4.4.2 Bulk-to-brane relative emissivities 

The aim of this section is to perform an analysis which provides an answer to the question 
of the relative bulk-to-brane emissivity. This requires the differential energy emission 
rates in the bulk and on the brane to be evaluated; the two quantities are then compared 
for different numbers of extra dimensions. 
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Equation ()4.2|1 for energy emission in the bulk may ahernatively be written, in terms 
of the absorption coefficient, as 

at ^ exp (cu/Th) — 1 ztt 

The above can be compared with the corresponding expression for the emission of brane- 
locahzed modes given, as in eq. ()4.14|1 . by 

at ^ exp (cu/Th) — 1 27r 

where iV^ = 2£ + 1. Since the temperature is the same for both bulk and brane modes, 
the bulk-to-brane ratio of the two energy emission rates will be simply given by the 
expression 
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n 





1 


2 


3 


4 


5 


6 


7 


Bulk/Brane 


1 


0.40 


0.24 


0.22 


0.24 


0.33 


0.52 


0.93 



Table 4.13: Relative bulk-to-brane energy emission rates for scalar fields. 

and will depend on the scaling of the multiplicities of states and the absorption coeffi- 
cients with n.' 

As is clear from eq. ()4.5H) . the multiplicity of bulk modes A^^ increases quickly for 
increasing n, while Ni remains the same. However, it turns out that the enhancement 
with n of the absorption probability \A.i\^ for brane emission is considerably greater than 
the one for bulk emission. This leads to the dominance of the emission of brane-localized 
modes over bulk modes, particularly for intermediate values of n. The behaviour of this 
ratio is shown in Figure 14.111 This figure shows that, in the low-energy regime, the 
ratio for large n is suppressed by many orders of magnitude, compared to the value of 
unity for n = 0. In the high-energy regime on the other hand, the suppression becomes 
smaller and the ratio seems to approach unity. A more careful examination reveals that 
bulk modes in fact dominate over brane modes in a limited high-energy regime which 
becomes broader as n increases.-" 

A definite conclusion regarding the relative amount of energy emitted in the two 
'channels' — bulk and brane — can only be drawn if the corresponding total energy emis- 
sivities are computed. By integrating the areas under the bulk and brane energy emission 
rate curves, the relative energy emission rates are determined. The results obtained, for 
values of n from 1 to 7, are given in Table 

'The absorption coefficients are related to the grey-body factors through eq. H4.3|l . This includes 
a multiplicative coefficient which depends on both wrh and n, and so the coefficients might have a 
completely different behaviour from the grey-body factors themselves. 

■Tor the very high energies for which this is true, it is not clear that the expressions for higher- 
dimensional Hawking radiation will remain valid. 
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From the entries of the Table 14.131 it becomes clear that the emission of brane- 
localized scalar modes is dominant, in terms of the energy emitted, for all the values 
of n considered. As n increases, the ratio of bulk-to-brane emission gradually becomes 
smaller, and is particularly suppressed for intermediate values, i.e. n = 2-5; in these 
cases, the total energy emitted in a bulk mode varies between approximately 1/3 and 
1/4 of that emitted in a brane mode. As n increases further, the high-energy dominance 
of the bulk modes mentioned above gives a boost to the value of the bulk-to-brane 
ratio — nevertheless, the energy ratio never exceeds unity. 

The above analysis provides exact results for the energy emission rates of brane and 
bulk scalar modes and gives considerable support to earlier, more heuristic, arguments 
[163], according to which a (4 + n)-dimensional black hole emits mainly brane modes. 
A complete confirmation would mean performing a similar analysis for the emission of 
gravitons, but there are still theoretical and numerical issues which prevent this. 

4.5 Rotating black holes 

We have already seen that a black hole created in a hadron collider is expected to have 
some angular momentum J about an axis perpendicular to the plane of parton collision. 
Using the same naive argument as in ^3.2.11 we assume J = bM-Q^L/'i (where b is the 
impact parameter of the colliding partons) and that b < 2rh if a black hole is to form. 
The maximum possible value of the rotation parameter defined in eq. ()3.7j) is then found 
to be 



max 



n + 2 



(4.59) 



2 



The equivalent of eqs. ()4.13|1 and ()4.14|1 are now 




(4.60) 
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and 

where the Hawking temperature is now given by 



(n+l) + (n-l)a^ 



and r2 is defined by 



(1 + a2)rh 



For the rotating case, the metric takes on a more comphcated form than previously. 
In general (4 + n) -dimensional objects arc described by (n + 3)/2 angular momentum 
parameters; however, it is assumed that here there is only one non-zero parameter (about 
an axis in the brane). This is reasonable because the partons which collide to produce 
the black hole are themselves on the brane. Hence the metric reduces to 



^2 (. \ , 2a//sin^g S . 



v^j/)2 / 2 I 2 I (l^lJ^Sm^O 1 •2/172 

— ijdB — r -I- a H — — sm Bdif , 



(4.64) 



where 

A = r^ + a^ - and S = + cos^ 9, (4.65) 

with a = a*rh and /x = rg"*"^ (rs is the Schwarzschild radius of a non-rotating black hole 
of the same mass) . 

Again the field equation is separable using a solution of the same form as in the 
non-rotating case, but now the radial equation is 

^~^Tr + (^^^^ + 4i.c.r + .A" - .A^) i?. = , (4.66) 
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where 

K ={r'^ + a^)uj - am and = .Ef + 2s + a^u^ - 2amuj . (4.67) 

The angular equation is now 



' sin6' + 



sin Ode V de 
2ms cot 9 m^ 



+ a^u"^ cos^ 6 - 2aujs cos ^ + s - cot^ + ^ET sTPie) = , 
sm sin ^ / 

(4.68) 

where e^"^'^ gT^^ (6) are known as spin- weighted spheroidal harmonics. 

The black hole horizon is given by solving A(r) = 0. Unlike the 4D Kerr black 
hole for which where are inner and outer solutions for rh, for n > 1 there is only one 
solution of this equation (as can easily be seen graphically — Figure Ii.l2|) . In the n = 
case the maximum possible value of a* is 1 otherwise there are no solutions of A = 
(this is forbidden as it would mean there is no horizon and hence there is a naked ring 
singularity at r = 0). For n > 1 there is no fundamental upper bound on a^ but only the 
bound given in eq. ()4.59p which was argued geometrically. For general n, the horizon 
radius rh is found to be 

rh = "^^^ . (4.69) 

(l + a2)n+i 



The solution of the radial equation to obtain the grey-body factor is carried out in 
exactly the same same way as for the non-rotating black hole. There is a slight difference 
in the asymptotic solution close to the horizon in that the u in eq. ()4.H()j) is replaced by 
k where 
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Figure 4.12: Number of solutions of A(r) = for n = 0, 1 and 4. On all three 
plots, the two functions shown are /(r) = r^ + a^ and g(r^ = fir^~"- for particular 
values of a and /i. 



and that the equivalent of eq. ()4.9|) is now 



dr* + 
dr A(r) 



(4.71) 



4.5.1 Calculation of the eigenvalue 

The only other complication is to determine from the angular equation. The 6 wave 
functions ^T™ are no longer like those of the spin-weighted spherical harmonics and so 
sE"^ 7^ sXi . The eigenvalues of the spin-weighted spheroidal harmonics are functions 
of auj and can be obtained by a continuation method [181]. This is a generalization 
of perturbation theory which can be applied for arbitrarily large changes in the initial 
Hamiltonian for which the eigenvalues are known. Making the aw-dependence of both 
T and E clear, eq. ()4.fj8|l can be written as 



{Ho + Hi) sTJ' 



, auo 



.E'^iau) sTr{e,aiu) 



(4.72) 
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where 

1 d f „ d \ ( 2ms cot 9 rn? ^ ^ \ , . „ . 

^0 = —^Ta ^i^^;^ + -T-n + ^ - s'^oee] , 4.73 

sm dO \ dv ) \ sm a* sm Q J 

and 

Til = (aW cos^ 9 - 2aujs cos 0) . (4.74) 
From these equations it is evident that 

,Tr(^,O) = ,Sr(0), (4.75) 

and that 

,£;r(o) = = ^{^ + 1) - s{s + 1) . (4.76) 

The continuation method was outhned in [182] although this paper unfortunately 
seems to contain a number of errors. First consider the case in which auj is small and 
normal perturbation theory can be used. This means we have 

sE]r{auj) = sXe - {sim\ni\sim) + . . . , (4.77) 



and 



where 



{s£'m\ni\sim) = J dn {,3"^)* ,8"^ Hi . (4.79) 



Care must be taken about the signs of the terms in eqs. ()4.77|1 and ()4.78|1 due to the 
sign of the eigenvalue in eq. ()4.72j) . This seems to account for some differing results 
in [182]. 
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The product in eq. ()4.79|) can be split into terms for which the following standard 
formulae are useful: 



1 2 / 2/ + 1 

(s£'m| cos^ e\s^m) = -Sw + - 77^7^ ) (£2m0|£'m)(^2(-s)0|£'(-s)), (4.80) 



and 



{s^'m\ cos9\slm) 



2i + l 



{£lmO\em){il(-s)0\e(-s)), (4.81) 



2i' + l^ 

where {lil2^ifT^2\LM) are Clebsch-Gordan coefficients. Using the above we obtain 



771m 



Xe - 2au^ + 0[{au) 



2m? +l~2i{i+l) 
(2£-l)(2^-3) 



if s ^ 



(4.82) 



Again this differs (in the s = case) from the result given in [182]. However it agrees 
with [183] where several discrepancies in the literature are helpfully clarified. 

To employ the continuation method we write the sT^{9^auj) functions in the basis 
of the 6'-parts of the spin-weighted spherical harmonics: 



(4.83) 



By differentiating eq. ()4.72j) and applying the same techniques as in perturbation theory 
it is possible to obtain the results 



d{aLu) 



{a 



(4.84) 



o,/3 



and 



d{auj) 



- E 



{a\f3)sB^e' 



(4.85) 
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10.5 




auj 

Figure 4.13: Comparison of the numerical result obtained in this work for the 
s = 1,^ = 3, m = —2 eigenvalue (solid line) with the polynomial approximation 
in [79] (dashed line). The lines are exactly on top of each other (the maximum 
difference between them is ~ 10~^). 

where = {sam\dl-ii/ d{auj)\s(3m) . The initial conditions are obtained from eqs. 

K7^ and K7^ which imply that ^5™,(0) = 5u,. By integrating eqs. (lOl) and (OKI) 
it is possible to obtain the eigenvalues of the spin-weighted spheroidal functions for 
any s, £ and m and for arbitrarily large values of au. This integration was performed 
by using the NDSolve package in Mathematica to implement a Runge-Kutta method. 
The computational effort involved increases significantly for larger values of £ as it is 
necessary to include a larger number of terms in the summations, increasing the number 
of differential equations to be solved. 

Although only the s = eigenvalues were required here, as a check the continuation 
method was also used to calculate the s = 1 eigenvalues as a function of auo. This allowed 
comparison with the polynomial approximations given in [79]. Excellent agreement was 
found — see Figure E. 131 for an example. 
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4.5.2 Calculation of the grey-body factor 

In order to calculate the grey-body factors in the rotating case, it is necessary to know 
the equivalent of eqs. ()4.H4|1 . ()4.H(jj) and ()4.H7|1 . Unfortunately, for fermions and bosons 
these expressions are not available in the literature. However for the scalar case, since 
eq. ()4.34j) only involves coefficients at infinity, exactly the same formula holds even when 
the black hole is rotating. This means that for scalars it has been possible to numerically 
calculate the grey-body factors for rotating black holes. 

Since the denominator in eqs. ()4.fi()|l and ()4.fiH) has m-dependence, it is less useful 
to plot an equivalent of the a used in the non-rotating case,"^ so only the power spectra 
are shown. They are presented as a function of curh for different values of a* and hence 
cannot be used to directly compare the spectra from two black holes of the same mass 
(but different J) since this is itself a function of Mbh- 

Figure I?. 141 can be directly compared to Figure 2 in [145] and shows markedly different 
behaviour for values of ur^i larger than ~ 0.2. Almost all of this effect is due to the 
low-energy expansion used by those authors to obtain the grey-body factors. Their 
approximation that ^-E™ = ^A^ = -|- 1) — s{s + 1) is much less significant except at 
the very highest values of curh. 

These numerical results for the grey-body factors also allowed the statement in [31] 
that emission is dominated by modes with £ = m to be tested. It was found that a plot 
like Figure 14.141 looks the same at the ~ 90% level if only the the i = m modes are 
included in the sum of eq. ()4.6Hl 

''The absorption cross section (Tabs summed over all angular momentum modes can still be calculated; 
the high-energy asymptotic value appears to increase with a,, apparently in contradiction with [156] 
which states that roughly the same value is expected as in the non-rotating case. 
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4.5.3 Super-radiance 

The numerical resuhs obtained allow the possibility of super-radiance (i.e. negative grey- 
body factors) to be investigated in an extra-dimensional context. Previously an ana- 
lytic approach [113] was used to confirm super-radiance for bulk scalars incident on 5D 
black holes, but there do not seem to be any other results in the literature for higher- 
dimensional black holes. To study this phenomenon it is useful to consider angular 
momentum modes individually and, for comparison with the 4D scalar results presented 
in [78], the quantity on the horizontal axis of the plots is 



rather than the curb used in previous plots. The vertical axis is the absorption probability 
(in fact — l^f.mP) expressed as a percentage so that it gives the percentage energy 
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Figure 4.15: Super-radiant scattering of scalars by maximally rotating (a* = 1) 
4-dimensional black holes. 

amplification of an incident wave. Figure I?. 151 shows excellent agreement with the results 
in Figure 1 of [78] which are for n = and a^, = 1 (equivalent to a = Mrh in the 4D case). 
Figure 14.161 is for n = 2 and 6 but still with the same value of the rotation parameter. 
It is clear that the peak amplification is more significant than previously and becomes 
comparable to that for gauge bosons in the 4D case. It is still found that the i = m = 1 
mode provides the greatest amplification. However Figure 14.171 shows that for larger 
values of the maximum amplification can occur in modes with larger values of i. This 
trend continues as n and a* increase further: for n = Q, equation ()4.59|1 shows that the 
geometric maximal rotation is a™'^^ = 4.0 and for this value of the rotation parameter 
the maximum energy amplification is found to be just under 9% (and occurs in the 
i = m = 7 mode). 

As in the 4-dimensional case, the region in which the grey-body factors are negative 
exactly corresponds to the region in which the thermal factor in the denominator of 
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equations like ()4.60p and ()4.61|) is negative; hence the quantum emission rate always 
remains positive. For all angular momentum modes the grey-body factor is positive for 
uj > mQ and the change from negative to positive aX u = mQ occurs such that the flux 
and energy spectra are smooth positive functions. Therefore although super-radiance 
can in principle be used to extract energy from rotating black holes, its significance when 
considering the Hawking radiation spectra is limited. 

4.6 Conclusions 

The revival of the idea of extra space-like dimensions in nature has led to the formulation 
of theories which allow gravity to become strong at a significantly lower energy scale. 
This has opened the way for the proposal of the creation of miniature higher-dimensional 
black holes during collisions of energetic particles in the earth's atmosphere or at ground- 
based particle colliders. The Hawking radiation emitted by such black holes is described 
by the relevant grey-body factors; a wide variety of these factors have been calculated 
numerically in the work described in this chapter and, at the same time, some open 
questions from previous analyses in the literature have been addressed. 

The exact numerical results obtained in ^4.31 for the grey-body factors and emission 
rates of scalars, fermions and gauge bosons confirm that previous analytic results [170, 
171] were only valid in limited low-energy regimes. The grey-body factors depend on 
both the dimensionality of space-time and the spin. They adopt the same high-energy 
asymptotic values (which decrease as n increases) for all particle species, confirming the 
geometrical optics limit of [163]. In the low-energy limit the grey-body factors are found 
to be enhanced as n increases for s = 1 and 1/2 and suppressed for s = 0. 

Using these grey-body factors it was possible to accurately calculate the flux and 
power emission spectra for the black holes. As the number of extra dimensions projected 
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on the brane increases there is a substantial (orders of magnitude) enhancement in the 
energy emission rate, mainly as a result of the increase in black hole temperature with 
n. However the effect of including the grey-body factors is that the increase in the rate 
depends on the spin of the particle studied. The computed relative emissivities show 
that scalar fields, which are the dominant degree of freedom emitted by n = black 
holes, are just outnumbered by gauge bosons for large values of with the fermions 
becoming the least effective emission channel. Therefore the emission spectra and the 
relative particle emissivities may possibly both lead to the determination of the number 
of extra dimensions. 

The emission spectra for modes of different spins were combined with the total num- 
ber of degrees of freedom for all the particle types which could be emitted from an 
extra-dimensional black hole. This made it possible to calculate the total flux and 
power emission and hence also the expected time for, and average multiplicity of, the 
black hole decay. These calculations suggest that some of the assumptions usually made 
when modelling black hole decay by sequential emission of Hawking radiation may not 
be valid for the larger values of n. 

In section |^3J the details of the emission of bulk scalar modes were investigated with 
the aim being to provide an accurate estimate for the amount of energy lost into the bulk. 
Comparing the total bulk and brane emission rates for scalar modes, integrated over the 
whole energy regime, the conclusion is that more of the black hole energy is emitted in 
brane modes than bulk modes, for all values of n considered. The total emissivity in the 
bulk is less than 1/4 of that on the brane for n = 2-4, while it becomes substantial for 
the extreme value of n = 7 without, however, exceeding the brane value. The accurate 
results presented in this chapter provide firm support to the heuristic arguments made 
in ref. [163]. 
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No results are presented in this work concerning the emission of gravitons either on 
the brane or in the bulk. The derivation of a consistent equation, which can describe 
the motion of gravitons in the induced black hole background on the brane, is still under 
investigation. Nevertheless it is expected that, as for the emission of fields with spin 
s = 0, 1/2 and 1 the graviton emission becomes enhanced as the dimensionality of 
space-time increases, but remains subdominant compared to the other species at least 
for small values of n. 

The work on rotating black holes in ^4.51 suggests that, at least for the emission 
of scalars, the power spectra can be significantly modified compared to those in the 
Schwarzschild-like case. The spectra are dominated by the angular momentum modes 
with i = m and as the rotation parameter a^, increases, the peak in the spectrum shifts 
to much larger values of ur^. The phenomenon of super-radiance was also observed for 
higher- dimensional rotating black holes. The peak energy amplification was found to 
be significantly larger than in the four-dimensional case and occurred in higher angular 
momentum modes for larger values of the rotation parameter. 

A final comment is appropriate here concerning the validity of the results obtained 
in this work. As pointed out in the text, the horizon of the black hole is an input 
parameter of the analysis. The results are applicable for all values of which are 
smaller than the size of the extra dimensions R, no matter how small or large R is 
(as long as it remains considerably larger than ip to avoid quantum corrections). The 
analysis therefore remains valid for all theories postulating the existence of fiat extra 
dimensions and a fundamental scale of gravity even a few orders of magnitude lower 
than Mp. If the emission of Hawking radiation from these black holes is successfully 
detected, either at next-generation colliders or in the more distant future, the distinctive 
features discussed here may help in the determination of the number of extra dimensions 
existing in nature. Chapter El describes a black hole event generator which incorporates 
these grey-body effects and would be a useful tool in such an analysis. 



Chapter 5 

CHARYBDIS: A Black Hole Event 
Generator 



5.1 Introduction 

The work on black holes in this thesis was motivated by extra dimension models which 
allow the fundamental Planck scale to be of order a TeV. These models are attractive 
because they solve the hierarchy problem, but they are also exciting because they mean 
that gravity is strong at these energy scales and so black holes can be investigated at 
the LHC and other next-generation machines. If these models are realized in nature, we 
would expect particle accelerators at TeV-scale energies to be able to produce miniature 
black holes. These would then decay rapidly by Hawking evaporation, giving rise to 
characteristic high-multiplicity final states. 

There has already been much discussion of this issue in the literature (for example, in 
many of the works cited in Chapters El and 0} but little work has been done in trying to 
realistically simulate black holes at the LHC. This chapter describes the implementation 
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of a simple model of black hole production and decay which can be interfaced to existing 
Monte Carlo programs using the Les Houches accord [184]. The major new theoretical 
input to the generator is the inclusion of the results from Chapter |3] — the grey-body 
factors for black holes in extra dimensions [2]. The recoil and change of temperature of 
the black hole during decay is also taken into account, and various models are provided 
for the termination of the decay process. 

Section reviews important aspects of the theory of extra-dimensional black holes. 
Then, in ^5.31 there is a description of the CHARYBDIS'^ event generator itself including 
a discussion of the theoretical assumptions involved, details of the Fortran code and 
instructions on how to use the program. P. Richardson was involved in the develop- 
ment of some of this code, specifically in ensuring baryon number and charge conserva- 
tion, defining the colour flow required by general-purpose event generators, and making 
CHARYBDIS compatible with the Les Houches accord. An event display and some sam- 
ple production cross sections, particle spectra and particle emissivities produced using 
the event generator are presented in §5.41 Details of some preliminary experimental 
studies are given in ^5.51 and there is a concluding discussion in ^5.61 



5.2 Black hole production and decay 

Most of the details of the production and decay of black holes in extra dimension models 
have already been discussed in Chapter El Here some of the key features of the theory 
are re-iterated and some comments made on those which are particularly relevant for 
the task of event generation. 

A fundamental Planck scale as low as ~ TeV means that it is possible for tiny black 
holes to be produced at the LHC when two partons pass within the horizon radius set 



^'In Greek mythology, Charybdis was a nymph daughter of Poseidon and Gaia who was turned into 
a monster by Zeus. She hved in a cave at one side of the Strait of Messina (opposite the monster Scylla) 
and sucked water in and out three times each day. 
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by their centre-of-mass energy, Using the results of [161], the horizon radius for a 



non-spinning black hole is found to be 



A/7rMp(4+„) \Mp(4+„) 




) 



n + 1 




) 



n+1 



(5.1) 



where Mbh is the mass of the black hole and, as throughout the rest of this chapter, 
Mp(4+„) is the fundamental (4 + n)-dimensional Planck scale in convention 'd' of Ta- 



there is enough energy to be in the trans-Planckian regime, the parton-level cross section 
is Fnirr^ where F„ ~ 1. 

Black holes can be produced with any gauge and spin quantum numbers so to de- 
termine the p-p or p-p production cross section it is necessary to sum over all possible 
quark and gluon pairings. Although the parton-level cross sections grow with black hole 
mass, the parton distribution functions (PDFs) fall rapidly at high energies and so the 
cross section also falls off quickly. 

These miniature black holes are expected to decay instantaneously on LHC detector 
time scales (typical lifetimes are ~ 10^^^ s). The decay is made up of four phases 
(discussed in more detail in Chapter IH]): the balding phase, the spin-down phase, the 
Schwarzschild phase (accounting for the greatest proportion of the mass loss in four 
dimensions [156]) and the Planck phase. 

A black hole of a particular mass is characterized by a Hawking temperature and 
as the decay progresses the black hole mass falls and the temperature rises. For an 
uncharged, non-rotating black hole the decay spectrum is described, as in eq. ()4.1H|1 . by 



ble 11.21 Numerical calculations, discussed in §3.2.11 show that Mbh < V^- Provided 
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where the Hawking temperature is given by 

^ n + 1 , , 

Tu = , 5.3 

47rrs 

and is the absorption probability from which the grey-body factor can be obtained 

(see eq. ()4.3|) in Chapter H}. Equation ()5.2|) can be used to determine the decay spectrum 
for a particular particle, by considering the number of degrees of freedom (Table 13.11 
shows the relevant numbers for all the particle types which the black hole can emit). It 
is relevant for predicting both the spectra and the relative emission probabilities of the 
different particle types. 

The grey-body factors modify the spectra of emitted particles from that of a perfect 
thermal black body [73]; their spin dependence means that they are also necessary to 
determine the relative emissivities of different particle types from a black hole. Until 
the work in Chapter EJ the full numerical grey-body factors were only available in the 
literature for the 4D case [156,166]. 

It is usually assumed that a quasi- stationary approach to the black hole decay is 
valid — that is, after each particle emission the black hole has time to come into equilib- 
rium at the new temperature before the next particle is emitted.^ 

It was argued in [163] that the majority of energy in Hawking radiation is emitted 
into modes on the brane (i.e. as Standard Model particles) but that a small amount is 
also emitted into modes in the bulk (i.e. as gravitons). The work described in Chapter |3] 
provides support for this statement based on the relative emissivities of scalars in the 
bulk and on the brane. 

'^Note that the results of Chapter 01 have cast doubt on this assumption for cases with n > 3. 
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5.3 Event generator 

5.3.1 Features and assumptions of the event generator 

There are a number of features of the CHARYBDIS generator which, within the un- 
certainties of much of the theory, allow rehable simulation of black hole events. Most 
notable is that, unlike other generators (e.g. [178]), the grey-body effects are fully in- 
cluded. The dependence on both spin and the dimensionality of space-time means that 
grey-body factors must be taken into account in any attempt to determine the number 
of extra dimensions by studying the energy spectra and relative emission probabilities 
of particles from a black hole. When studying black hole decay, other experimental 
variables may also be sensitive to these grey-body effects. 

The generator also has an option to allow the black hole temperature to vary as the 
decay progresses and is designed for simulations with either p-p or p-p. 

CHARYBDIS only attempts to model in detail the Hawking evaporation phase which 
is expected to account for the majority of the mass loss. To provide a further sim- 
plification only non-spinning black holes are modelled. This is perhaps a less good 
approximation but comparison with the 4D situation suggests that most of the angular 
momentum will be lost in a relatively short spin-down phase [185]. The balding phase is 
difficult to model and is neglected — this is equivalent to the assumption that Mbh = 
in spite of the evidence that this will not necessarily be the case. A related assumption 
is that the cross section calculation assumes F„ of equation ()3.6|) is equal to unity (i.e. 
the parton-level cross section is assumed to be vrrg). 

It is hoped that the way in which the black hole decay is terminated will provide a 
reasonable approximation to the Planck phase of the decay. The generator also assumes 
that energy lost from the black hole in graviton emission can be neglected, based on the 
work in [163] and in Chapter |3] 
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As discussed in ^3.2.41 it is possible that black hole decay does not conserve baryon 
number. However the treatment of processes which do not conserve baryon number is 
complicated, in both the QCD evolution and hadronization, and has only been studied 
for a few specific processes [186-189]. At the same time, the violation of baryon number 
is extremely difficult to detect experimentally and so the effect of including baryon num- 
ber violation is not expected to be experimentally observable. Therefore CHARYBDIS 
conserves baryon number in black hole production and decay; lepton number, however, 
is not conserved. 

5.3.2 General description 

The black hole event generator developed attempts to model the theory as outlined in 
^5.21 within the assumptions of ^5.3.11 There are several related parameters and switches 
which can be set in the first part of the Les Houches subroutine UPINIT [184]. No other 
part of the charybdislOOO.F code should be modified. 

Firstly the properties of the beam particles must be specified. IDBMUP(l) and 
IDBMUP(2) are their Particle Data Group (PDG) codes (only protons and anti-protons 
are allowed) and the corresponding energies are EBMUP(l) and EBMUP(2). 

As already discussed, the geometric parton-level cross section (cr = Trrg) is used but 
the parameters MINMSS and MAXMSS allow the mass range for the black holes produced to 
be specified. This means that it is possible to adjust the lower mass limit at which this 
cross section for the trans-Planckian production of black holes is thought to become valid. 
The details of the Monte Carlo (MC) cross section calculation are given in Appendix O 

Three other parameters which must be set before using the event generator are 
TOTDIM, MPLNCK and MSSDEF. The total number of dimensions in the model being used 
is given by TOTDIM (this must be set between 6 and 11). There are a number of different 
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MSSDEF 


Conversion 


Convention 


1 


MPLNCK = (2"-V"-i)^Mp(4+„) 


a 


2 


MPLNCK = Mp(4+„) 


d 


3 


MPLNCK = (2"-37r"-i)^Mp(4+„) 


b 



Table 5.1: Definitions of the Planck mass. 



definitions of the Planck mass (set using MPLNCK) but the parameter MSSDEF can be set to 
three different values to aUow easy interchange between the three conventions outlined 
in Appendix A of [31]. The conversions between these conventions are summarized in 
Table IKTT] (in fact options 1,2 and 3 correspond to MPLNCK being set in the conventions 
a, d and b of Table ITT^ respectively) . 

It has been suggested that since black hole formation is a non-perturbative process, 
the momentum scale for evaluating the PDFs should be the inverse Schwarzschild radius 
rather than the black hole mass. The switch GTSCA should be set to . TRUE . for the first 
of these options and .FALSE, for the second.'^ The PDFs to be used are set using the 
Les Houches parameters PDFGUP and PDFSUF. 

As discussed in ^5.21 the Hawking temperature of the black hole will increase as the 
decay progresses so that later emissions will typically be of higher energy. However to 
allow comparison with other work which has ignored this effect, the TIMVAR switch can 
be used to set the time- variation of the Hawking temperature as on (.TRUE.) or off 
(.FALSE.). The indication in ^4.31 that a quasi-stationary approach to the black hole 
decay might not be valid for n > 3 means that the . FALSE . option might in fact give a 
better description of the decay for higher values of n. 

The emission probabilities of different types of particles are set according to the new 
theoretical results of Chapter HI (see particularly Table . Heavy particle production 

■^As confirmed in [190], tlic cross sections quoted in reference [31] were actually calculated with the 
latter PDF scale. 
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is allowed and can be controlled by setting the value of the MSSDEC parameter to 2 for 
top quark, W and Z production, or 3 to include the Higgs boson as well (MSSDEC=1 gives 
only light particles) . The production spectra for heavy particles may be unreliable if the 
initial Hawking temperature is below the rest mass of the particle in question. 

If GRYBDY is set as . TRUE . the particle types and energies are chosen according to the 
grey-body modified emission probabilities and spectra. If instead the .FALSE, option 
is selected, the black-body emission probabilities and spectra are used. The choice of 
energy is made in the rest frame of the black hole prior to the emission, which is isotropic 
in this frame. As overall charge must be conserved, when a charged particle is to be 
emitted the particle or anti-particle is chosen such that the magnitude of the black hole 
charge decreases. This reproduces some of the features of the charge-dependent emission 
spectra in [191] whilst at the same time making it easier for the event generator to ensure 
that charge is conserved for the full decay. 

Although the Planck phase at the end of decay cannot be well modelled as it is not 
well understood, the Monte Carlo event generator must have some way of terminating 
the decay. There are two different possibihties for this, each with a range of options for 
the terminal multiplicity. 

If KINCUT= . TRUE . termination occurs when the chosen energy for an emitted particle 
is ruled out by the kinematics of a two-body decay. At this point an isotropic NBODY 
decay is performed on the black hole remnant where NBODY can be set between 2 and 5. 
The NBODY particles are chosen according to the same probabilities used for the first part 
of the decay. The selection is then accepted if charge and baryon number are conserved, 
otherwise a new set of particles is picked for the NBODY decay. If this does not succeed in 
conserving charge and baryon number after NHTRY attempts the whole decay is rejected 
and a new one generated. If the whole decay process fails for MHTRY attempts then the 
initial black hole state is rejected and a new one generated. 
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In the ahernative termination of the decay (KINCUT=. FALSE.) particles are emitted 
according to their energy spectra until Mbh falls below MPLNCK; then an NBODY decay as 
described above is performed. Any chosen energies which are kinematically forbidden 
are simply discarded. 

In order to perform the parton evolution and hadronization the general-purpose event 
generators require a colour flow to be defined. This colour flow is defined in the large 
number of colours (A^^c) limit in which a quark can be considered as a colour line, an 
anti-quark as an anti-colour line and a gluon both a colour and anti-colour line. A simple 
algorithm is used to connect all the lines into a consistent colour flow. This algorithm 
starts with a colour line (from either a quark or a gluon) and then randomly connects 
this line with one of the unconnected anti-colour lines (from either a gluon or an anti- 
quark). If the selected partner is a gluon the procedure is repeated to find the partner 
for its colour line; if it is an anti-quark one of the other unconnected quark colour lines is 
selected. If the starting particle was a gluon the colour line of the last parton is connected 
to the anti-colour line of the gluon. Whilst there is no deep physical motivation for this 
algorithm it at least ensures that all the particles are colour-connected and the showering 
generator can proceed to evolve and hadronize the event. 

After the black hole decay, parton-level information is written into the Les Houches 
common block HEPEUP to enable a general-purpose event generator to fragment all emit- 
ted coloured particles into hadron jets, and perform the decays of any unstable particles. 

5.3.3 Control switches, constants and options 

Those parameters discussed in ^5. 3. 21 which are designed to be set by the user are sum- 
marized in Table 
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Name 


Description 


Values 


Default 


IDBMUP(2) 


PDG codes of beam particles 


±2212 


2212 


EBMUP(2) 


Energies of beam particles (GeV) 




7000.0 


PDFGUP(2) 


PDFLIB codes for PDF author group 




-1 


PDFSUP(2) 


PDFLIB codes for PDF set 




-1 


MINMSS 


Minimum mass of black holes (GeV) 


< MAXMSS 


5000 


MAXMSS 


IVTaximum mass of black holes I GeV i 


<^ c m erierffv 


c m erierffv 


MPLNCK 


Planck mass (GeV) 


< MINMSS 


1000.0 


MSSDEF 


Convention for MPLNCK (see Table 15. 1|) 


1-3 


2 


TOTDIM 


Total number of dimensions (4 + ti) 


6-11 


6 


GTSCA 


Use r^^as the PDF momentum scale 


LOGICAL 


. FALSE . 




rather than the black hole mass 






TIMVAR 


Allow Th to change with time 


LOGICAL 


. TRUE . 


MSSDEC 


Choice of decay products 


1-3 


3 


GRYBDY 


Include grey-body effects 


LOGICAL 


. TRUE . 


KINCUT 


Use a kinematic cut-off on the decay 


LOGICAL 


.FALSE. 


NEDDY 


Number of particles in remnant decay 


2-5 


2 



Table 5.2: List of parameters with brief descriptions, allowed values and default 
settings. 



5.3.4 Using charybdislOOO.F 

The generator itself only performs the production and parton-level decay of the black 
hole. It is interfaced, via the Les Houches accord, to either HERWIG [62,63] or PYTHIA 
[192] to perform the parton shower evolution, hadronization and particle decays. This 
means that it is also necessary to have a Les Houches accord compliant version of either 
HERWIG or PYTHIA with both the dummy Les Houches routines (UPINIT and UPEVNT) 
and the dummy PDFLIB subroutines (PDFSET and STRUCTM) deleted. For HERWIG, 
the first Les Houches compliant version is HERWIG6.500 [193]; for PYTHIA, version 
6.220 [194] or above is required.*^ 

"^Versions of PYTHIA above 6.200 support the Les Houches accord but cannot handle more than 7 
out-going particles, which is necessary in black hole decays. 
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The black hole code is available'^ as a gzipped tar file which includes 

• charybdislOOO.F (code for the black hole generator), 

• dummy.F (dummy routines needed if not using PDFLIB), 

• mainpythia.f (example main program for PYTHIA), 

• mainherwig.f (example main program for HERWIG), 

• charybdislOOO.inc (include file for the black hole generator). 

The Makefile must specify which general-purpose event generator is to be used (i.e. 
GENERATOR=HERWIG or GENERATOR=PYTHIA) and also whether PDFLIB is to be used 
(PDFLIB=PDFLIB if required, otherwise PDFLIB= ). The name of the HERWIG or PYTHIA 
source and the location of the PDFLIB library must also be included. 

If the code is extracted to be run separately then the following should be taken into 
account: 

• charybdislOOO.F will produce the HERWIG version by default when compiled; the 
flag -DPYTHIA should be added if the PYTHIA version is required. 

• dummy.F will by default produce the version for use without PDFLIB; the flag 
-DPDFLIB should be added if PDFLIB is being used. 

5.3.5 List of subroutines 

Table EISl contains a list of all the subroutines of the generator along with their functions. 
Those labelled by HW/PY are HERWIG/PYTHIA dependent and are pre-processed ac- 
cording to the GENERATOR flag in the Makefile. Many of the utility routines are identical 
to routines which appear in the HERWIG program. 

Ihttp : //www . ippp . dur . ac . uk/montecarlo/leshouches/generators/ cha rybdis/| 
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Name 


Description 




Les Houches routines 


UPINIT 


Initialization routine 


UPEVNT 


Event routine 




Particle decays 


CHDFIV 


Generates a five-body decay 


CHDFOR 


Generates a four-body decay 


CHDTHR 


Generates a three-body decay 


CHDTWO 


Generates a two-body decay 




Hard sub-process and related routines 


CHEVNT 


Main routine for black hole hard sub-process 


CHFCHG 


Returns charge of a SM particle 


CHFMAS 


Returns mass of a SM particle (HW/PY) 


CHHBHl 


Chooses next particle type if MSSDEC=1 


CHHBH2 


Chooses next particle type if MSSDEC=2 


CHHBH3 


Chooses next particle type if MSSDEC=3 


CHPDF 


Calculates the PDFs (HW/PY) 




Random number generators 


CHRAZM 


Randomly rotates a 2-vector 


CHRGEN 


Random number generator (HW/PY) 


CHRLOG 


Random logical 


CHRUNI 


Random number: uniform 




Miscellaneous utihties 


CHUBHS 


Chooses particle energy from spectrum 


CHULB4 


Boost: rest frame to lab, no masses assumed 


CHULOB 


Lorcntz transformation: rest frame lab 


CHUMAS 


Puts mass in 5th component of vector 


CHUPCM 


Ccntrc-of-mass momentum 


CHUROB 


Rotation by inverse of matrix R 


CHUROT 


Rotation by matrix R 


CHUSQR 


Square root with sign retention 


CHUTAB 


Interpolates in a table 




Vector manipulation 


CHVDIF 


Vector difference 


CHVEQU 


Vector equality 


CHVSUM 


Vector sum 



Table 5.3: List of CHARYBDIS subroutines with brief descriptions. 
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MPLNCK 


MINMSS 


TOTDIM 


Cross section (fb) 


1000.0 


5000 


8 


155 X 10^ 


1000.0 


7000.0 


8 


6.0 X 10^ 


1000.0 


10000.0 


8 


6.9 


1000.0 


5000.0 


10 


0.233 X 10^ 


1000.0 


7000.0 


10 


8.7 X 10^ 


1000.0 


10000.0 


10 


9.7 


2000.0 


2000.0 


11 


0.511 X 10^ 


6000.0 


6000.0 


7 


0.12 



Table 5.4: Black hole production cross sections for various different param- 
eters. Those in the upper part of the table can be directly compared with 
values quoted in [31] (they are calculated using MSSDEF=1, PDFGUP(I)=4 and 
PDFSUP(I)=48) whereas those in the lower part can be compared with [34] (they 
are calculated using MSSDEF=2, PDFGUP(I)=3 and PDFSUP(I)=34). In all cases 
MAXMSS=14000.0. 



5.4 Event generator output 



5.4.1 Sample cross sections 



In this section some indicative black hole production cross sections are shown for various 
different settings of the CHARYBDIS parameters. These are chosen so as to allow direct 
comparison with cross sections quoted in [31] and [34]. The results in Table IF!^ are 
found to be in good agreement with the previously published values. 



5.4.2 Sample event display 

Figure I3?T] is an example black hole event display produced from the CHARYBDIS output 
using the ATLANTIS program^ which is being developed for the ATLAS experiment. The 
black hole event used was generated with an earlier version of CHARYBDIS which did 
not fully take into account the grey-body effects; however it can still be considered as a 



jhttp : //atlant is . web . cern . ch/atlantis/| 



138 
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Figure 5.1: Sample event display (details in text). 
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'typical' black hole event which displays many of the expected features. The event shown 
is for a black hole mass of ~ 8 TeV in a scenario with T0TDIM=10 and MPLNCK=1000 . 
(MSSDEF=l). The event includes 9 quarks, a gluon, an electron, a positron and three 
neutrinos. 

5.4.3 Sample particle spectra and emissivities 

Figures I5.2H5.4I show the results, at parton level, of neglecting the time-variation of 
the black hole temperature (TIMVAR= . FALSE . , dashed line) or the grey-body factors 
(GRYBDY=. FALSE., dot-dashed line) for initial black hole masses in the range from 
MINMSS=5000.0 to MAXMSS=5500 . (with the defauh values for the other parameters). 
The solid line is for simulations with the default parameter settings (but with the same 
reduced range of initial black hole masses used in the other two cases). 

The effect of time-variation is to harden the spectra of all particle species. However, 
the effect of the grey-body factors depends on the spin, in this case slightly softening 
the spectra of scalars and fermions but hardening the spectrum of gauge bosons. 

Relative numbers of primary emitted particles were also calculated and are shown as 
percentages in Table The same parameters were used as for the spectra, except that 
in both cases the time- variation of the temperature was included (TIMVAR=.TRUE.). For 
comparison, the theoretical values used by the event generator (and easily calculable 
using Tables 13.11 and 14. 5|) are also shown. The differences are due to any kinematic 
constraints on the massive particles and also the imposed conservation of charge and 
baryon number; however since the average particle multiplicity per event is relatively 
high in these examples, the decay of the black hole is not significantly constrained. 
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5000 



I 4000 - 
> 




500 1000 1500 2000 2500 

Energy of Higgs boson (GeV) 

Figure 5.2: Parton-level energy spectra of Higgs bosons, ttih = 115 GeV. Solid: 
predicted energy spectrum of Higgs bosons from decay of black holes with initial 
masses 5.0-5.5 TeV. Dashed: neglecting time-variation of temperature. Dot- 
dashed: neglecting grey-body factors. 



12000 



10000 - 




500 1000 1500 2000 2500 

Energy of electron (GeV) 

Figure 5.3: Parton-level energy spectra of electrons and positrons. As Fig- 
ure Ell but for electron and positron spectra. 
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6000 - 




500 1000 1500 2000 2500 

Energy of photon (GeV) 

Figure 5.4: Parton-level energy spectra of photons. As Figure 15.21 but for 
photon spectra. 



Particle type 


Particle emissivity (%) 


GRYBDY= . 


TRUE. 


GRYBDY= . FALSE . 


Generator 


Theory 


Generator 


Theory 


Quarks 


63.9 


61.8 


58.2 


56.5 


Gluons 


11.7 


12.2 


16.9 


16.8 


Charged leptons 


9.4 


10.3 


8.4 


9.4 


Neutrinos 


5.1 


5.2 


4.6 


4.7 


Photon 


1.5 


1.5 


2.1 


2.1 


ZO 


2.6 


2.6 


3.1 


3.1 


W+ and W" 


4.7 


5.3 


5.7 


6.3 


Higgs boson 


1.1 


1.1 


1.0 


1.1 



Table 5.5: Comparison of relative emissivities of primary partons with their 
theoretical values. 
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5.5 Experimental studies 

Results of a full experimental study of black hole production and decay at the LHC will 
be presented in [4]. However some of the experimental issues will be highlighted here 
with reference to the theoretical work of Chapter 0] 

If extra dimensions do exist, then a priority will be to try to determine how many 
there are. The black hole generator described in this chapter was constructed with the 
aim of being able to more realistically examine ways of doing this experimentally. The 
grey-body factors of the previous chapter were incorporated for the same reason. 

One method for determining n is described in [34]: Black hole events are isolated 
using the experimental cuts mentioned in ^3.2.3[ events are put into 500 GeV bins based 
on the estimated black hole mass and then, for each mass bin, photon or electron energy 
spectra are produced. If Hawking temperatures can be determined by fitting the spectra, 
then the correlation between Th and black hole mass can be used to try to determine 
the number of extra dimensions (and potentially also the value of Mp(4+„)). 

It was argued in [34] that this method is quite successful but a number of unjustified 
approximations were made. These include a failure to take account of the temperature 
increase as the decay progresses, the grey-body effects described in Chapter lU and the 
secondary photons and electrons. Isolation cuts can be used to reduce the problem of 
secondary particles, but using the CHARYBDIS generator (with GETJET providing a 
very simple detector simulation) it was found that the time-dependence of the Hawking 
temperature significantly reduces the viability of such a method. The effect of the black 
hole not always being at rest is also found to have a significant effect on the energy 
spectra (although this is itself a matter of some uncertainty since the way in which a 
black hole recoils against high-energy emissions may be complicated [109, 111, 195]). 
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Attempts were made to perform a more sophisticated version of this method by using 
a more reahstic spectrum, obtained by numerically integrating the Hawking spectrum, 
for the fitting procedure. This approach was significantly more successful for n = 2 but 
less promising for larger values of n. This is because, as the emitted particle multiplicity 
decreases, it is difficult to know what cut-off to use when calculating the time-integrated 
spectrum, particularly without making assumptions about the value of Mp(4_|_„). This 
highlights a problem which was found repeatedly when dealing with the larger values of 
n: the expected particle multiplicity at the LHC is not particularly high (see calculations 
in §4.3.5j) and so it is difficult to extract information in a way which is not affected by the 
kinematic constraints on the decay, and the way in which the Planck phase is modelled. 

Experimentalists from the ATLAS collaboration have subsequently spent a signifi- 
cant amount of time trying to find event shape variables and other distributions which 
can be used to extract the number of extra dimensions from black hole events; however 
results are as yet inconclusive. It has hoped that by using the different event generator 
options for modelling the Planck phase it will be possibly to confirm that any proposed 
method is relatively insensitive to this final phase of the decay. However most of meth- 
ods investigated have been limited by the problem discussed above — that is, the low 
multiplicity at the LHC for large numbers of extra dimensions (and with Mp(4_|.„) as- 
sumed to be ~ TeV). A higher energy collider in the more distant future would alleviate 
many of these problems. The rate of fractional change of n decreases with n which also 
makes it progressively more difficult to distinguish n from n + 1 extra dimensions. It 
is unfortunate that these problems are most significant for large values of n since these 
have the best theoretical motivation and are the least constrained astrophysically. 
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5.6 Discussion 

The CHARYBDIS program described in this chapter would appear to be the most so- 
phisticated black hole event generator available. Figures f5.2H5.4l and Table IH7H1 confirm 
that the inclusion of the grey-body factors can have a significant effect on the particle 
energy spectra and emissivities. 

However it will have become clear to the reader that there are many assumptions and 
approximations made in the way the production and decay of black holes is modelled. 
Various improvements could be made, but most of them would be impossible without 
further theoretical work. 

A major advance would be to correctly model stages of the black hole decay other than 
the Schwarzschild phase. The balding phase could be taken into account if the amount 
of energy trapped by the black hole horizon was more accurately known. It was recently 
pointed out [151] that this might have a significant effect on the black hole discovery 
potential of the LHC Further numerical calculations to estimate the typical amount of 
trapped energy could also provide information about the expected distribution of initial 
angular momenta. This information, together with more theoretical work on the grey- 
body factors for rotating black holes, would allow the spin-down phase of the decay to 
be simulated in the event generator. If attempting to model the spin-down phase, the 
assumption of isotropic particle emission from the black hole would no longer be a good 
one. Grey-body factors are calculated for individual values of I and m and so it might 
be possible for the angular distribution of emitted particles to be modelled correctly 
(although to do this fully would also require the spin-weighted spheroidal harmonics of 
equation ()4.68p to be calculated). 

The assumption that the emission of energy into the bulk as gravitons (and possibly 
scalars) can be neglected is not ideal. The emission of gravitons from the black hole 
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could be taken into account but a full implementation of this would again require more 
theoretical work on the grey-body factors. 

Other possible improvements would be to more carefully account for the modifications 
in the energy spectra of massive gauge bosons and other heavy particles, or to allow the 
black holes to emit supersymmetric particles (since extra dimension models motivated 
from string theory are expected to be supersymmetric). 

One relatively simple modification to CHARYBDIS would be to try to simulate black 
holes in a separated fermions model [24, 25] which would avoid some of the proton decay 
issues (see discussion in ^3.2.4j) . A model like this would be implemented by altering the 
number of degrees of freedom available for decays as described in [35]. However, since the 
values of Mp(4_|_„) required in such models make black hole production experimentally 
inaccessible at the next generation of colliders, there is little motivation for such a 
modification. 

It would be desirable to have a better understanding of both the validity of the quasi- 
stationary approach to the decay of higher-dimensional black holes, and of the Planck 
phase of the decay. In the latter case, it seems unlikely that this will be possible without 
a full theory of quantum gravity and in the meantime the best approach would seem to 
be to try to extract information which is independent of the way in which the Planck 
phase is modelled. 



Chapter 6 
Conclusions 



The individual chapters of this thesis contain their own conclusions, so here there is 
only a brief overview, and a few comments about the future of the field. Two different 
possibilities for new physics beyond the Standard Model have been examined: supersym- 
metry and extra dimensions. Although the phenomenology of both models is diverse the 
particular aspects addressed were the detection of heavy leptons in intermediate scale 
supersymmetric models, and the production and decay of black holes in extra dimension 
scenarios. In both cases it was shown that the next generation of particle colliders will 
be able to investigate these models. 

The study of the detection of charged heavy leptons at the LHC using a time-of- 
flight technique showed that, in the particular model considered, the discovery mass 
range extends up to 950 GeV. In addition it was found that it would be difficult to use 
their differing angular distributions to distinguish the heavy leptons from scalar leptons 
if the lepton mass was greater than 580 GeV. 

The majority of the thesis addressed various aspects of the production and decay 
of miniature black holes. The grey-body factors which were calculated numerically in 
Chapter m are useful in that they allow a more realistic Monte Carlo event generator to 
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be constructed, as described in Chapter El Whilst the experimental signatures produced 
by black hole decays are expected to make events unmissable, it is important to model 
grey-body effects correctly if information (like the number of extra dimensions) is to 
be accurately extracted. It was found that the grey-body factors modify the particle 
energy spectra and emissivities and so might be important in such a task. The work on 
grey-body factors also allowed some of the usual assumptions and approximations made 
in black hole decay to be examined more carefully. The assumption that the relative 
emission of particles into the bulk is small was shown to be correct, but some doubt 
was cast on the usual approximation that a higher- dimensional black hole decay can be 
treated in a quasi-stationary manner. 

One conclusion from this work was that if the number of extra dimensions is more 
than three or four, the centre-of-mass energy available at the LHC is barely sufficient 
for the semi-classical approximations used to be legitimate. Therefore it might need a 
higher energy collider to investigate black hole production in detail. The LHC would 
still be expected to provide the first indications of extra dimensions through graviton 
production, Planckian effects and the first signs of black hole production. 

The question of how likely it is that any particular model for new physics is realized 
in nature is inevitably impossible to quantify. Every physicist has their own prejudices 
and in general supersymmetric models seem to be the favoured option. Extra dimension 
models are 'wacky' enough to attract the attention of the popular science media and are 
quickly dismissed by some physicists. However, as outlined in Chapter ^ they can be 
well motivated theoretically. 

Almost all high-energy physicists are convinced that the LHC will discover new 
physics of some kind. Leaving aside all the theoretical arguments, it would be un- 
precedented in the history of the field if the order-of-magnitude increase in available 
energy did not reveal something new. It may well be that nature surprises us and re- 
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veals physics that nobody has so far suggested. Even if extra dimensions do exist the 
compactification scheme may make the physics much more complicated than the simple 
cases discussed here. 

With half a decade still to wait before CERN's Large Hadron Collider has produced 
a substantial amount of data, there is good reason for theorists to continue investigating 
these models, and refining and extending their phenomenological studies. 



Appendix A 

Cross Section for Heavy Lepton 
Pair Production 



A.l 2^2 body scattering in centre-of-mass frame 



The differential cross section for 2 — > 2 scattering in tlie centre-of-mass frame is given 

by 

'M^fil >■ (A.l) 



In the above expression, the asterisk indicates the centre-of-mass frame, subscripts i and 
f refer to initial- and final-state particles, and s is the square of the total centre-of-mass 
energy. 

Therefore the major work in calculating the cross section is to determine the spin 
averaged square of the transition matrix element ( |Map), evaluated in the centre-of- 
mass frame. 
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A. 2 Definitions and conventions 



A. 2 Definitions and conventions 



The 4-momenta of the particles involved are defined in this section. Subscripts 1, 2, 3 
and 4 refer to the quark, anti-quark, lepton and anti-lepton respectively, and the angle 
6* is between the out-going lepton and the incoming quark in the centre-of-mass frame. 
Hence the {E,p) 4-momenta are as follows (note that the incoming quarks are treated 
as massless): 

Pi = (£;,o,o,£;); 

P2^{E,Q,Q,-E)- 

(A.2) 

Pa = (£;3,p3sin6l*,0,p3cos6l*) ; 
P4 = (-^4, -pa sin 6*, 0, -pa cos 0*) . 

Later we will require the Lorentz dot products of several pairs of these 4-vectors. 
The useful ones are listed below: 

Pi.p2 = 2E^ ; 

Pi.pa = £^£^3 — Ep3 cos 9* ; 

P2.p4 = EEi - Eps cos e* ; (A.3) 
Pi.p4 = EE4 + Eps cos 9* ; 
P2.pa = EE3 + Eps cos 9*. 
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A. 3 Matrix element calculation 



Using the Feynman rules, the matrix element for Z boson exchange is 



iMl 



v{p2)- - igzl^]^ (cV - 47^) u{pi) 



-1'9^lv 



s — my + iirtzTz 



, (A.4) 



which gives 



Mi 



9l 



s — m| + imzTz 



^{P2)l^- (cV - ck7') u{pi] 



U{P3)Y^ (Cy - 47^) V{p4 



. (A.5) 



Note that the matrix element for 7-exchange is easily obtained from this by setting mz 
and ca to zero, cy to 2, and g'^ to q^q^e^. 

Taking into account Z- and 7-exchange the full matrix element may be written as 
Mfi = M| + Ml Hence 



|M|p + \Ml\^ + 23fJ{ (M|)*M^} 



(A.6) 



We start by calculating the first term in this expression (from which the second term is 
easily obtained using the substitutions mentioned above). Finally the third term which 
contains the Z-7 interference information is calculated. 

To obtain the spin-averaged square of the transition matrix element which appears 
in equation ()A.1|) . we must average over the spins of the incoming quarks and sum over 
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A. 3 Matrix element calculation 



the spins of out-going leptons such that 



(A.7) 



spins 



A. 3.1 Z- exchange 



Firstly we want to calculate X^gping l-ZV^I This is given by 



E 



Xs-mlY + mlTl 

X Vr'{pA) (Cy + Ca7^) y^i^'bs) Us{pi) (cV + Ca7^) Y^^r(p2) ^^'^^ 

- Iv ■ ■ 1 _ 7^ f f 

X ^r{P2)^ (cV - Ca7^) li5(Pl)J li5'(P3)y (cv - Ca7^) ^^r'(P4)- 



spins 



To obtain the above expression the following have been used: 

{uYv)* = (MV7"t;)T = v^YViiyu = t;^(7°7V07°« = ; 
(uYl^v)* = (mV7Vv)^ = ^;t(75)t(7-)t(/)t„ ^ ^t^5(^o^.^O)^o^ (A.9) 
= ti^7^7°7''ii = -t'^7°7^7''ii = -v^Y'^ ■ 

We now make use of familiar trace identities {uJPaisUjS = Tr (uuT) etc. with a, [5 
Dirac matrix indices) as well as the following spinor identities: 




spins 



(A.IO) 




spins 
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Using these identities (and temporarily ignoring the propagators) eq. ()A.8|) becomes 



>7 ,f f 5^ 



Tr 



7.Pi (cV + ck7') ^l-P2^{c'v-cW) 



X Tr 



7 / 
(7-P4 - (cy + Ca7^) y(7-P3 + "^3) Y (cy - 47^) 



(All) 



We now evaluate the traces using the properties that the trace of the product of an 
odd number of gamma matrices is zero, Tr (7'^7'^) = 4:g'^'^, Tr (^z^^''^-'^^'^) = A{g^^g^^ — 
gf^Xgua^gf^aguX-^^ r^.^ ^ 5 ^ ^ ^ ^ ^ r^.^ ^ 5 ^ ^ ^ ^ A ^ ^ ^ ^ 4ie^''^'". Therefore the second 

trace becomes 



(cy)^ Tr (7.^47^7-^37'' - m3m4'j''Y) 
- (4)^ Tr (7.P47^7''7-P37^ - m3m47^7''7''7^) + 244Tr (7.P47^7''7-P37' 



(A12) 



(4) ' + (4) ' (p>3 + plp'3 - p^-p^r) - (4) ' - (4) 



24ci^e"^^'^(p4)a(p3)/3 



where the properties of 7^ (it squares to 1 and anti-commutes with the other 7-matrices) 
have also been used. The first trace is similar, but simpler as we are treating the incoming 
quarks as massless. It becomes 



[cyY + {c\Y [{pi)^,{p2)u + (pi)z.(P2)m -Pi-P2^Mz.] - 2cyc\ie^^suPlpi ■ (A.13) 
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A. 3 Matrix element calculation 



Finally we must multiply these two terms together, noting which terms do not con- 
tribute because of the anti-symmetry/symmetry of the indices. We obtain 

(4)' + (4)1 f(cV)' + (4)1 2 [ (Pl.P4)(P3.P2) + (Pl.P3)(P2.P4) ] 



+ 



(cy) - (ci) (cV) + (ck) '^rn3m4:Pl.p2-4:CyC^j^CyC'j^ej^Sueaixl3uPlpiip4)aip3)l3- 

(A.14) 



To complete the calculation we use the identity ep,,^a/3e^i/7(5 = '^{9i9p~9Z9p)^ which means 
that 



^'1ti5ueatil3uPlP2{Pi)a{Pz)l3 = ^tiual3^tiv'ySP\P2{P^)a{Pi) (3 

- H9i9} - 9^) pUipMPsh 
= 2 [ {pi.P3){p2-Pa) - {Pi-Pa){P2-P?,) 



(A.15) 



Putting all of the above together gives us that 



ml?) 



9i 



4 (s - m|)2 + m|r| 

X ([(4)' + (4)'] [(cV)' + (4)'] 2 [(P1.P4)(P3.P2) + (P1.P3)(P2.P4); 



+ 



(4)' - (4)' (4)^ + (4) 



2m3m4Pi.p2 



+ 



8Cv444 [ (Pl-P4)(P2-P3) - (Pl-P3)(P2-P4) ]) ■ 



A. 3. 2 7- exchange 



Using the replacements outlined above, it is immediately obvious that 



(A.16) 



mi?)- 



4 s2 



-32 [ (P1.P4)(P3-P2) + (Pl-P3)(P2-P4) + m3m4Pi.p2) ] ■ (A.17) 
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A. 3. 3 Interference term 

The final part of tlie squared matrix element to be calculated is K{(Mg )*Mg}. Neglect- 
ing the propagators for now, the spin sum is 



E 



r,r',s,s' 



= Tr 



7-Pi (cV + ck7^) y7-P27i^ 



X Tr 



(7.P4 - m4) (cy + Ca7^) y(7-P3 + m3)Y 



(A.18) 



Now the second trace is 



[cyTr {j.p4j^j.p2Y - m^rriAYY) + c^Tr (7-P47^7''7-P37'')] 

= 24 {p/ipl + p^p^ - P3.P45'^" - ^37714^'^'^) - 24ze-'^/^^(p4)a(P3)/3 , (A.19) 



and the first 



2cV [ {.Pl)^,{.P2)y + (Pi)i/(P2)m - Pi-P29i,A - 2c\ie^^suPlpi 



(A.20) 



Multiplying these together and using the same identity as before for the e-tensor we 
obtain 



2 2 f i 



4 s{s — m| — imzFz) 



X 8(cyCy [ (P1.P4)(P3-P2) + (Pl-P3)(P2-P4) + ^37714^1.^2] (^-21) 
+cick[(Pl-P4)(p2.P3) - (Pl-P3)(P2.P4)]). 
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A. 4 Differential cross section 



A. 3. 4 Complete matrix element squared 



Making the change of variables 



cv9z 
2e ' 



(A.22) 



and remembering that the real part of the interference term contributes twice, we obtain 



8e^ 

(|Mfip) = — [Ci [{pi.Pi){p3-P2) + {Pl-P3){P2-Pi)] + C2m3miPi.p2 



+C3 [ (P1.P4)(P2.P3) - (Pl.P3)(P2.P4) ] , (A.23) 



where 



Co 



(s-m|) +m|r| 



Z) "T "^z-^ z 

f J{ ^2 



C, = 2 



q^q'd\d\s {s - m|) 
(s-m|)Vm|r| ' (s-m|)Vm|r| 



+ 



(A.24) 

(A.25) 
(A.26) 



A. 4 Differential cross section 



Substituting the squared matrix element into equation (jA.lj) we have 



da 1 p3 



dn* 87r2 s3 E 



Ci [ {pi-Pi){P3-P2) + {Pi-P3){P2-Pa) ] + C2m3m4pi.p2 



+C3[{pi.Pi){p2.p3) - (P1.P3)(P2.P4)] . (A.27) 
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From ^A.2I we see that 



{Pi-Pa){P2-P3) = E'^lE^Ei + {Es + E^)p3 cosO* + pj cos^ 6*' 



(A.28) 



{PI.P3){P2.PA) = E''[EsE^ - {Es + E^)ps cos e* +pI cos^e*' 



and so, noting that s = AE^ and E^ + E4 = 2E, 



da e* 1 p3 



dVL* 167r2 s2 E 



[Ci^E^Ei + pI cos^ e*) + C2m3m4 + C^{2Ep^ cos 9*) ] . (A.29) 



The final resuh for the parton-level cross section (indicated by the 'hat') includes a 
factor of 1/3 since the quarks must both be of the same colour if they are to annihilate. 
So we have 



(qq ^ L-L+) = [C^ [E^E^ + pj cos^ 6*) + C^m^m^ + 2C^Ep^ cosO*] 



487r2 P E 



(A.30) 

Although the notation is different, this expression is found to agree with the equivalent 
result in [196]. 

A. 5 Full cross section 

To integrate this we use dQ* = 2TTd{cos6*), and so obtain 



<7(qq L-L+) 



4 1 

e I P3 
Un P E 



Ci ( EsE, + C^m.m, 



(A31) 
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A. 6 Slepton Drell-Yan cross section 



A. 6 Slepton Drell-Yan cross section 

The slepton calculation is in many ways similar to the lepton one above, although it is 
made slightly simpler by the scalar coupling of the sleptons to Z and 7. Left-handed and 
right-handed scalar states can mix but we just consider the cross section for production 
of one or the other (L below). 

For Z the vertex factor is then i^ghiPs — Pa)^, and hence 



-iMi = 



which gives 



1 • • 

v{P2)- - Wil^i^ (cV - Ca7^) u{pi) 



s — rriy + imzTz 



my(P3 -P4) 



(A.32) 



9l 



S — 777.2 + if^Z^Z 



v{P2)y\ (cV - 47^) u{pi) 



91./ \v 
y(P3 -P4) 



(A.33) 



Note that the first square bracket is exactly as previously, so contributes the same factor 
when we evaluate X^gpins l-^fi The second square bracket contributes 



9l 



iP3 - P^rips - pd" ■ 



(A.34) 



Multiplying these together we obtain 



9i 



4 (s - miy + 777|r| 



X 



2 

{cyY + (ca)^ ^ [2Pl-(P3 -P4)P2-(P3 -P4) -Pl-P2(P3 -P4)-(P3 -P4)])- 

(A.35) 
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This time to obtain the 7-exchange version it is necessary to make the additional 
substitution of = 2. Therefore 



4 s2 



4 [2pi.(p3 - P4)P2.(P3 - P4) - P1-P2(P3 - P4).(P3 " P4) ] ■ 

(A.36) 



The final term to be calculated is once again the interference term. Neglecting the 
propagators, the spin sum is 



Tr 



7-Pi (cV + ck7') |^7-P27.] f (P3 - P4)^(P3 - P^y. (A.37) 



Multiplying these terms (once again noting that the trace is the same as appeared 
in the lepton cross section) we obtain 



1 gle'^qW 



4 s{s — m| — imzFz) 



X CygL [2pi.(p3 - P4)P2-(P3 - P4) - Pl-P2(P3 " P4)-(P3 " P4) ] • (A. 38) 



Making the additional change of variables hj^ = gi^gz/^e, we arrive at the total matrix 
element squared: 



|Mfi|^) = —D[2pi.{p3 -Pi)p2-{P3 -Pi) - Pl-P2iP3 - Pi)-{.P3 - P4)] , 



(A.39) 



where 



hi 



D 



(s — m|) + mlr^ 



(s-m|) +m|r| 
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A. 6 Slepton Drell-Yan cross section 



Therefore the differential cross section is 
dcT 1 

dn* " 647r2s3 [2P1.(P3 -P4)P2-(P3 -P4) - Pl-P2{P3 - P4)-{P3 ~ Pi)] ■ (A.41) 

The 4-momenta in ^A.2I give us 

Pi-(P3 - Pi)p2.{p3 - Pa) =E^ [ {Es - E^f - Apl cos^ 9*] , 

(A.42) 

PiMps - Pa)-{P3 - Pa) =2E^ [ {E^ - E^f - 4^^] . 

Therefore 

2pi-(P3 - Pa)P2-{P3 - Pa) - Pi-P2{P3 - Pa)-{P3 - Pa) = 8E^pl sin^ 9* , (A.43) 
and the differential cross section is 



dn* 327r2 s2 E 



D-^pism'9*. (A44) 



The final result for either left- or right-handed sleptons, including the 1/3 colour factor, 
is therefore 



Appendix B 

Black Hole Master Equation 



In this appendix, Kanti's derivation of a master equation is reproduced from [2]. The 
Newman-Penrose formahsm is used to obtain this equation describing the motion of a 
particle with spin s in the background of a higher-dimensional, non-rotating, neutral 
black hole projected onto a 3-brane. The corresponding 4-dimensional metric tensor is 
given in eq. ()4.4p . 

We first need to choose a tetrad basis of null vectors {i'^, n'^, m^, m'^), where i and n 
are real vectors and m and m are a pair of complex conjugate vectors. They satisfy the 
relations 1 • n = 1 and m • m = — 1, with all other products being zero. Such a tetrad 
basis is given by 



r=(i 1,0.0). „"=(i.-|,o,o 



m'-IO, 0. 1. -^)-=-, m" = (O, 0, 1, ^ . (B.l) 



The \ahc coefficients, which are used to construct the spin coefficients, are defined by 



where stands for each one of the null vectors and 
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(i, j) denote the components of each vector. Their non- vanishing components are found 
to be 

\ - ^' \ - ^ \ - ^ \ - ^"^^ m 9^ 

^^122 — — TT 1 ^^134 — - , ^234 — —-pT i ^^334 — "7= : T ■ K^-'^) 

2 r 2r V2rsm6' 

The above components must be supplemented by those that follow from the symmetry 
Aa&c = —^cba and the complex conjugates obtained by replacing an index 3 by 4 (or vice 
versa) or by interchanging 3 and 4. 

We may now compute the spin coefficients defined by ^ahc — {^abc + ^cah — ^bcajf^. 
Particular components, or combinations, of the spin coefficients can be directly used in 
the field equations [173, 174]. They are found to have the following values: 

K — a — X — u — T — TT — e — 0; 



1 h h' ^ cote 



In what follows, we will also employ the Newman-Penrose operators 



Id d I Id h d s 1/(9 i d \ 



^ hdt^dr' ^ 2dt 2dr' ^ y^r \d9 ^ sinO dip) ' ^^'^^ 
and make use of the field factorization 

^s{t, r, 9, if) = e-'^' e*'"'^ R,{r) sSr{9) , (B.5) 

where sY^'"(^, (f) ~ e*"**' sS'^{9) are the spin- weighted spherical harmonics [175]. We will 
now consider each type of field separately. 
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B.l Gauge bosons (s = 1) 

In the Newman-Penrose formalism, there are only three 'degrees of freedom' for a gauge 
field (namely $0 = -^is? ^1 = (-^12 + -^43)72 and $2 = -F42) in terms of which the different 
components of the Yang- Mills equation for a massless gauge field are written as 



{D - 2p) $1 - {S* - 2a) $0 =0 , 
<5<l>i - (A + /i-27) $0 =0, 
{D - p) $2 - 5* $1 =0 , 
{5 + 2(3) $2 - (A + 2/i) $1 =0 , 



(B.6) 
(B.7) 
(B.8) 
(B.9) 



where 5* stands for the complex conjugate of 5. Rearranging eqs. ()B.6|) and ()B.7|) . 
we can see that $1 decouples leaving behind an equation involving only $o- Using the 
explicit forms of the operators and spin coefficients, as well as the factorized ansatz of 
eq. ()B.5|1 . this can be separated into an angular equation. 



1 d f . d^S^f 
— sm ^ 



sin 6' dO 



dO 



+ 



2m cot 9 
sin 6' 



m 



+ 1 - cot^6' + lA^ 



,ST{e) = Q, (B.IO) 



with eigenvalue ^A^ = £ (£ + 1) — s (s + 1), and a radial equation. 



1 d 
A dr 



^ 2 dRi 
dr 



+ 



h 



iuor'^h' , ^ ,, 
+ 2iujr — + (A" - 2) - iA| 



RAr) = 



(B.ll) 



where A = hr"^. 
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B.2 Fermion fields (s = 1/2) 



B.2 Fermion fields (s = 1/2) 



For a massless two-component spinor field, the Dirac equation can be written as 



{5* - a) xo -p)xi, 
(A + /^-7)Xo ={^ + (3) Xi- 



(B.12) 
(B.13) 



Performing a similar rearrangement as in the case of bosons, we find that xi is decoupled 

and that the equation for xo reduces to the following set of angular. 



— I sm 6* — 4^ I + 



sin^ de 



and radial. 



de 



moot 9 1 1 



sin^ sin^e 2 4 



+ - - - cot' 9 + i/2A^ 



(B.14) 



A dr\ dr 



+ 



h 



+ ^^r -^ + - (A" -2) -1/2 A, 



i?i/2(r) = 0, (B.15) 



equations, with the same definitions for A and ^A^ as before. 



B.3 Scalar fields (s = 0) 

For completeness the equation of motion for a scalar field propagating in the same 
background is included here. This equation can be determined quite easily by evaluating 
the double covariant derivative g^^^D^Di, acting on the scalar field. It finally leads to 
this pair of equations: 



J_d_( . doSr 
sm9d9V'''^ d9 



+ 



m 



sin' 9 



qSJ" — ; 



(B.16) 
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Here, oYP'[9,ip) = e*™*^ 0'S'™(^) are the usual spherical harmonics Yl"{6,Lp) and qA^ = 
£{£ +1). The above equations were used in [170] for the analytic determination of 
grey-body factors for the brane emission of scalar particles by higher-dimensional black 
holes. 



B.4 Master equation for a field with arbitrary spin 

Combining the equations derived above for bosons, fermions and scalar fields, we may 
now rewrite them in the form of a master equation, valid for field types. The radial 
equation then takes the form 



A-^ ^ ( A 
dr 



s+i dRg 
dr 



+ 



h 



j CI.) nn^ h ' 

^2tsu;r- ^ + s (A" - 2) - , ) i?,(r) = , (B.18) 



while the angular equation reads 



1 d 



sm9 d9 



-r^ sine 



dS 



si 



de 



2ms cot 6 



+ 8-8^ cot^ .A, 



sin 9 sin 9 



^"^, = 0. (B.19) 



The latter equation is identical to the one derived by Teukolsky [172] in the case 
of a non-rotating, spherically-symmetric black hole. However the radial one differs by 
the extra factor s (A" — 2) because for the metric tensor here this combination is not 
zero (unlike for the 4-dimensional Schwarzschild and Kerr metrics). The A"-term can 
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B.4 Master equation for a field with cLrbitrEiry spin 



be removed if we make the redefinition Rg — A ^ Pg. Then, we obtain 

where now ^A^ — sXe + — i {i + 1) — s {s — 1). The above form of the radial equation 
was used in [171] to determine analytically the grey-body factors and emission rates for 
fermions and gauge bosons on the brane. 
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C.l Monte Carlo calculation of the cross section 

The main features of the Monte Carlo (MC) method can be described by the following 
mathematical relationship: 

"'^1 i=i 
where the f{xi) are values of f{x) at the N randomly chosen values of x e [xi, X2]. 
In this case the integral we require is for the black hole production cross section a: 

(7= / a{E)p{E)dE, (C.2) 

where E is the centre-of-mass energy of the partons involved and a is the parton-level 
cross section. Using the substitution y{E) — dY/dE we obtain 

.^r~^(^dr. (C.3) 

Jy^=y{e^^) y{E) 
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C.l Monte Carlo calculation of the cross section 



The MC procedure will be most efficient if the integrand is as flat as possible as 
a function of E. The parton-level geometrical cross section for black hole production 
means that o"(i?) ~ where (3 = 2/{n + 1). In addition we expect that p{E) ~ E~^ 
because of the behaviour of the parton distribution functions (PDFs). Hence a sensible 
choice for y{E) is 

y{E) = aE'^-^ Y{E) = E", (C.4) 

where a = P — 7. Returning to eq. (jdlll we have 

,J criE)piE) \ E^,^- E^J a{E)p{E) \ 
a = iY{E^.^) - [-^^) = \-^^) ' (C-5) 

where Y{E) is chosen at random from a uniform distribution between the maximum and 
minimum values. 

p{E) is calculated from the PDFs and is more commonly written as a function of 
T = s/s (hence r oc E"^). Therefore p{E) can be re-written as 

2r dr 

/ ax ...... . 



/ 2r , , 2r /" dx ^ , . ^ /t 



This integral can also be evaluated by a MC procedure, and once again a change of 
variable is useful. We may conveniently write p(r) as 

P{r) = - ['-xMx)-h (-) = - ['-h{x,T), (C.7) 

T X X \X/ T X 

and the substitution z{x) = dZ/dx gives us that 

l^M^^^^ (C.8) 

- 'z(r) XZ{X) 



r 
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The assumption that h{x, r) is fairly smooth (since f{x) is usually of the approximate 
form l/x) means that the MC procedure should be efficient if 

z(x) = - ^ Z(x) = \jix. (C.9) 

X 

So, using eq. ()C.1|) . we conclude that 

P{r) = ^^^'^-^^^^^ {h{x, r)) = -i^ {h{x, r)) , (C.IO) 

T T 

where Z{x) is chosen at random from a uniform distribution between the maximum and 
minimum values. 

Combining these two MC procedures, a good approximation to the total cross section 
a is given by 



^^^^^^ , (C.ll) 

where firstly Y{E) is chosen at random (— ^ E and hence r) and then, given this value 
of r, Z{x) is chosen at random (— > x). 



C.2 Monte Carlo generation of unweighted events 

As well as estimating the total cross section, the MC program must generate events 
drawn from this distribution. This can be done using essentially the same procedure — 
E and then x are chosen and the cross section value is computed. The calculated weight 
is then accepted or rejected against the maximum weight found by an initial search. A 
weight Wi is accepted with probability Wi/wumx- 
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